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Abstract 

Inverse scattering and spectral one-dimensional problems are discussed systematically in a self- 
contained way. Many novel results due to the author are presented. The classical results are often 
presented in a new way. Several highlights of the new results include: 

1) Analysis of the invertibility of the steps in the Gel'fand-Levitan and Marchenko inversion procedures, 

2) Theory of the inverse problem with /-function as the data and its applications; 

3) Proof of the property C for ordinary differential operators, numerous applications of property C; 

4) Inverse problems with "incomplete" data; 

5) Spherically symmetric inverse scattering problem with fixed-energy data: analysis of the Newton- 

Sabatier (NS) scheme for inversion of fixed-energy phase shifts is given. This analysis shows that 
the NS scheme is fundamentally wrong, and is not a valid inversion method. 

6) Complete presentation of the Krein inverse scattering theory is given. Consistency of this theory is 

proved. 

7) Quarkonium systems; 

8) A study of the properties of /-function; 

9) Some new inverse problems for the heat and wave equations are studied. 

10) A study of inverse scattering problem for an inhomogeneous Schrodinger equation; 
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Chapter 1 

Introduction 



1.1 Why is this paper written? 

There are excellent books ^M, and |Lj, where inverse spectral and scattering problems are discussed in 
detail. The author decided to write this paper for the following reasons: 1) He gives a new approach 
to the uniqueness of the solutions to these problems. This approach is based on property C for Sturm- 
Liouville operators; 2) the inverse problem with /-function as the data is studied and applied to many 
inverse problems; 3) a detailed analysis of the invertibility of the steps in Marchenko and Gel'fand- 
Levitan (GL) inversion procedures is given; 4) inverse problems with "incomplete" data are studied; 5) 
a detailed presentation of Krein's inversion method with proofs is given apparently for the first time; 
6) a number of new results for various inverse problems are presented. These include, in particular, a) 
analysis of the Newton-Sabatier (NS) inversion scheme for finding a potential given the corresponding 
fixed-energy phase shifts: it is proved that the NS scheme is fundamentally wrong and is not an inversion 
method; b) a method for finding confining potential (a quarkonium system) from a few experimental 
data; c) solution of several new inverse problems for the heat- and wave equations; d) a uniqueness 
theorem for finding a potential q from a part of the corresponding fixed-energy phase shifts; and many 
other results which are taken from [^, |R1| - |R29| . 

Due to the space limitations, several important questions are not discussed: inverse scattering on the 
full line, iterative methods for finding potential q: a) from two spectra [H ] . |R5| . b) from S'— matrix alone 
when q is compactly supported pR9) . approximate methods for finding q from fixed-energy phase shifts 
| R14| . [RT5] . property of resonances 0, jR29| . inverse scattering for systems of equations, etc. 

1.2 AuxiUary results 

Let q{x) e Li,i, Li,„ = {q : q{x) = q{x), {l + x)'^\q{x)\dx < oo, and q G Lf^^{R+)}, where Lf^^ (R+) 
consists of functions belonging to L^(0,a) for any a < oo, and overline stands for complex conjugate. 

Consider the differential expression £u — —u" + q{x)u with domain of definition D{lo) = {u : u(0) = 
0, u G Cq (0, oo)}, where Cq (0, oo) is the set of C^(R+)-functions vanishing in a neighborhood of infinity, 
K+ := [0, oo). U H is the Hilbert space L^(R+), then £o is densely defined symmetric linear operator 
in H, essentially self-adjoint, that is, the closure £ of £o in H is selfadjoint. It is possible to construct a 
selfadjoint operator £ without assuming that q S Lfg^{R+). Such a theory is technically more difficult, 
because it is not even obvious a priori that the set D{£q)u := {u : u G Cq(M+), £u G L^(M_|_)} is dense 
in H (in fact, it is dense). Such a theory is presented in jNai| . If one drops the assumption q G Lf^c^ 
then D{£q) is not a domain of definition of £ since there are functions u G 'D{£q) for which £u ^ _L^(R+). 
In the future we mean hy £ a self-adjoint operator generated by the differential expression £ and the 
boundary condition u{0) = 0. 
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This operator has absolutely continuous spectrum, which fills (0,oo), and discrete, finite, negative 
spectrum {~fc|}i<j<j, where — fcj are the eigenvalues of i, all of them are simple, 

% + =-fc>„ ^,(0)=0, ^;(0) = 1, (1.2.1) 

where (pj are corresponding eigenfunctions which are real-valued functions, and 

1 



c. 



ip^^dx. (1.2.2) 

The functions </3(x, k) and 0{x, k) are defined as the unique solutions to the problems: 

^ip = k^ip, x>0; p{0, k) = 0, (^'(0, k) = 1, (1.2.3) 



ie^k^e, x>0; 9{0,k)^l, e'{0,k) = 0. (1.2.4) 

These functions are well defined for any q{x) G Ll^^{S.^). Their existence and uniqueness can be proved 
by using the Volterra equations for if and 0. If q G then the Jost solution f{x,k) exists and is 

unique. This solution is defined by the problem: 

if := -/" + qf = k^f, fix, k) = cxp(ifca;) + o(l) as 2: ^ +oo- /(O, k) := f{k). (1.2.5) 

Existence and uniqueness of / is proved by means of the Volterra equation: 

f{x,k)^exp{ikx)+ r ^^"[^(^~^)] q(t)f{t,k)dt. (1.2.6) 

If g G ii.i then this equation implies that /(x, k) is an analytic function of k in C+ = {k : Imk > 0}, 
f{x,k) — f{x,—k) for fc > 0. The Jost function is defined as f(k) := f{0,k). It has exactly J simple 
roots ikj, kj > 0, where —k'j, 1 < j < J, are the negative eigenvalues of £. The number fc = 
can be a zero of /(fc). If /(O) = 0, then /(O) 7^ 0, where /(fc) := Existence of /(O) is a fine 
result under the only assumption q G £1,1 ( see Theorem 3.1.3 below, and ^) and an easy one if 
9 € ii,2 '■= {q '■ q = /o°°(l + x'^)\q{x)\dx < 00}. The phase shift S{k) is defined by the formula 

fik) = \fik)\exp{-t5{k)), <5(oo) = 0, /(^) = 1, (1.2.7) 



where the last equation in H1.2.18(l follows from p.2.6|l . Because q{x) = q{x), one has 5{—k) = —S{k) 
for fc G K. One defines the S'-matrix by the formula 

S{k):^^jr^, fcGM. (1.2.8) 

The function S{k) is not defined for complex fc if q G ii.i, but if \q{x)\ < ciexp{—C2\x\'^), 7 > 1, then 

/(fc) is an entire function of fc and S{k) is meromorphic in C. If q{x) — for x > a, then /(fc) is an 
entire function of exponential type < 2a (see Section l5.1|l . 

If g G il l, then at fc^ = —k'^,kj > 0, the Jost solution fj{x) := f{x,ikj) is proportional to 
'Pjix) :— if{x,ikj), fj and (pj both belong to L^(R_|_). The integral equation for (p is: 

, ,, sin(fcx) sin k(x — s) ,,,,,, „ „n 

^{x, fc) = — ^ + ^ -q{s)p{s, k)ds. (1.2.9) 

One has: 

, f{x,k)f{-k)- f{x~k)f{k) 

ip{x,k)^ ^ , (1.2.10) 
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because the right-hand side of p.2.1UI) solves equation (|1.2.5|l and satisfies conditions (|1.2.3|l at x = 0. 
The first condition 11.2.31) is obvious, and the second one follows from the Wronskian formula: 



/'(O, k)f{-k) ~ /'(O, -fc)/(fc) = 2ik. (1.2.11) 

If A; = ikj then /j(x) S i^(R+), as one can derive easily from equation (|1.2.6|l . In fact, |/j(a:)| < ce~''^^, 
X > 0. If fc > 0, then f{x,—k) — f{x,k). li q = q then fj(x) is a real-valued function. The function 
f{x,k) is analytic in C+ but is, in general, not defined for k G C_ := {k : Imk < 0}. In particular, 
H1.2.11|l . in general, is valid on the real axis only. However, if \q{x)\ < ciexp{—C2\x\'^), 7 > 1, then /(fc) 
is defined on the whole complex plane of fc, as was mentioned above. Let us denote /(x, k) := f+{x, k) 
for k E C+ and let f-{x, k) be the second, linearly independent, solution to equation (|1.2.5|l for k e C+. 
If /+ e i^(M+), then /_ ^ L^(K+). One can write a formula, similar to H1.2.1()|l . for fc e C+: 

^{x, k) = c(fc)[/_(0, k)f{x, k) ~ /(O, k)f^[x, fc)], (1.2.12) 

where c{k) — const 7^ 0. For Lp{x,ikj) e L^(M_|_), it is necessary and sufficient that f{ikj) = 0. In fact 

f{ik,)=0, f{ik,)^0, l<j<J, (1.2.13) 

^ f - ^ 

gets 



where f — j^- To prove the second relation in (|1.2.13l) . one differentiates (|1.2.5|) with respect to k and 



f" + k^f^qf^-2kf. (1.2.14) 

Existence of the derivative / with respect to k in C+ follows easily from equation H1.2.6(l . Multiply 
H1.2.14|l by / and l|1.2.5|l by /, subtract and integrate over R+, then by parts, put k = ikj, and get: 

/"OO 

-2ik, / ffdx = (//' - /'/)|(r = /'(O, ^k,)n^k,). 
Jo 

Thus 

f.,,^mMm..^l^o. (1.2.15) 



s 



It follows from H1.2.15|l that f{ikj) 7^ 0. The numbers sj > are called the norming constants: 



/'(o,^fc,)/(^fc,; 

Definition 1.2.1. Scattering data is the triple: 



l<j<J. (1.2.16) 



S := {S{k), k„ s„ 1 < J < J}, S{k) := k, > 0, s, > 0. (1.2.17) 

The Jost function f{k) may vanish at fc = 0. If /(O) = 0, then the point fc = is called a resonance. 
If \q{x)\ < ciexp{—C2\x\'^), 7 > 1, then the zeros of /(fc) in C_ are called resonances. As we have seen 
above, there are finitely many zeros of /(fc) in C+, these zeros are simple, their number J is the number 
of negative eigenvalues —kj, I < j < J, of the selfadjoint Dirichlet operator £. If q £ Li^i then the 
negative spectrum of £ is finite 'W . 

The phase shift (5(fc), defined in (|1.2.18|l . is related to 5(fc): 

5(fc) = e^**^'^'), (1.2.18) 
so that S'(fc) and 6{k) are interchangeable in the scattering data. One has fj{x) — f'{0,ikj)Lpj{x), 

f'{0,tk,) 



because {fe'^-t-'.'l solves (|1.2.3f) . Therefore 



1 1 
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Thus 



2ikjf'{0,ikj) 



1<J<J- 



(1.2.20) 



In Section l¥!n the notion of spectral function p{X) is defined. It wiU be proved in Section I^TI for q E Li.i 
that the formula for the spectral function is: 

VXdX ^ ^ ^ 

^I/(Va)P' - ' 

dp{X) = { J (1.2.21) 

CjSiX + k^)dX, A < 0, 

where Cj are defined in (|1.2.19|I - I|1.2.2Q|I . The spectral function is defined in Section [4.11 for any q E 

L;^Q^(M+), q = q. Such a q may grow at infinity. On the other hand, the scattering theory is constructed 
for q E 

Let us define the index of S{k): 



J := ind S{k) 



-i^ARarg^(fc) = ^ 



dlnS'(fc). 



This definition implies that indS'(fc) = ind/(— fc) — ind/(fc) = — 2ind/(fc). Therefore: 

-2Jif/(0)^0, 



ind S{k) = 



-2 J - 1 if /(O) = 0, 



(1.2.22) 



(1.2.23) 



because a simple zero fc = contributes ^ to the index, and the index of an analytic in C-|_ function 
/(fc), such that /(oo) = 1, equals to the number of zeros of /(fc) in C+ plus half of the number of its 
zeros on the real axis, provided that all the zeros are simple. This follows from the argument principle. 

In Section 14.21 and Section 15.21 the existence and uniqueness of the transformation (transmutation) 
operators will be proved. Namely, 

sin(fca;) 



ip{x,k) ^ ipo{x,k) + K{x,y)ipo{y)dy:= {I + K)ipQ, (po 



and 



fix, fc) = e*'=- + / A{x, y)e'^ydy := (/ + A)h, /„ := e 



ikx 



(1.2.24) 



(1.2.25) 



and the properties of the kernels A{x,y) and K{x,y) are discussed in Section [5.21 and Section [4.21 
correspondingly. The transformation operator I + K transforms the solution ipo to the equation (|1.2.3() 
with q = into the solution ip of (|1.2.3|l . satisfying the same as ipo boundary conditions at a; = 0. The 
transformation operator I + A transforms the solution /o to equation H1.2.5|) with q ^ into the solution 
/ of H1.2.5|) satisfying the same as /o "boundary conditions at infinity" . 
One can prove (see [M] and Sec. 5.7) the following estimates 



\A{x,y)\<ca(^^^y c = const > 0, '^i^) ■= k{t)\dt, 



x + y 



1 (x + y 
79 



< ca{x)a 



x + y 



and A[x, y) solves the equation: 



qds 



ds 



dtq{s - t)A{s - t, s + t). 



(1.2.26) 
(1.2.27) 

(1.2.28) 



By H™' = if'"(R-)-) we denote Sobolov spaces M^™-^. The kernel A{x, y) is the unique solution to l|1.2.28|l , 
and also of the problem H5.1.1(l - H5.1.3|l . 
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1.3 Statement of the inverse scattering and inverse spectral 
problems. 

ISP: Inverse Scattering problem (ISP) consists of finding q G Li.i from the corresponding scattering 
data S (see (jTTBll ). 

A study of ISP consists of the foUowing: 

1) One proves that ISP has at most one solution (see Theorem 5.2.1). 

2) One finds necessary and sufficient conditions for S to be scattering data corresponding to a q E Li^i 

(characterization of the scattering data problem) . 

3) One gives a reconstruction method for calculating q e Li i from the corresponding S. 
In Chapter |31 these three problems are solved. 

ISpP: Inverse spectral problem consists of finding q from the corresponding spectral function. 
A study of ISpP consists of the similar steps: 

1) One proves that ISpP has at most one solution in an appropriate class of q: if qi and q2 from this 

class generate the same p{X), then qi = q2- 

2) One finds necessary and sufficient conditions on p(A) which guarantee that p(A) is a spectral function 

corresponding to some q from the above class. 

3) One gives a reconstruction method for finding q{x) from the corresponding p{X). 



Denote by im operators £ corresponding to potentials g„i G ii.i, and by fmix, k) the corresponding Jost 
solutions, m — 1,2. 

Definition 1.4.1. We say that a pair {ii,£2} has property C+ iff the set {/i(a;, A:)/2(a;, fc)}vfc>o is 
complete (total) in L^{M.-f-). 

This means that ii h E i^(R+) then 



We prove in Section ITU that a pair does have property C+ if qm E ^14. Let £ip := —ip" + 

q{x)tp, and let ipj correspond to g = g^, 

£(^-fcV = 0, v9(0,/c) = 0; (^'(0,fc) = l; £9 - k^O = 0, e{0,k)^l, e'{0,k)^0. (1.4.2) 

Definition 1.4.2. We say that a pair {£1, £2} has property C-^- iff the set {(pi{-,k)ip2{-,k)} is complete 
in L^{0, b) for any & > 0, b < 00. 

This means that \i h E i^(0, b), then: 



1.4 Property C for ODE. 




h{x)fi{x,k)f2{x,k)dx = Vk>0} ^ h^O. 



(1.4.1) 




(1.4.3) 



In Theorem 12 . 2 . 21 we prove that there is a ft, 7^ for which 




Vfc > 
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for a suitable qi q2, ?i, 92 G ii.i- Therefore Property with b = 00 does not hold, in general. 

Property Cg is defined similarly to Property C<^, with functions 0j{x,k) replacing (pj(x,k). 

In Chapter 121 we prove that properties C+, and Ce hold, and give many applications of these 
properties throughout this work. 

1.5 A brief description of the basic results. 

The basic results of this work include: 

1) Proof of properties C+, and Cg. Demonstration of many applications of these properties. 

2) Analysis of the invcrtibility of the steps in the inversion procedures of Gel'fand-Levitan (GL) for 

solving inverse spectral problem: 

p^L^K^q, (1.5.1) 

where 

dK(x,x) 

, = 24-i, (1.5.2) 

the kernel L — L(x, y) is: 

^oix,X)ipoiy,X)da{X), da{X) := d[p{X) - poiX)], (1.5.3) 

-00 

, I , A>0, 

dpo = < TT 

[ 0, A < 0. 

Pq = po is the spectral function of £ with q = 0, and K solves the Gel'fand-Levitan equation 

K{x,y)+ ( K{x,s)L{s,y)ds + L{x,y) ^0, 0<y<x. (1.5.4) 
Jo 



Our basic result is a proof of the invcrtibility of all the steps in l|1.5.1(l : 

p^L^K^q, (1.5.5) 
which holds under a weak assumption on p. Namely, assume that 

peg, (1.5.6) 

where G is the set of nondecreasing functions p of bounded variation on every interval (—00, 6), b < 00, 
such that the following two assumptions, Ai) and A2) hold. 

Denote Lq(M_|_) the set of L^(R+) functions vanishing in a neighborhood of infinity. Let h E Lq(M+) 
and H{X) :— h{x)ipo{x, X)dx. 

Assumption Ai) is: 

/oo 

H'^{X)dp{X) ^ 0, then h = 0. (1.5.7) 
-00 

Let 

poo 

H:={H{X) ■.heC^{R+)}, H{X):= h{x)ipo{x, X)dx, (1.5.8) 
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Pi and p2 belong to V, and v := pi — p2 (see Section 4.2). 
Assumption A2) is: 

/CO 
H^{X)dv = {) VH eH, then = Q. (1.5.9) 
'CO 

In order to insure the one-to-one correspondence between spectral functions p and selfadjoint op- 
erators £, we assume that q is such that the corresponding £ is "in the limit point at infinity case" . 
This means that the equation (£ — z)u — 0, Imz > has exactly one nontrivial solution in i^(R+), 
£u = —u" ^ q{x)u. If q S Li^i then £ is "in the limit point at infinity case". 

3) Analysis of the invcrtibility of the Marchcnko inversion procedure for solving ISP: 

S ^ F ^ q, (1.5.10) 

where 

1 ^ 
Fi^)-=7r [l~S{k)y''^dx + ys,e-''^-:=Fs{x)+Fd{x), (1.5.11) 
27r J -00 f-{ 

, , dA(x,x) 
g(a;) = -2— i-^, (1.5.12) 

and y) solves the Marchenko equation 

A{x,y)+ A{x,s)F{s + y)ds + F{x + y) ^0, 0<x<y<oo. (1.5.13) 

J X 

Our basic result is a proof of the invcrtibility of the steps in p. 5. 10(1 : 

S ^ F ^ A^ q (1.5.14) 
under the assumption q e We also derive a new equation for 



AiO.y), y>0, 
0, y<0. 



A{y) := 
This equation is: 

F{y) + A{y)+ A{s)F{s + y)ds = A{-y), -00 < y < 00. (1.5.15) 
Jq 

The function A{y) is of interest because 

poo 

f{k) = l+ A{y)eJ''ydy -.^l + Aik). (1.5.16) 



Therefore the knowledge of A{y) is equivalent to the knowledge of f{k). 

In Section r5.5l we give necessary and sufhcient conditions for S to be the scattering data corresponding 
to q € Li^i. We also prove that if 

|g(x)| <ciexp(-C2|xr), 7 > 1, (1.5.17) 

and, in particular, if 

g(x) = for a: > a, (1.5.18) 

then S{k) alone determines q{x) uniquely, because it determines fcj, sj and J uniquely under the as- 
sumption (|1.5.17|l or p. 5. 18(1 . 
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4) We give a very short and simple proof of the uniqueness theorem which says that the /-function, 



Hk):^^^^, Vfc>0, (1.5.19) 

determines q € Li.i uniquely. The /-function is equal to Weyl's m-function if (/ G Li.i. 

We give many applications of the above uniqueness theorem. In particular, we give short and simple 
proofs of the uniqueness theorems of Marchenko which say that S determines q G Li i uniquely, and p{X) 
determines q uniquely. We prove that if (|1.5.18|l (or Hl.S.lTfl l holds, then either of the four functions 
S{k), S{k), f{k), /'(O, fc), determines q{x) uniquely. This result is applied in Chapter [TUl to the heat and 
wave equations. It allows one to study some new inverse problems. For example, let 

utt = Uxx - q{x)u, X > 0, i > 0, (1.5.20) 

u^Ut^O &tt = Q. (1.5.21) 

w(0, t) = 5{t) or u'(0, t) = 5{t). (1.5.22) 

Assume 

q==Ofora;>l, q^q, q e L\0,1), (1.5.23) 
and let the extra data (measured data) be 

u{l,t)^a{t) yt>0. (1.5.24) 

The inverse problem is: given these data, findq{x). 
Another example: Let 

ut = Uxx - q{x)u, 0<a;<l, <> 0, q e L^[0,l], (1.5.25) 



u(a;,0)=0 (1.5.26) 
M(0,i) = 0, u(l,t)=a(t), a{t)eL^{R+), a ^ 0. (1.5.27) 

The extra data are 

Uxil,t) = b{t) yt>0. (1.5.28) 
The inverse problem is: given these data, findq{x). 

Using the above uniqueness results, we prove that these two inverse problems have at most one 
solution. The proof gives also a constructive procedure for finding q. 

5) We have already mentioned uniqueness theorems for some inverse problems with "incomplete data" . 
"Incomplete data" means the data which are a proper subset of the classical data, but "incomplete- 
ness" of the data is compensated by the additional assumptions on q. For example, the classical 
scattering data are the triple (|1.2.17|l . but if (|1.5.18|l or H1.5.17|l is assumed, then the "incomplete 
data" alone, such as S{k), or 6{k), or /(fc), or /'(O, fc), Vfc > 0, determine q uniquely. Another 
general result of this nature, that we prove in Chapter is the following one. 
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Consider, for example, the problem 

iifj^Xjipj, 0<a;<l; (pj{0) ^ ipj{l) = 0. (1.5.29) 

Other boundary conditions can also be considered. 
Assume that the following data are given. 

{A™o)Vj; 5 < x < 1}, q{x) e L'[0,1], q ^ q, (1.5.30) 

where < 6 < 1, and 

j 

m{j) = —(1 + €j), \ej\ < 1, ^ as J — !■ oo, a = const, < cr < 2. (1.5.31) 

Assume also 

C30 

^|e,|<oo. (1.5.32) 

We prove 

Theorem 1.5.1. Data (|1.5.30|l - (|1.5.31|l determine uniquely q{x) on the interval 0<x<bifa> 2b. If 
H1.5.32|l is assumed additionally, then q is uniquely determined if cr > 2b. 

The a gives the "part of the spectra" sufficient for the unique recovery of q on [0, b]. For example, if 
b — ^ and (|1.5.32() holds, then ct = 1, so "one spectrum" determines uniquely q on [0, If 6 = i, then 
a = i, so "half of the spectrum" determines uniquely q on [0, -j]. If 6 = i, then "| of the spectrum" 
determine uniquely q on [0, If 5 = 1, then cr = 2, and "two spectra" determine q uniquely on the 
whole interval [0,1]. The last result belongs to Borg [B]. By "two spectra" one means {Aj} 
where /ij are the eigenvalues of the problem: 

(uj^fijUj, Uj(0)=0, u'^{l) + huj{l)^0. (1.5.33) 

In fact, two spectra determine not only q but the boundary conditions as well [M| . 

6) Our basic results on the spherically symmetric inverse scattering problem with fixed-energy data are 

the following. 

The first result: If g = q{r) = for r > a, a > is an arbitrary large fixed number, r := \x\, x M.^ , 
q = q, and r'^\q{r)\'^dr < oo, then the data {Se}\/i^c determine q{r) uniquely. Here Si is the phase 
shift at a fixed energy fc^ > 0, £ is the angular momentum, and C is any fixed set of positive integers 
such that 

eec 

The second result is: If g = q{x), x G M^, q ~ for |a;| > a, q G L'^{Ba), where Ba '.= {x : |a;| < a}, 
then the knowledge of the scattering amplitude A{a' ,a) at a fixed energy fc^ > and all a' £ Sj 
determine q{x) uniquely |3, |R7| . Here Sj, j — 1,2, are arbitrary small open subsets in S"^ and is 
the unit sphere in R^. The scattering amplitude is defined in Section imi 

The third result is: The Newton-Sabatier inversion procedure (see |US| . El) fundamentally wrong. 

7) Following |R16| we present, apparently for the first time, a detailed exposition (with proofs) of the 

Krein inversion theory for solving inverse scattering problem and prove the consistency of this 
theory. 

8) We give a method for recovery of a quarkonium system (a confining potential) from a few experimental 

measurements. 

9) We study various properties of the /-function. 

10) We study an inverse scattering problem for inhomogeneous Schrodinger equation. 
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Chapter 2 

Property C for ODE 



2.1 Property 



By ODE in this section, the equation 

{£ - k^)u := -u" + q{x)u - k^u = (2.1.1) 

is meant. Assume q G Li ^. Then the Jost solution f(x,k) is uniquely defined. In Section ll.4l Defini- 
tion II. 4.'T1 property C+ is explained. Let us prove 

Theorem 2.1.1. If q E j = 1,2 then property C+ holds. 

Proof. We use (|1.2.25|) and H1.2.2()|) . Denote A{x,y) := Ai{x,y) + A2{x, s). Let 



= 



dxh{x)fi(x, k)f2{x, k) 



dxh(x) 



OO /'OO 



2ikx 



A{x,y)e'^ydy 

+ / / dydzA^{x,y)A2{x,zY''^y^^'^ 
for some h € L-'^(]R+). Set y-\-z = s,y — z — a and get 



(2.1.2) 



where 



^i(x,?/)A2(a;,z)e*'=(^+^)dydz = / T{x,sy^'ds, 

Jlx 



T{x, s)^^ J ^ ^Ai (^x, j A2 {^x, ' ^'^^ 



(2.1.3) 



(2.1.4) 



Thus, /1/2 = (/+V"*)e2*'=^, where V* is the adjoint to a Volterra operator, V* f 2 /^"^^ A(a;, 2s)/(s)ds+ 

2XrT(a;,2s)/(s)ds. 

Using (|2.1.3|l and (|2.1.4|l one rewrites (|2.1.2|l as 







dse 



2iks 



h{s) + 2 / A{x,2s- x)h{x)dx + 2 / T{x,2s)h{x)d: 



,Vfc > 0. (2.1.5) 



The right-hand side is an analytic function of k in C+ vanishing for all fc > 0. Thus, it vanishes identically 
in C+ and, consequently, for fc < 0. Therefore 

h{s) + 2 [ A{x,2s-x)h{x)dx + 2 [ T(x, 2s)ft(a;)da; = 0, Vs > 0. (2.1.6) 
Jo Jo 
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Since A{x,y) and T{x,y) are bounded continuous functions, the Volterra equation (|2.1.6|l has only the 
trivial solution h — 0. □ 

Define functions g± and f± as the solutions to equation (|1.2.5|l with the following asymptotics: 

g± = exTp{±ikx) + o(l), x ~> —oo, (2.1.7) 

/± = exp(±zA;a;) + o(l), a; ^ +00, (2.1.8) 
Let us denote /+ = / and 5+ = g. 

Definition 2.1.2. The pair {li,l2\ has property C_ iff the set {(7i52}vfe>o is complete in _L^(M_). 

Similar definition can be given with {g-,j) replacing gj, j = I, 2. 
As above, one proves: 

Theorem 2.1.3. If qj e Li_i(M_), j = 1,2, then property C_ holds for 

By Li_i(M) we mean the set 

/oo 
{l + \x\)\q{x)\dx <(x}. (2.1.9) 
-00 

2.2 Properties Cip and Co. 

We prove only property C^. Property Ce is proved similarly. Property is defined in Section ll. 41 
Theorem 2.2.1. If qj E = 1,2, then property holds for {^1,^2}- 

Proof. Our proof is similar to the proof of Theorem 12. 1.11 Using H1.2.24|l and denoting (p — k(p, K :— 
Ki + if 2, one writes 

</'i02 = sin^(fca;) + / K{xTy)sa\(kx)sm{ky)dy 







X i^x 



Jo 



(2.2.1) 



Ki{x,y)K2{x, s){cos[k{y - s)] - cos[fc(j/ + s)]}dyds. 



Assume: 



Then 



0= / h{x)(j>i{x,k)(f>2{x,k)dx Vfc>0. (2.2.2) 







= / dxh{x) — i dxh{x) cos(2fca;) 
Jo JO 



+ / dscos(fcs) / dxh{x)K{x,x - s) (2.2.3) 







26 /'min(6,s) 

dscos{ks) / dxh{x)K{x, s — x) + I , 







where 



I ■= dxh{x) i i Ki{x,y)K2{x,s){cos[k{y — s)] — cos[k{y A- s)]}dyds 
Jo Jo JO 

Let y — s := t, y + s :— V. Then 

pX pX fX 

/ / KiK2C0s[k{y — s)]dyds = / dscos(fcs)-Bi(a;, s), 
^0 "'0 Jo 
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where 
Bi (x, s) 



2x-\s\ 



S + V 



2 

2x 



K2 X, 



V + s 



dv. 



f-X f-X t^^X 

I I K1K2 cos[k{y + s)]dyds — / B2{x, s) cos{ks)ds, 
Jo Jo Jo 



B2{x, s) 



Ki ( X, ^ \K2\x, 



s-t 



dt. 



and w = s if < s < x; w = 2.T — s if a; < s < 2x. 
Therefore 



2b 



dscos{ks) I dxh{x)Bi{x, s) — / dscos{ks) / dxh{x)B2{x, s) 



(2.2.4) 



From 12.2.3|l and l|2.2.4|l . taking k — > 00, one gets: 



h(x)dx — 0, 



and (using completeness of the system cos(fcs), < fc < 00, in L^(0, h)) the foUowing equation: 



+ / K{x^x — s)h{x)dx 



min(6.s) 



s/2 



dxh{x)K{x, s — x) 



+ I dxh{x)Bi{x, s) — / dxh(x)B2{x, s). 

J s/2 



(2.2.5) 



The kernels K,Bi, and B2 are bounded and continuous functions. Therefore, ii b < 00 and h{x) = for 
X > b, H2.2.5I) implies: 

j*b 

\h{y)\<c \h{x)\dx + c \h{x)\dx, 

J2y Jy 

where c > is a constant which bounds the kernels 2K, 2Bi and 2,62 from above and 2y ~ s. From the 
above inequality one gets 

. max |/i(y)| < ce max |/i(y)|, (2.2.6) 

b—€<y<b b—€<y<b 

where e, < e < 5, is sufficiently small so that ce < 1 and b — e < 2b — 2e. Then inequality (|2.2.6|l implies 
h{x) — Oifb— e<x<b. Repeating this argument, one proves, in finitely many steps, that h{x) — 0, 
< a: < 5. 

Theorem 12.2. II is proved. □ 

The proof of Theorem l2 . 2 . II is not valid if 6 = 00. The result is not valid either if b = 00. Let us give 
a counterexample. 

Theorem 2.2.2. There exists qi,q2 G and an h ^ 0, such that 

h{x)ipi{x,k)(p2{x,k)dx ^0 Vfc > 0. (2.2.7) 
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Proof. Let qi and q2 are two potenetials in L14 such that Si{k) = S2{k) Vfc > 0, £1 and £2 have one 
negative eigenvalue — fc^, which is the same for £1 and €2, but si 7^ S2, so that qi ^ q2- Let h :— q2 — qi- 
Let us prove that H2.2.7|) holds. One has — k^fi, £2^2 — k'^V2- subtract from the first equation 
the second and get: 

—if" — k^ip + qi^p ~ h(p2, ip :— ipi — (p2, ^{0, k) = (y3'(0, fc) — 0. (2.2.8) 
Multiply H2.2.8|) by ipi, integrate over (0, 00) and then by parts to get 

/•CO 

/ hip2Pidx = {ipip'-^-ip'ipi)\^ = Vfc>0. (2.2.9) 
Jo 

At X = we use conditions 12.2.8|) . and at x = cx) the phase shifts corresponding to qi and q2 are the 
same (because 5*1 (fc) = S2{k)) and therefore the right-hand side of H2.2.9|l vanishes. Theorem 12.2.21 is 
proved. 

□ 



15 



Chapter 3 



Inverse problem with /-function as 
the data 



3.1 Uniqueness theorem 

Consider equation (|1.2.5|l and assume q G Li.i Then f{x, k) is analytic in C+. Define the /-function: 

m = qj^- (3.1.1) 

From (|3.1.1|l it foUows that I{k) is meromorphic in C+ with the finitely many simple poles ikj, I < j < J- 
Indeed, ikj are simple zeros of f{k) and /'(O, ikj) 7^ as follows from H1.2.4() . Using (|1.2.20|l . one gets 

a,:=Res..,,,/(fc) = ffliM._^, fc^. > 0; «o = (3.I.2) 

where Ima^ > 0, 1 < j < J, and Imap > 0, oq 7^ iff /(O) 0. We prove that if g e and /(O) — 
then /(O) exists and /(O) 7^ fTheorem 13. 1 . 31 below) . This is a fine result. 

Lemma 3.1.1. The I{k) equals to the Weyl function m{k). 

Proof. The m{k) is a function such that 9{x, k) + in{k)ip{x, k) E L^(R+) if Imfc > 0. Clearly 

/(x, k) = c{k)[6{x, k) + m{k)ip{x, k)], 

where c{k) 7^ , Imfc > 0. Thus I{k) = '^gl^f},lZ(k)i'(s>i) = "^('^) because of (ITT!^ and (ITOIl . □ 

Our basic uniqueness theorem is: 
Theorem 3.1.2. If qj G ii.i, j — 1,2, generate the same I{k), then qi — q2- 

Proof. Let p '■— q2 — qi, fj be the Jost solution H1.2.5|l corresponding to q^, w := /i — /2- Then one has 
—w" + qiw — k'^w — pf2, \w\ + \w' \ — o{l) , X +00. (3.1.3) 
Multiply H3.1.3|) by /i, integrate over R+, then by parts, using H3.1.3|) . and get 







pf2fidx^w'{Q)fi{Q)~w{0)f[{0) 

= /{(O, fc)/2(fc) - /^(0,fc)/i(fc) = fi{k)f2{k)[h{k) - hik)] = 0. (3.1.4) 
By property C+ fTheorem 12. l.l|l . p(x) = 0. □ 
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Remark: If e j = 1,2, and (*) - l2{k)\ < ce-2'»I^^^ where k = |fc|e*argfe^ v|fc| > 

0, < argfc < TT, then qi{x) = q2{x) for almost all x € (0, a). This result is proved in |GSlj for qj G Lj^^. 
Our proof is based on (|3.1.4|l . from which, using (*), one gets (**) pf2fidx = ©(e'^'^ '=). Note 
that /1/2 = {I + V*)e'^^^^ , where V* is the adjoint to a Volterra operator (see the formula below (2.1.4)). 
Thus, (**) can be written as (***) pie'^'''^dx = Oie''^''^^''), where pi {I + V)p. Formulas (6.5.6) 
(see Chapter 6 below) and (***) imply pi = for almost all x G (0, a). Since F is a Volterra operator, 
it follows that p = for almost all x £ (0, a), as claimed. 

Theorem 3.1.3. If q e and /(O) = 0, then /(O) exists and /(O) ^ 0. 

Proof. Let us prove that f{k) = ikAi{k), li(0) ^ 0, Ai e'''* Ai{t)dt, and Ai e L\R+). Let 

Ai{t) := A{s)ds, A := f{k) - 1, and A{y) = A(0, y), where A(x, y) is defined in (|1.2.25|l and A(?/) e 
L^(R+) by (|1.2.26l) . Integrating by parts, one gets A{k) = -exp{ikt)Ai{t)\'^ + zfcAi = ikAi - 1. Thus 
/(A;) = ikAi. The basic difficulty is to prove that Ai e L^(M_|_). If this is done, then limfc^o = /(O) 
exists and /(O) = iAi{Q). To prove that /(O) ^ 0, one uses the Wronskian formula (|3.2.2|l with x — Q: 
f{—k)f'{0, fc) — /(fc)/'(0, — fc) = 2ik. Divide by k and let A; 0. Since existence of /(O) is proved, one gets 
-/(0)/'(0,0) = i, so /(O) ^ 0. We hav e used here the existence of the limit Muik^o f'{0,k) = /'(0,0). 
The existence of it follows from H1.2.25() : 

/■oo 

f{0,k) = ik-A{0,0)+ I A,(0,y)e^'=^dy, (3.1.5) 

and 







/'(0, 0) = -A(0, 0) + / A,(0, y)dy. (3.1.6) 







From H1.2.27|l one sees that Ax{0,y) € L^(R+). Thus, to complete the proof, one has to prove Ai e 
L^(R+). To prove this, use (|1.5.13|l with a; = and l|l. 5.1111 . Since f{ikj) — 0, one has A{ikj) = -1. 
Therefore H1.5.13|l with x ^ yields: 







A{y)+J A{t)Fsit + y)dt + Fsiy) = 0, y > 0. (3.1.7) 
Integrate (|3.1.7|l over {x,oo) to get: 

/•oo 

Ai{x)+J A{t) j Fs{t + y)dydt+ I Fs{y)dy = 0, x > 0, (3.1.8) 
where Fs{y) G L^(M+). Integrating by parts yields: 

/>oo 

A{t) Fs{t + y)dydt^A^{0) I F,{y)dy - I A,{t)Fs{x + t)dt. (3.1.9) 



Because = /(O) = 1 + A{y)dy, one has yli(O) = -1. Therefore (ISXH|l and imply: 

Ai{x)- Ai{t)Fs{x + t)dt^O, x>0. (3.1.10) 



"'0 

From this equation and from the inclusion Fs{t) G L^(]R+), one derives Ai E L^(R+) as follows. Choose 
a T(t) G C^(M+) such that \\F, - r||ii(R^) < 0.5, and let Q ■= F^ - T. Then ()j.l.lO|l can be written 
as: 

/•oo /"OO 

M{x)- I Q{x + t)Ai{t)dt = a{x) := T{x + t)Ai{t)dt, x>0. (3.1.11) 
Jo Jo 

Since T G C^(R+) and A G ii(M+), it follows that Ai is bounded. Thus a G ii(M+). The operator 

QAi := q{x + t)Ai{t)dt has norm ||Q||li(r^)^li(r^) < 0.5. Therefore equation (|3.1.3|l is uniquely 

solvable in iHK+) and Ai G Li(M+). Theorem IXOl is proved. 

□ 
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3.2 Characterization of the /-functions 



One has 



T Td.^ ^ i^M mm k 



where the Wronskian formula was used with x ~ 0: 



fix, k)f{x, k) - fix, k)f'{x, k) = 2tk. (3.2.2) 
From p.2.21|l and H3.2.1|) with k = VA, one gets: 

- lmI{V\)dX = dp, A > 0. (3.2.3) 

TT 

The I{k) determines uniquely the points ikj, 1 < j < J, as the (simple) poles of I{k) on the imaginary 
axis, and the numbers Cj by H3.1.2|l . Therefore /(fc) determines uniquely the spectral function /o(A) by 
formula (|1.2.21() . The characterization of the class of spectral functions p(A), given in Section lTBl induces 
a characterization of the class of /-functions. 

The other characterization of the /-functions one obtains by establishing a one-to-one correspondence 
between the /-function and the scattering data S p.2.17|) . Namely, the numbers kj and J, I < j < J, 
are obtained from /(/c) since ikj are the only poles if I{k) in C+, the numbers Sj are obtained by the 
formula (see (|1.2.16l) and ^T^): 

^ 'Mk, 

if fik) is found from /(fc). Finally, /(fc) can be uniquely recovered from /(fc) by solving a Riemann 
problem. To derive this problem, define 

ii'(fc) := n T:-V if /(0)<oo, (3.2.5) 
k + ikn 

and 

k 

woik):^- -wik), if/(0) = oo, K^kj Vj. (3.2.6) 

k ~\~ tK, 

Assumption H3.2.5|) . means that /(O) ^ 0, and H3.2.6|) means /(O) = 0. 
Define 

hik) ■.= w-\k)f{k), /(0)<oo (3.2.7) 



hoik) ■.= Wo\k) fik), /(0) = oo. (3.2.8) 



Write as /(fc) = 



or 



(fe) /(-fe) 

h+ik) ^ gik)h^ik), -oo < fc < oo, (3.2.9) 

where h+ik) :— /i(fc) is analytic in C+, h+ik) ^ in C+, the closure of C+, /i(oo) = 1 in C+, ft._(fc) := 
hi—k) has similar properties in C_, 

gik) > for fc > 0, gik) is bounded in a neighborhood of fc = and has a finite limit at fc = 0. From 
l|3.2.9|l and the properties of h, one gets: 
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and 

f{k) ^w{k)h{k), lnifc>0. (3.2.12) 

In Section we prove: 

- / [I{k) - ik]e~'^'dk = - y r,e'=^*, t < 0, (3.2.13) 



where = —iaj- Taking t —> -co in 1)3.2.13(1 . one finds step by step all the numbers rj, kj and J. If 
/(O) < oo, then ro = 0. Thus the data ((1.2.17() are algorithmically recovered from I{k) known for all 
fc > 0. 

A characterization of S is given in Section 15.51 and thus an implicit characterization of I{k) is also 
given. 

3.3 Inversion procedures. 

Both procedures in Section lX^ which allow one to construct either p(A) or S from I{k) can be considered 
as inversion procedures I ^ q because in Chapter 0| and Chapter |S1 reconstruction procedures are given 
for recovery of q{x) from either p(A) or S. All three data, I{k), p(A) and S are equivalent. Thus, our 
inversion schemes are: 

/(fc) ^ p(A) ^ g(x), (3.3.1) 

I{k) ^ S ^ q{x), (3.3.2) 
where ()1.5.1|l gives the details of the step p(A) q{x), and (|1.5.10|) gives the details of the step S =5> q{x). 

3.4 Properties of I{k) 

In this section, we derive the following formula for /(fc): 
Theorem 3.4.1. One has 

^ ci ■ r°° 

I(k)=ik + Y^ \ +S(fc), a{k)= a{t)e^*dt, (3.4.1) 

where fco, ImaQ > if and only if f(0) = 0, aj are the constants defined in (|3.1.2() . Imaj > 0, 1 < j < J, 
a{t) e Li(R+) if fiO) ^0 andqe Li,i, a(t) S L^(K+) j//(0) = and g e ii,3(M+). 

We prove this result in several steps which are formulated as lemmas. Using ((1.2.25|l one gets 

ik-A{0,0)+ j;^ A,{0,y)e'''ydy 



I{k) 

1 + A{k) 

A(y):=A(0,y), I(fc):= / A{y)e^''ydy. 

One has (cf. ^T^) 



(3.4.2) 







f(k) = l + Aik):=f„{k)wik)-^,wik):=Y[^—^,K^k, Vj (3.4.3) 

k + IK fc + ikj 

j=i J 

/o(fc)^OinC+, /o(c^) = l, (3.4.4) 

/o(fc) is analytic in C+, the factor in (|3.4.3|l is present if and only if /(O) = 0, and w{k)-£^ wo(fc). 
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Lemma 3.4.2. ///(O) ^ and q e then 

POO 

/o(fc) = l + 6o(fc), bo{x)eW''\R+), ||&olki.i(R+) := / (|6o| + l^oDdx < oo. (3.4.5) 

Jo 

Proof. It is sufficient to prove that, for any 1 < j < J , the function 

k + ik 



f{k) = 1 + / 5, (i)e"=* dt, g, e M/i-i(R+). (3.4.6) 



Since = 1 + and since ^(y) e W^i'i(R+) provided that g e Lia(M+) (see (ll.l>.:jHl - (ll.ll>7H . 

it is sufficient to check that 



/(fc) 

Jz 2 /u J 



git)e'^'dt, geW'-\R+). (3.4.7) 



^ - ^°°e^'=*[<5(t)-2fc,e-^-^*^(t)]dt, ^(t) := | J' (3.4.7') 



Note that 

k + ikj~J_^'^ "-J- ^vvj-, "^"^■^■[0, t<0 

One has f{ikj) = 0, thus 



/l 2 kj k "hk j i/Q k zkj 



ry foo 
Aiy)e-''^i'i e'^'"-'^^^" ds = e'^'hjis) ds, 



hj{s):^i A{y)e'''^'^y-'Uy^i A{t + s)e^''^* dt. (3.4.8 



where 



From one obtains l|n7|) since A{y) e 1^^^^^+). 

Lemma [3.4. 21 is proved. □ 

Lemma 3.4.3. ///(O) = and g e Li.2{R+), then (|TTK|) Wds. 

Proof. The proof goes as above with one difference: if /(O) — then fco = is present in formula (|3.4.1() 
and in formula (|3.4.8|l with ko = one has 



hois) = i I A{t + s)dt. (3.4.9) 

Thus, using p.2.2t)|l . one gets 

|ft-o(s)| (is < c / ds / dt I \q(u)\du 



OO /'OO 



OO />00 /"OO /"OO /"OO 



2c / ds dv \q{u)\du<2c / ds / \q{u)\udu 
Jo J§ Jv Jo 

Ac u^\q{u)\du < OO if g G Li,2(IR+), 



where c > is a constant. Similarly one checks that /io(s) e i-'^(R+) if g G Li,2{M.+ ). 

Lemma r-i.4.3l is proved. □ 

Lemma 3.4.4. Formula H3.4.1|) holds. 
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Proof. Write 



Clearly 



1 



k+i ttJ k+ikj 
k llj = l k—ikj 



k + i Yj k + ik. 
h 11 



/(fc) /o(fc) 
J 



3=0 



— , /co := 0, A;,- > 0. 



By the Wiener-Levy theorem |GRSt §17], one has 

1 



/o(fc) 



1+ / b{t)e"'*dt, b{t) €W''\R+), 



(3.4.5') 



where fo{k) is defined in (3.4.3). Actually, the Wiener-Levy theorem yields b{t) G -L^(R+). However, 
since Bq G iy^'-'^(M+), one can prove that b{t) S VF^'^(R4.). Indeed, h and 6o are related by the equation: 



which implies 



(l + 6o)(l + fc) = 1, Vfce 

h = -bo - bob, 



or 



5(i) = -5o(i) - / feo(t - s)b{s) ds := -5o - 6o * 6, 



(3.4.10) 



where * is the convolution operation. 

Since b'o e ^^1^+) and b € Li(M+) the convolution fog * 6 G Li(M+). So, differentiating l|3.4.1()|l one 
sees that b' £ L^{R+), as claimed. 

From the above formulas one gets: 



I{k) = {^k - A{0) + A,){1 + b)\l + J2 T^, 



J 

ik + c + y^ - — h a, 

j=0 J 



(3.4.11) 



where c is a constant defined in (|3.4.13|l below, the constants aj are defined in (|3.4.14|) and the function 
a is defined in H3.4.15|l . We wiU prove that c = (see (|3.4.17f) '). 
To derive (|3.4.11|l . we have used the formula: 



ikb = ik 



^ikt 



-m 



1 

ik 



e"'*b'{t)dt 



MO) ~ b', 



and made the following transformations: 



I{k) ^ik- A{0) - 6(0) -b' + Ai~ A{0)b + AibY^ 



Cjzk 
k zk-j 



^ c,[^(o)-f 6(0)] ^ ' m~mj) ^ . ^ 



3=0 



3=0 



3=0 



k xk'i 



(3.4.12) 



where 



g{k) := -b' + Ai~ A(0)6 + A^b. 
Comparing H3.4.12(l and H3.4.11|l one concludes that 

J 



c := -A(0) -6(0) -Hi^, 



(3.4.13) 



3=0 
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■=-Cj[k,+A{Q) + b{Q)~g{ik,)]. 



(3.4.14) 



m ■.= g{k)+Y. 



3=0 



g{k) ~ gjikj) ^ 



(3.4.15) 



To complete the proof of Lemma [3.4.41 one has to prove that c = 0, where c is defined in H3.4.13|l . 
This follows from the asymptotics of I{k) as k ^ oo. Namely, one has: 



A{k) 



^(0) _ 1^,^ 

ik ik 



(3.4.16) 



From 13.4. and H3.4.2I) one gets: 

I{k) = {ik^A{0) + Ai) 



m 

ik 



{ik-A{0)+Ai) (1 + ^ + 
ik 



"^1 
1 

k 



ik + o{l), k ^ +00. 



(3.4.17) 



From H3.4.17|l and H3.4.11|l it follows that c — 0. Lemma IS.4.41 is proved. 

Lemma 3.4.5. One has aj = irj, rj > 0, 1 < j < J, and ro — if f{0) ^ 0, and rg > if f{0) = 0. 
Proof. From H3.1.2() one gets: 



□ 



Tpir = iTTTT ^Tj, rj := > 0, j > 0. 



If j = 0, then 



flQ = Res/(fc) 



/'(0,0) 



fc=o ' - /(O) 



(3.4.18) 



(3.4.19) 



Here by ReSfe=o I{k) we mean the right-hand side of (|3.4.19|l since I{k) is, in general, not analytic in a 
disc centered at fc = 0, it is analytic in C+ and, in general, cannot be continued analytically into C_. 
By Theorem 13 . 1 . 31 the right-hand side of l|3.4.19(l is well defined and 



From H1.2.25|l one gets: 



/(0)=z/ A{y)ydy. 



(3.4.20) 



(3.4.21) 



Since A[y) is a real-valued function if q{x) is real-valued (this follows from the integral equation 
p.2.28|l . formula H3.4.21|l shows that 



/(O) 



<0, 



and (|3.4.2n|l implies 



Lemma |3.4.5l is proved. 



ro >0. 



(3.4.22) 

(3.4.23) 
□ 



One may be interested in the properties of function ait) in H3.4.1() . These can be obtained from 
H3.4.15|l and H3.4.5(l as in the proof of Lemma l3 . 4 . 21 and Lemma Hi 4. 31 
In particular, the statements of Theorem 13.4. II are obtained. 
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Remark 3.4.6. Even if q{x) ^ is compactly supported, one cannot claim that a(t) is compactly 
supported. 

Proof. Assume for simplicity that J = and /(O) 0. In this case, if a(t) is compactly supported then 
I{k) is an entire function of exponential type. It is proved in ^ p. 278] that if q{x) ^ is compactly 
supported, q E i^(M+), then f{k) has infinitely many zeros in C. The function /'(O, z) 7^ if f{z) = 0. 
Indeed, if f{z) = and /'(O, z) = then f{x, z) = by the uniqueness of the solution of the Cauchy 
problem for equation (|1.2.5II with k = z. Since f{x, z) ^ 0, one has a contradiction, which proves that 
/'(O, z) 7^ if f{z) = 0. Thus /(fc) cannot be an entire function if q{x) ^ 0, q{x) e L-'^(M+) and q{x) is 
compactly supported. □ 

Let us consider the following question: 

What are the potentials for which a{t) — in Ijtj. 4.1(1 ? 

In other words, let us assume 

^W^^fc + Eri^' (3.4.24) 



3=0 



and find q{x) corresponding to /-function (|3.4.24|) . and describe the decay properties of q{x) as a; ^ +00. 
We give two approaches to this problem. The first one is as follows. 
By definition 

/'(0,fc)=/(fc)/(fc), /'(0,-fc) = /(-fc)/(-fc), fceM. (3.4.25) 
Using i3.4.2,^ll and i^.'2.^^\^ one gets [/(fc) - I{-k)]f{k)f{-k) = 2ik, or 

= 7;;^- VfceM. (3.4.26) 

By (13.4. 18|l one can write (see (11.2. 21|l ) the spectral function corresponding to the /-function (|3.4.24|l 
{^/\ = k)■. 

*(A) = (T'"' , (3.4.27) 

\llUi'^kjrj5{\ + ^)d\, A<0, 

where ^(A) is the delta-function. 

Knowing dp{\) one can recover q(x) algorithmically by the scheme 1(1. 5. 
Consider an example. Suppose /(O) 7^ 0, J = 1, 

/(fc) = ^fc + ^ ^^k+ ^IlitLlhll = Jk + (3.4.28) 



Then (|3.4.27|l yields: 



fc - iki k"^ + kl \ k'^ + kfj fc^ + kf 



dp(A) (3.4.29) 
[2kiri6{X + k1)dX, A < 0. 



Thus yields: 



1 /"^ tiVA sin^/Axsin^/Ay sh{kix) sh{kiy) /o ^ Qn^ 

= "^ A + fcf va — TT^^""'''^. kT' ^'-'-'"^ 



TT 
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and, setting A = fc^ and taking for simplicity 2fciri = 1, one finds: 



2ri r°° dkk'^ sin(fca;) sin(fcy) 
' ~ Jo + kf 
2ri dk siii{kx) siTi{ky) 



k^ + kf 

-ki{x+y) 



k^ + kf 

ri f°° c?fc[cos k{x — y) — cos k{x + y)] 
Jo 

Il_ (p-ki\x-y\ 



(3.4.31) 



2fci 



, fci > 0, 



where the known formula was used: 



cos kx 
TT Jo fc^ + a? 



dk 



-a\x\ 



2a 



a > 0, 



X e 



Thus 



^~ki\x-y\ _ ^-ki{x+y) 



sh{kix) sh{kiy) 



ki ki 

Equation (|1.5.4|l with kernel H3.4.33|l is not an integral equation with degenerate kernel: 

-e-ki\t-y\ _ e-'^-i(*+2') sh{kit) sh{kiy)' 



Kix,y)+ / Kix,t) 



2fci/ri ki ki 

e-ki\x-y\ „ ^-k,ix+y) sh{kix) sh{kiy) 



(3.4.32) 
(3.4.33) 

(3.4.34) 



2fci/ri 



ki ki 



This equation can be solved analytically jRa) . but the solution is long. By this reason we do not give 
the theory developed in jRaj . but give the second approach to a study of the properties of q{x) given 
I{k) of the form (|3.4.28l) . This approach is based on the theory of the Riemann problem [H]. 

Equations H3.4.2()|l and (|3.4.28ll imply 



f{k)f{-k) = 



P + kl 
k'^ + I/? 



2 7 2 I 

i/^ := k^ + ri. 



The function 



Write t\AMl as 



/o(fc):=/(fc) Ji^^O in C+. 
k — iki 



J. -.k - iki k + iki k'^ + kj 

fo{k) , , /o(-fc) 



(3.4.35) 
(3.4.36) 



Thus 



k + iki ' 



k^ + 



k — iki k"^ + ul' 
./o(-fc) 



(3.4.37) 



The function fo{^k) 7^ in C_, /o(oo) = 1 in C_, so h := fa{-k) analytic in 
Consider H3.4.37(l as a Riemann problem. One has 



k'^ 



Therefore (see |^) problem (|3.4.37() is uniquely solvable. Its solution is: 



k + iki 
k + ivi 



h{k) = 



k — iui 

k — iki ' 



(3.4.38) 



(3.4.39) 
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as one can check. 
Thus, by ^3.4.36|l . 

^ (3.4.40) 

The corresponding S-matrix is: 
Thus 

/ [1 - S'(fc)]e''=^dfc = O (e-'^i^) for x > 0, (3.4.42) 

and 

Fix) = Fs{x) + Fdix) = O (e-'^i") . (3.4.43) 
Equation 1)1. 5. 13(1 implies A{x,x) = O {e^'^^^^), so 

q{x) = O (e-^''!^) , a; ^ +oo. (3.4.44) 

Thus, if /(O) 7^ and a(i) = then q[x) decays exponentially at the rate determined by the number 
fci, ki = min fc,-. 

1<3<J 

If /(O) = 0, J = 0, and a(i) = 0, then 

I{k) = ik+'-^, (3.4.45) 



f{k)fi-k) = ro > 0. (3.4.46) 



Let /o(fc) = i^tiiiW . Then equation i'iAAti^ implies: 

A:2 + l 



fo{k)fo{-k) = 7:2—3' ^0 '^o, (3.4.47) 

and /o(fc) 7^ in C+. 

Thus, since ind^^-^^ — 0, fo{k) is uniquely determined by the Riemann problem 13.4.47|l . 



One has: 



Mk) = ±^, M-k) ^ ' 



and 



k + ivQ ' k — ivq ' 

-2iiyo f°° e^^^'dk 

" ' 2iyoe~''"'^, a; > 0, 



and Fd{x) = 0. 
So one gets: 



2ti J-oo k ~~ Wo 



F{x) = Fsix) = 21^06-"°'', x>0. (3.4.49) 
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Equation p.5.13|) yields: 

A{x,y) + 2iyo I A{x,t)e"'°^*+y^dt = -2voe~''°''''+y\ y>x>0. (3.4.50) 



Solving H3.4.5U|I yields: 

Aix,y) = -2.oe--(^+'^)^^:^. (3.4.51) 

The corresponding potential (|1.5.12|) is 

q{x) = O [e-^""'') , x^oo. (3.4.52) 

li q{x) — O (e^*^^), fc > 0, then a{t) in (|3.4.1|l decays exponentially. Indeed, in this case 6'(t), ^i(y), b{t)^ 

Ai *h decay exponentially, so g{t) decays exponentially, and, by H3.4.15|l . the function ^'^^]rs'ii~^'^ '■= h 
with h{t) decaying exponentially. We leave the details to the reader. 
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Chapter 4 

Inverse spectral problem 



4.1 Auxiliary results 
4.1.1 Transformation operators 

If Ai and A2 are linear operators in a Banach space X, and T is a boundedly invertible linear operator 
such that AiT — TA2, then T is called a transformation (transmutation) operator. If A2/ = A/ then 
AiTf = AT/, so that T sends eigenfunctions of A2 into eigenfunctions of Ai with the same eigenvalue. 
Let £j = —-^ + qjix), j = 1,2, be selfadjoint in H L^(0, cx)) operators generated by the Dirichlet 
boundary condition a,t x = 0. Other selfadjoint boundary conditions can be considered also, for example, 
■u'(O) - hu{0) = 0, /i = const > 0. 

Theorem 4.1.1. Transformation operator for a pair {€1,^2} exists and is of the form Tf = {I + K)f , 
where the operator I + K is defined in (|1.2.24|l and the kernel K{x,y) is the unique solution to the 
problem: 

K^^{x,y) ~ qi{x)K{x,y) = Kyy - q2{y)K, (4-1.1) 



K{x,Q) = 



(4.1.2) 



1 f-' 

K{x,x) = 2 J ("^1 ^ q2)dy. 



(4.1.3) 



Proof. Consider for simplicity the case q2 = 0, qi = q. The proof is similar in the case q2 ^ 0. If 52 = 0, 
then (|4.1.3|l can be written as 

,(,)^2^^i^, 7^(0,0)^0. 



dx 



If hTf = mf and Tf = f + /; K{x, y)f dy, then 

-r+q{x)f + qTf - [K{x, x)f]' - ^^^^f - f dy 

= -/"-/ K{x,y)fyydy ^ - Kyyfdy-K{x,y)f' + Kyf 



(4.1.4) 



(4.1.5) 



Since / e D{£i), /(O) = 0, and / is arbitrary otherwise, H4.1.5|) imphcs H4.1.1|) . H4.1.2|) and (|4.1.4|l . 
Conversely, if K{x,y) solves ()4.1.1|l . H4.1.2I) and H4.1.4|) . then I + K is the transformation operator. To 
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finish tiie proof of Theorem l4.1.1l we need to prove existence of the solution to H4.1.1II . (|4. 1.211 and (|4.1.4|l . 
Let ^ = x + y,ri — x~y, K{x, y) := B{^, rj). Then (I4.1.1|l . (|4.1.2|l and H4.1.4I) can be written as 

Bin = ll (4^) B, 0)^1 q{s)ds, $) = 0. (4.1.6) 

Integrate (|4.1.6|l to get 

B^{i,^)^-a(^^+\ f q(i±A\B{i,r)dT. (4.1.7) 



(4.1.8) 



Integrate (|4.1.7|l with respect to ^ over (77,^) and get 



This is a Volterra integral equation which has a solution, this solution is unique, and it can be obtained 
by iterations. 

Theorem 14 .1.1 1 is proved. □ 
4.1.2 Spectral function 

Consider the problem 11.2.1(1 . The classical result, going back to Weyl, is: 

Theorem 4.1.2. There exists a monotone increasing Junction p{X), possibly nonunique, such that for 
every h E i^(0,oo), there exists h{X) G L'^{R;dp) such that 



\h\'^dx^ \hfdp{X), h{X) ■= lim / f{x)ip{x,Vx)dx, (4.1.9) 
J -00 Jo 

where the limit is understood in L^(R, dp) sense. If the potential q in (|1.2.1|l generates the Dirichlet 
operator i in the limit point at infinity case, then p{X) is uniquely defined by q, otherwise p{X) is defined 
by q nonuniquely. The spectral function of I has the following properties: 



■0 9X3/2 



e'"' dp(A) < c», p(A) = — — + o(Ai/2)^ A ^ +00. (4.1.10) 



-00 Stt 

Theorem 4.1.3. (Weyl). For any X, ImX ^ 0, there exists m(A) such that 

W{x, A) ip{x, A) + m{X)ip{x, A) G L'^{R+). (4.1.11) 

The function m(A) is analytic in C+ and in . 

The function m(A) is called Weyl's function, or m-function, and W is Weyl's solution. Theorem 14. 1.21 
and Theorem 14 . 1 . 31 are proved in |M| . 

4.2 Uniqueness theorem 

Let p(A) be a non-decreasing function of bounded variation on every compact subset of the real axis. Let 
h G iQ(R+), where L'^(R+) is a subset of L^(R+) functions which vanish near infinity. Let ipo := '^'"(^^) 
and 

/•oo 

H{X)^ I h{x)Lpo{x,X)dx. (4.2.1) 
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Our first assumption Ai) on p(A) is: 

H^{X)dp{X) =0, => h{x) = 0. (4.2.2) 



This implication sliould hold for any h S Lq(R+). It holds, for example, if dp{X) 7^ on a set which has 
a finite limit point: in this case the entire function of A, H(X), vanishes identically, and thus h = 0. 

Denote by P a subset of p(A) with the following property: if pi,p2 & V , v :— pi — p2, and H :— 
{H{X) : h e C^(R+)}, where H{X) is defined in gTII, then 

I J H'^{X)dv{X) =0 MH eH^^ v{X) = 0. (4.2.3) 

Our second assumption A2) on p{X) is: 

pCV. (4.2.4) 

Let us start with two lemmas. 

Lemma 4.2.1. Spectral functions p{X) of an operator £q = ^357 + in the limit-point at infinity 
case belong to V. 

Proof. Let pi, p2 be two spectral functions corresponding to £1 and £2, £j — £qj, j — 1,2, v ~ pi~ p2 and 
(*) H^{X)dv = V/i e Lo(M-|_). Let I + V and I + W be the transformation operators corresponding 
to £1 and £2 respectively, such that 

ipo = iI + V)^i = iI + W)ip2, (4.2.5) 
where (pj is the regular solution (|1.2.1|l corresponding to qj. Condition (*) implies 

\\{I + V*)h\\^\\{I + W*)h\\ \fheL\0,b), (4.2.6) 

where, for example, 

Vh= f V{x,y)h{y)dy, V*h^ f V {y , x)h{y)dy . (4.2.7) 
It follows from (|4.2.6(l that 

I + V* ^U{I + W*), (4.2.8) 

where U is & unitary operator in L^{Q, b). Indeed, U is an isometry and it is surjective because I + V* 
is. 

To finish the proof, one uses Lemma 14.2.21 below and concludes from (|4.2.8|l that V* = W* , so 
V = W, ipi = ip2, and qi = q2 := q. Since, by assumption, q is in the limit-point at infinity case, there 
is only one spectral function p corresponding to q, so pi = p2 = P- ^ 

Lemma 4.2.2. If U is unitary and V and W are Volterra operators, then (|4.2.8|) implies V — W. 

Proof. From one gets I + V + {I + W)U* . Since U is unitary, one has [I + V){I + V*) = 

(/ + W){I + W*). Because V \s a. Volterra operator, {I + V)^^ = / + Vi, where Vi is also a Volterra (of 
the same type as V in (|07t V Thus,(/ + Vi){I + W) = {I + V*){I + W^), or 

Vi+W + ViW = V* +Wl + V*W'^ (4.2.9) 

The left-hand side in 1)4.2.9(1 is a Volterra operator of the type V in 14.2.7(1 . while its right-hand side is 
a Volterra operator of the type V* . Since they are equal, each of them must be equal to zero. Thus, 
Vi{I + W) = -W, or (/ + Vy^{I + W)^I,oyV ^W. □ 

Theorem 4.2.3. (Marchenkoj The spectral function determines £q uniquely. 
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Proof. If iq^ and £q^ have the same spectral function p(A), then 

/OQ nOO 
\H,{X)\^dp^ \H,{X)\^dp yheLl{Q,b), (4.2.10) 
-C30 J — oo 

where 

Hj{X) := i h{x)(pj{x,k)dx, k = V^, j = 1,2. 
Jo 

Let I + K he the transformation operator Lp2 = {I + K)ipi, and g := {I + K*)h. Then H2 = (/i, (^52) = 
{h,(I + K)'.pi) = {g,(pi). From (|4.2.10|) one gets ||/i|| = + K*)h\\. Thus I + K* is isometric, and, 
because is a Volterra operator, the range of / + K* is the whole space L^{0, b). Therefore / + K* 
is unitary. This imphes K* = 0. Indeed, (/ + K*)-^ = I + K (unitarity) and (/ + K*)-^ = I + V* 
(Volterra property of A'*). Thus K = V* , so K = V* = 0. Therefore ip2 — fi and qi =52, so iq-^ — iq^. 
□ 

Remark 4.2.4. If pi = cp2, c = const > 0, then the above argument is applicable and shows that c 
must be equal to I, c = 1 and qi =52- Indeed, the above argument yields the unitarity of the operator 
\fc(I + K*), which implies c — 1 and K* — 0. 

Here the following lemma is useful: 

Lemma 4.2.5. If bl + Q = 0, where b — const and Q is a compact linear operator, then b — and 
Q = 0. 

A simple proof is left to the reader. 



4.3 Reconstruction procedure 

Assume that p{X), the spectral function corresponding to £q, is given. How can one reconstruct £q, that 
is, to find q{x)l We assume for simplicity the Dirichlet boundary condition at a; = 0, but the method 
allows one to reconstruct the boundary condition without knowing it a priori. 

The reconstruction procedure (the Gel'fand-Levitan or GL procedure) is given in (|1.5.1|l - (|1.5.4|) . Its 
basic step consists of the derivation of equation 11.5.4|) and of a study of this equation. 

Let us derive ()1.5.4|l We start with the formula 

^{x, VA)^(y, Vx)dp{X) = 5{x - y), (4.3.1) 

and assume that L{x, y) is a continuous function of x, y in [0, b) x [0, b) for any b G (0, 00). 
li < y < X, one gets from (|4.3.1|) the relation: 



ip{x, VX)ipiy, VX)dp{X) =0, Q<y <x. (4.3.2) 
Using Ijl. 2.24(1 . one gets Lpo = {I + A:)"V- Applying (/ + to (^(y, VA) in (|4.3.2|l . one gets 

0= / (^(a:,\/A)(^o(y,\/A)dp:=/(a^,y), Q<y<x. (4.3.3) 
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The right-hand side can be rewritten as: 

poo 

I{x,y)= I {ipo + Kipo{x)ipo{y,VX)d{p-- po) 

— L{x,y)+ I K{x, s)L{s,y)ds + 5{x — y) + / K{x, s)S{s — y)ds 
Jo Jo 

= L{x,y)+[ K{x,s)L{s,y)ds + K{x,y), 0<y<x. (4.3.4) 



From H4.3.3|) and H4.3.4|) one gets, using continuity a,t y = x, equation H1.5.4|) . 

In the above proof the integrals H4.3.2|l - (|4.3.4|l are understood in the distributional sense. If the first 
inequality (|4. 1.10(1 holds, then the above integrals over (—00, n) are well defined in the classical sense. 
If one assumes that the integral in 1(4. 3. 5() converges to a function L{x) which is twice differentiable in 
the classical sense: 

L{x) :— lim Ln{x) := lim [ ^—^^^—^da{X), 

n~too n-^oo 2A 

then the above proof can be understood in the classical sense, provided that (*) sup„ ^gj-^ {,) |L„(a;)| < 
c(a, b) for any —00 < a < b < 00. If p(A) is a spectral function corresponding to i, then the sequence 
Ln{x) satisfies (*). It is known (see L ) that the sequence 



^n{x,y) 



^ix,V-xMy,V-x)dpix)- r "'^^"^\""(^^) dpo(A) 



satisfies (*) and converges to zero. 

Lemma 4.3.1. Assume ()4.2.2|l and suppose that the function L{x) G Hl^^{^+), 

Hx):=[J^^^^d.iX). (4.3.5) 

Then equation ((1.5.4|l has a solution in L^{0,b) for any b> 0, and this solution is unique. 
Proof. Equation H1.5.4|) is of Fredholm-type: its kernel 

L{x,y)^L{x + y)-L{x-y), L{x, x) ^ L{2x), L(0) = 0, (4.3.6) 

is in H^{0,b) x H^{Q,b) for any b e (0, 00). Therefore Lemma [4.3.11 is proved if it is proved that the 
homogeneous version of (|1.5.4|l has only the trivial solution. Let 

h{y)+ [ L{s,y)h{s)ds = 0, 0<y<x, heL^{0,x). (4.3.7) 

Because L{x,y) is a real-valued function, one may assume that h{y) is real-valued. Multiply l|4.3.7|) by 
h{y), integrate over (0, x), and use (|4.2.1II . (|1.5.3|l and Parseval's equation to get 

/OO /"OO /"OO 

H^{X)da^\\hf+ H^{X)dp~\\hf^ H^{X)dp. (4.3.8) 
00 J — 00 J — 00 

From and it follows that h = 0. □ 

If the kernel K{x,y) is found from equation l(1.5.4|l . then q{x) is found by formula H4.1.4() . 
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4.4 Invertibility of the reconstruction steps 



Our basic result is: 



Theorem 4.4.1. Assume KT^ . and suppose L{x) e iJ;Q|,(M+). Then each of the steps in 

(|1.5.1() is invertible, so that (|1.5.5|l holds. 



Proof. 1. Step, p ^ L \s done by formula H1.5.3() . Let us prove L 
sponding to the same L{x,y), and v := pi — p2, then 



Multiply BXT|) by h{x)h{y), heC^{R+), use (jOTTl and get 



/OQ 
H^{\)dv{\) \JH e n. 
-00 



By (14.2.3(1 it follows that v = Q, so pi = p2. Thus L 



p. If there are pi and /92 corre- 

(4.4.1) 

(4.4.2) 
□ 



2. Step. L ^ K is done by solving H1.5.4|l . Lemma [4.3.11 savs that K is uniquely determined by L. Let 
us do the step K ^ L. Put y = x in (|1.5.4|l . use H4.3.5|l and (|4.3.6|l and get: 



L{2x)+ I K{x,s)[L{x + s) - L{x - s)]ds = -K{x,x), 



or 



L{2x) 



2x 



K{x,y - x)L{y)dy 



K{x, X - y)L{y)dy = -Kix, x). 



(4.4.3) 



(4.4.4) 



This is a Volterra integral equation for L{x) which has a solution and the solution is unique. Thus the 
step if => L is done. The functions L{x) and K(x^x) are of the same smoothness. □ 



3. Step. K ^ q\s done by formula (|4.1.4|l . q{x) is one derivative less smooth than K{x, x) and therefore 
one derivative less smooth than L{x). Thus q E Lf^^{M.^). The step q ^ K is done by solving the 
Goursat problem 1(4. 1.1(1 . I(4.1.2() . I(4.1.4() (with 92 = 0), or, equivalently, by solving Volterra equation 
1(4.1.8(1 . which is solvable and has a unique solution. The corresponding K{x,y) is in Hl^^{M.+ x 



q e LL(K+). 

Theorem 14.4. II is proved. 



if 
□ 

□ 



Let us prove that the q obtained by formula 14.1.4(1 generates the function Ki{x,y) identical to the 
function K obtained in Step 2. The idea of the proof is to show that both K and Ki solve the problem 



((4.1.4(1 with the same qi = q and (72 = 0. This is clear for Ki. In order to prove it for K, 
it is sufficient to derive from equation lll.5.4() equations ((4.1.1(1 and 1(4.1.2(1 with q given by 1(4.1.4(1 . Let 
us do this. Equation l(4.1.2() follows from ((1.2.5(1 because L{x, 0) — 0. Define D :— ^ — ^ := d1 — dy. 
Apply D to ((1.5.4(1 assuming L{x,y) twice diffcrentiable with respect to x and y, in which case K{x,y) 
is also twice diffcrentiable. (See Remark (4. 4.311 . By 1(4. 3. 6() . DL = 0, so 

DK + -^[K{x, x)L{x, y)] + K^{x, x)L{x, y) 
dx 



Kxxix,s)L{s,y)ds 



K{x, s)Lyy{s,y)ds = 0. 
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Integrate by parts the last integral, (use (|4.1.2ll 'l. and get 



{DK){x,y)+ {DK){x,s)L{s,y)ds 

Jo (4.4.5) 

+ KL + {K, + Ky)L{x, y) + K{L.Ax, y) - L,{s, y)|,=,) =0, < y < x, 

where K = K{x,x), L ^ L{x ,y), K = K^ + Ky = K, and Lx{x,y) - Ls{s,y)\s=x = 0. Subtract 

from <|4.4.5|) equation 11.5.4|) multiplied by q{x), denote DK{x,y) — q{x)K{x,y) :— v{x,y), and get: 



provided that —q(x)L(x, y) + 2KL(x, y) = 0, which is true because of 14.1.4|) . Equation (|4.4.6|l has only 
the trivial solution by Lemma [4.3. II Thus = 0, and equation (|4.1.1|l is derived. 
We have proved 

Lemma 4.4.2. If L(x,y) is twice differentiable continuously or in L^-sense then the solution K{x,y) 
of p.5.4|l solves H4. 1.1(1 . I|4.1.2|l with q given by 14.1.4(1 . 

Remark 4.4.3. // a Fredholm equation 



in a Banach space X depends on a parameter x continuously in the sense lim^^g 11^(2; + h) — ^(a;)|| = 
0, lim^^o 11/(2; + h) — f{x)\\ = 0, and at x — xq equation 1(4.4.71) has N{I + A{xo)) ~ {0}, where 
N{A) = {u : An = 0}, then the solution u{x) exists, is unique, and depends continuously on x in some 
neighborhood of xq, \x — xo\ < r. If the data, that is, A{x) and f{x), have m derivatives with respect to 
X, then the solution has the same number of derivatives. 

Derivatives are understood in the strong sense for the elements of X and in the operator norm for 
the operator A{x). 

4.5 Characterization of the class of spectral functions of the 
Sturm-Liouville operators 

From Theorem on it follows that if holds and L{x) G if/„^(R+), then q £ Lf^^(R+). Condition 

1(4.2.31) was used only to prove L =^ p, so if one starts with a q G Lf^^M.^), then by diagram ((1.5.5() 
one gets L{x,y) by formula ((4.3.6(1 . where L{x) £ Hl^^R^). If ((4.2.3(1 holds, then one gets from L{x) a 
unique p(A). 

Recall that assumption Ai) is (4.2.2). Let ^3) be the assumption L{x) £ i/j™+^(K+). 

Theorem 4.5.1. If Ai) holds, and p is a spectral function of £q, q £ HJ^^M.^), then assumption A3) 
holds. Conversely, if assumptions Ai) and A3) hold, then p is a spectral function of £q, q £ 11"^^^-^). 

Proof If ^1) holds and q £ i/;™ (M+), then L{x) £ i7;;'+i(R+) by (4.1.4). If Ai) and A3) hold, then 
q £ if;™^(R+) by (1.5.2), because equation (1.5.4) is uniquely solvable, and (1.5.5) holds by Theorem 
4.4.1. □ □ 

4.6 Relation to the inverse scattering problem 

Assume in the section that q £ Li.i. Then the scattering data S are ((1.2. 17() and the spectral function 




(4.4.6) 



{I + A{x))u^ fix) 



(4.4.7) 



is ((1.2.21(1 . 
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Let us show how to get dp, given S. If S is given then Sj, kj and J are known. If one finds f{k) then 
dp is recovered because 



4fc| 1 



as follows from (|1.2.20|) and H1.2.1()|) . To find /(fc), consider the Ricmann problem 



(4.6.1) 



/(fc) = 5(-fc)/(-fc), fceR, = (4.6.2) 

which can be written as (see H3.4.3|l ^: 

/o(fc) = 5(-fc)^^^/o(-fc) if ind5(fc) = -2J, (4.6.3) 

fo{k) = 5(-fc)^^^ ^±l^!/o(_fc) if ind S{k) = -2 J - 1. (4.6.4) 

Note that w{—k) — if k E R. The function fo{k) is analytic in C+ and has no zeros in C+j and 
fo{—k) has similar properties in C_. Therefore problems (|4.6.3|l and (|4.6.4|l have unique solutions: 

/o(fc)=exp(i- r log[^MV-^W]^n ifind^(fc) = -2J, Imfc > 0, (4.6.5) 



27ri t — fc 



. 27ri t-k J (4.6.6) 

if indS'(fc) = -2J- 1, lmfc>0, 

and 

f(k) = fo{k)w{k) if indS'(fc) = -2 J, Imfc > 0, (4.6.7) 



f{k) = fa{k)w{k)- if indS'(fc) = -2J- 1, lmfc>0. (4.6.8) 

k + IK. 

One can calculate f{x) for fc > by taking fc = fc + iO in 14.6.7|) or (|4.6.8|l . Thus, to find dp, given S, 
one goes through the following steps: 1) one finds J, Sj, kj, 1 < j < J; 2) one calculates ind5(fc) := J . 
If J = —2J, then one calculates /(fc) by formulas (|4.6.5II . H4.6.7I) . where u;(fc) is defined in (|3.4.3|l . and 
Cj by formula 1)4. 6. and, finally, dp by formula 11.2.21|l . 

If J' = —2J— 1, then one calculates /(fc) by formulas 14.6.61) and 1)4. 6. 8|) . where k > is an arbitrary 
number such that k, ^ kj, 1 < j < J. If /(fc) is found, one calculates Cj by formula H4.6.1I) . and then dp 
by formula H1.2.21|l . Note that /o(fc) in l|4.6.6() depends on k, but /(fc) in H4. 6.8(1 does not. 

This completes the description of the step S ^ p. 

Let us show how to get S given dp{X). 

From formula jl.2.2lll one finds J, kj, Cj and |/(fc)|. If |/(0)| ^ 0, then |/o(fc)| = |/(fc)| if fc G R. 
Thus, if 1/(0)1 7^ 0, then log/o(fc) is analytic in C+ and vanishes at infinity. It can be found in C+ from 
the values of its real part log |/o(fc)| by Schwarz's formula for the half-plane: 

log/o(fc) = - r^^T^^^' I-^>0- (4.6.9) 

ITT t-k 

If /(O) ^ 0, then / = fow, so 

/(fc)=cxp(l ri^^LM^MU(fc), in,fc>o. (4.6.10) 
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If 1/(0)1 = 0, then the same formula (|4.6.1UI) remains vahd. One can see this because is analytic in 

C+, has no zeroes in C4., tends to 1 at infinity, and = |/(fc)| if fc S R. 

Let us summarize the step dp =4» S: one finds J, kj, Cj, calculates f{k) by formula (|4.6.10|) . and then 
S{k) — ^fVfA , and sj are calculated by formula H4.6.1(l . To calculate f{k) for fc > one takes k = k + iO 



in (|4.6.10() and gets: 



/(fc) = exp(l /" 



^og\f{t)\dt . 
— 7 ; + log \f{k)\ } w{k) 



\f{k)\w{k)expi-p r 

I ITT 



\og\m\dt]^ 

t-k J 



(4.6.11) 

fc > 0. 
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Chapter 5 

Inverse scattering on half-line 



5.1 Auxiliary material 
5.1.1 Transformation operators 

Theorem 5.1.1. If q E Li.i, then there exists a unique operator I + A such that H1.2.25|l - (|1.2.28|l 

hold, and A{x, y) solves the following Goursat problem: 

Axx - q{x)A^ Ayy,Q <x <y <oo, (5.1.1) 



A{x, x) ~ 




(5.1.2) 



lim A{x,y) = lim Ax{x,y) = lim Ay{x,y) = Q. (5.1.3) 

x+y — >oo x+y — ^oo ' x+y — ^oo 

Proof. Equations (|5.1.1|l and (|5.1.2|l are derived similarly to the derivation of the similar equations for 
K{x,y) in Theorem l4. 1.11 Relations H5. 1.3(1 follow from the estimates (|1.2.26|l - (|1.2.27|l . which give 
more precise information than ((5.1.3|1 . Estimates H1.2.26|l - l(1.2.28|l can be derived from the Volterra 
equation ((1.2.28|l which is solvable by iterations. Equation 1(1.2.28(1 can be derived, for example, similarly 
to the derivation of equation 1(4.1.8(1 . or by substituting 1(1.2.25(1 into ((1.2.5(1 . 

A detailed derivation of all of the results of Theorem 15 . 1 . II can be found in |M]. □ 

5.1.2 Statement of the direct scattering problem on half-cixis. Existence and 
uniqueness of its solution. 

The direct scattering problem on half- line consists of finding the solution -0 = ^j(r, k) to the equation: 

ip" + k^Tp - q{r)TP = 0, r > 0, (5.1.4) 
satisfying the boundary conditions at ?' = and at r — oo: 

V'(0) = 0, (5.1.5) 

^/,(r) = e'^sin(fcr + (5) + o(l), r ^ +oo, (5.1.6) 
where 5 = S{k) is called the phase shift, and it has to be bound. An equivalent formulation of 1(5.1.6(1 is: 

^ = l[e-»'='-_5(fc)e*^'n+o(l), r->oo, (5.1.7) 
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where S{k) = = e^*^^'^). Clearly 

4'ir,k)^'-[f{r,-k)-S{k)f{r,k)]=a{k)ip{r,k), a(fc) (5.1.8) 
2 j{k) 

where (p{r, k) is defined in (|1.2.1|l . see also (|1.2.10|l . From H5.1.8|) . \1.2.7\ and H5.1.6|) one gets 

^(r,fc) = l^^^sin(/fcr + (5(fc)) +o(l), r ^ oo. (5.1.9) 
k 

Existence and uniqueness of the scattering solution ijjir^k) follows from (|5.1.8|l because existence and 
uniqueness of the regular solution ^pir, k) follows from (|1.2.1|l or from 11.2.9|) . 

5.1.3 Higher angular momenta. 

If one studies the three-dimensional scattering problem with a spherically-symmetric potential q{x) — 
ci{r), X e M^, |x| r, then the scattering solution il>{r,a,k) solves the problem: 

[V^ + k^ - q{r)]ij = in (5.1.10) 

ip = e'''"-'' +A(a',a,k)- — + o ( ,r := \x\ oo,a' := -,a e 5"^ (5.1.11) 

r \r J r 

Here is the unit sphere in M'^ a g S'^ is given, A{a',a,k) is called the scattering amplitude. If 
q = q{r), then A[a' , a, k) = A{a' ■ a, k). The converse is a theorem of Ramm 0, p. 130. The scattering 
solution solves the integral equation: 

V' = e^''"-"- / g{x,y,k)q{y)ij{y,a,k)dy, g :^ f ^'^ . (5.1.12) 

Jk3 - y\ 



It is known that 



e'""" ^ } ^ —i" ' YtXa')Ye{a), a := -, Ui:= \ —Jg+i, 5.1.13 
— ' k r r V 2 2 

Yi{a) are orthonormal in LF'{S'^) spherical harmonics, Yi — Ygm, —£ <m<l, and summation over m in 
H5.1.13|l is understood but not shown, and Jt.{r) is the Bessel function. 
If q = q{r), then 

47r .ni^fXr, k) 



oo 

V = E ^i^'^^^Y,{a')Yda), (5.1.14) 

^ — ' k r 



where 



iP'; + k^t^i - q{r)i;i - ^^^^t^i ^ 0, (5.1.15) 

ije = e'^'sin (^kr- Y + Se^ +o{l), r ^ oo, (5.1.16) 

i/;^ = 0(/+i), r^O. (5.1.17) 

Relation (|5. 1.16(1 is equivalent to 

p»f (^+1) 

i^i = z [e^'"'' - e^'^^^.e^^T -t- o(l), r ^ oo, (5.1.18) 
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similar to IjS. 1.811 . which is (|5.1.18|) with £ = 0. li q = q{r), then the scattering ampHtude A{a',a,) — 
A{a' ■ a, k) can be written as 



A{a' ■a,k)=Yl Ae{k)Yi{a')Ye{a), (5.1.19) 
while in the general case — q{x), one has 

oo 

A{a', a,k) = Y^ Ai{a, k)Ye{a'). (5.1.20) 



li q = q{r) then Se in H5.1.18|l are related to Ae in H5.1.19|) by the formula 

Se^l- ^Ae. (5.1.21) 
In the general case q — q{x), one has a relation between S'-matrix and the scattering amplitude: 

S^I-—A, (5.1.22) 

2TTi ^ ' 

so that (|5.1.21|l is a consequence of H5.1.22|l in the case q = q{r) : Se are the eigenvalues of S in the 
eigenbasis of spherical harmonics. Since S is unitary, one has — 1, so Se = e^"^^^ for some real 
numbers 6i, which are called phase-shifts. These numbers are the same as in H5.1.17|) (cf. H5.1.18II 1. 
From H5.1.21|l one gets 

Ae{k)^^e'^'sixY{6A. (5.1.23) 
k 



The Green function ge{r, p), which solves the equation 

k^- ^'" , ^' ]g, = -5{r-p), ^ - ikge -> 0, (5.1.24) 



-T^ + k 2 — \ge = -S{r-p), — ^ 

ar r J or r^+oo 

can be written explicitly: 

g,(r,p) = l ^o/ Voi{kp)foi{kr) r>p, F,,, := ^ ^^^^ 25) 

i F^f^ foe{kr)foe{kp), r < p, ipot{kr) = ^fr+r^, 

and the function i/'f (r, A;) solves the equation: 

/•oo 

%l)i(r,k) ^ ue{kr) - I ge{r, p)q{p)ipi(p,k)dp. (5.1.26) 



The function Foe{k) is the Wronskian W[foe, 'Poi], Voe{kr) is defined in (|5.1.25|l and foe is the solution 
to (|5.1.15|1 (with q = 0) with the asymptotics 

foe =e''''' + o(l), r -> +oo, foeikr) = e'^^^ iue{kr) + ive{kr)), 

(5.1.27) 

«f := \l Y^£+i(M- 

Let <^^(r, k) be the regular solution to H5.1.15|l which is defined by the asymptotics as r ^ 0: 

Vi{rM= (21+1)!! ^^-^-^^^ 

Then 

ipe{r,k) = ae{k)ipe{r,k), 

|/,(0,fc)| . / V (5-1-29) 
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Lemma 5.1.2. One has: 

sup |a£(fc)| < oo, (5.1.30) 

£=0,1,2,... 

where k > is a fixed number. 

We omit the proof of this lemma. 

5.1.4 Eigenfunction expansion 

We assume that q € Li,i and h e Co°°(M+), eh = -h" + q{x)h, A = let g = v>i^-^Wy.V\) ^ y>x>0, 

be the resolvent kernel of ^ : {i — X)g ~ 6{x — y), gh g{^}h :— g{x, y, X)hdy, and fj — f{y, ikj). 
Then ^ = —gh + jigh- Integrate this with respect to A G C over |A| = iV and divide by 27rj to get 

^ f ghdX + — f ^^d\:^h^h. 



27r« J\\\=N 27ri J\x\=n A 

The function gh is analytic with respect to A on the complex plane with the cut (0, oo) except for the 
points A = 1 5: J 5: which are simple poles of gh, and limAr^oo I2 ~ 0, because \igh\ — o(l) as 

N 00. Therefore: 

1 —1 f 
h = / [5(A + iO)/i-g(A-iO)/i]dA + V * ghdX. (5.1.31) 

One has (cf. H1.2.10|l ): 

g(A + ^0)-.9(A~^0) /(-fc)/(y, fc) - /(gy, -fc)/(fc) 

2^ = ^^7(fc)P 

A: 



(p{x, k)(p{y, k), k — > 0. 



Also 



(//idA — — Res_)j^_j,2 g/i 

|A+fc?|=<5 

00 



Sj are defined in (|1.2.16() . and 



fjiy)Hy)dy ^^r^^i^ = Sjfj{x)hj, hj := / fjhdy, 



_ fix,ikj) ^ fjjx) 

no,ik,)- no,tk,)- 



'^(^'»fcj) = 7777f-rT-7777r-7rT- (5-1-32) 

Therefore 



H^) = (^y^ ^{y, k)hiy)dyj ip{x, k)--^jj^ + Sjfj{x)hj. (5.1.33) 
This implies (cf. I|1.2.21|l . (|1.2.20|l . (|1.2.16|l ^: 



(5(x -?/) = - ^ix,k)ip{y,k)jjj^ +22sjfjix)fjiy) 



\fikW 

^' J=i (5.1.34) 



ip{x,VJ)ip{y,VX)dp{X). 



We have proved the eigenfunction expansion theorem for h E C^(R+). Since this set is dense in i^(R+), 
one gets the theorem for h G L^(R+). 
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Theorem 5.1.3. If q E Li i, then H5.1.33|l holds for any h G i7^(R-)_) and the integrals converge in 
L^(M-l-) sense. Parseval's equality is: 

\\h\\l^(^^)=Y.'j\hj\^ + - h:^j^ h{yMy,k)dy. (5.1.35) 

5.2 Statement of the inverse scattering problem on the half- 
hne. Uniqueness theorem 

In Section [l. 31 the statement of the ISP is given. Let us prove the uniqueness theorem. 
Theorem 5.2.1. If qi,q2 G generate the same data (|1.2.17l) . then qi = q2. 

Proof. We prove that the data p.2.17|l determine uniquely I{k), and this imphes qi = q2 by Theo- 
rem |XII2 

Claim 1. If H1.2.17|l is given, then f{k) is uniquely determined. 

Assume there are fi{k) and /2(fc) corresponding to the data (|1.2.17|l . Then 

7§ = 7TTT' — <^<- (5-2.1) 

72 (fc) 72 (-fc) 

The left-hand side of H5.2.1(l is analytic in C+ and tends to 1 as |fc| ^ oo, /c G C+, and the right-hand 
side of H5.2.1|) is analytic in C_ and tends to 1 as oo, k G C_. By analytic continuation is 

an analytic function in C, which tends to 1 as |fc| oo, /c G C. Thus, by Liouville theorem, = 1, 

so A = /2. ' □ 

Claim 2. If H1.2.17|l is given, then /'(0,fc) is uniquely defined. 

Assume there are /{ (0, k) and /2(0, k) corresponding to H1.2.17|) . By the Wronskian relation H1.2.11|l . 
taking into account that fi{k) — f2{k) := f(k) by Claim 1, one gets 

[/{(O, k) - /^(O, k)]fi-k) - [/{ (0, -k) - /^(O, -k)]fik) = 0. (5.2.2) 

Denote w(fc) := fiiO,k) - f^{0,k). Then: 

w{k) w{—k) 



fik) f{-k)' 



keR. (5.2.3) 



The function jj^ is analytic in C+ and tends to zero as |fc| — > oo, fc G C+, and has similar 

properties in C_. It follows that = 0, so /((0,fc) = /2(0, /c). Let us check that is analytic in 
C+. One has to check that w{ikj) = 0. This follows from (|1.2.16(l : if f{k), Sj and kj are given, then 
f'{0,ikj) are uniquely determined. 

Let us check that w{k) — > as |fc| — > oo, k G C+. Using (|3.1.5|l it is sufficient to check that A{0,0) 
is uniquely determined by /(fc), because the integral in (|3.1.5|l tends to zero as |fc| — > oo, k G C+ by the 
Riemann-Lebesgue lemma. From H1.5.16|l . integrating by parts one gets: 



ik ik 

Thus 



OO 



f{k) = 1 - 1 / e^'^yAyiO, y)dy. (5.2.4) 







A(0, 0) = - lim [ik{f{k) - 1)]. (5.2.5) 

k — >oo 

Claim 2 is proved. □ 

Thus, Theorem l5.2.1l is proved. □ 
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5.3 Reconstruction procedure 

This procedure is described in (|1.5.10(l . 

Let us derive equation H1.5.13|l . Our starting point is formula (|5.1.34fl : 

r°° 2k^dk 

^{x,kMy,k)--^j^ +J2sjf,{x)f,iy) ^0 y>x>0. (5.3.1) 



From l)l.!^.10|l and l)l.!^.i>5|l one gets: 

k(p{x, k) 



= sin(fca; + 5) + I A[x^ y) sm{ky + 5)dy 



\fm ' , ^ . ^532) 

= {I + A)s\u{kx + 5), 5 = 6{k). 
Apply to H5.3.1|) operator (/ + A)"^, acting on the fimctions of y, and get: 

^ ^^^j^smiky + 6)dk + J2sjf,{x)e-''-y ^0, y>x>0. (5.3.3) 

From IjS. 3.2(1 . 1(5. 3. 3(1 . and (|1.2.25|l with k — ikj, one gets: 



{I + A) J sm{kx + 6) sin(fcy + 5)dk 



(/ + ^)^Sje-'=j(^+^) ^0,y>x>0. 



One has 



2 1 

— / sin(fca: + ^) sin(fcy + (5)dfc = — / cos[k{x — y)]dk 

Jo ' Jo 



1 r°° 

= S{x -y) + — [1 - s{k)]e'''^'-'+y'>dk. 

27r J„oo 



(5.3.4) 



1 1 

-/ cos\k(x + y) + 26(k)]dk^S(x,y) / (e^'^^'''> - l)e'''^'-'+y'>dk (5.3.5) 



From (|1.5.11|) . H5.3.4I) and H5.3.5|) one gets (|1.5.13|l . By continuity equation (|1.5.13|l . derived for y > x > 
0, remains valid for y > x > 0. □ 

Theorem 5.3.1. // q G Li.i and F is defined by (|1.5.11|) then equation (|1.5.13|l has a solution in 
L^{Wx) n L°°(Rx), := [x,oo), for any a; > 0, and this solution is unique. 

Let us outline the steps of the proof. 

Step 1. If q £ Li^i, then F{x), defined by ((1.5.11|l satisfies the following estimates: 

\F{2x)\ < ca{x)\, \F{2x) + A{x,x)\ < ca{x), |F'(2x) - ^| < ca^ix), (5.3.6) 
where a{x) is defined in p.2.2ti|) . and 

\\F\\lhr+) + l|i^llLi(R+) + l|i^llL-(K+) + \\xF'ix)h^iR^) < oo, (5.3.7) 
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oo /'OO 



\F{s + y)\dsdy < oo. (5.3.8) 

Jo Jo 
Step 2. Equation 

/•oo 

{I + F^)h h{y) + h{s)F{s + y)ds^O, y>x>0 (5.3.9) 



is of Fredholm type in L^(Rx), L'^i^x) o.nd in i°°(R^). It has only the trivial solution h = 0. 

Using estimates ()5.3.6|l - H5.3.8|l and the criteria of compactness in i^(R^), p — 1, 2, oo, one checks 
that Fx is compact in these spaces for any a; > 0. The space fl L°° C because ||ft.||2 < ||^||i||^||oo, 
where \\h\\p := ||/i||ip(R^). We need the following lemma: 

Lemma 5.3.2. Let h solve l^TT^ . If h e L^{Rx), then h £ L'^ -.^ i^(M^). If h £ , then 

h £ L°°. Ifh£ L^, then h £ 

Proof If h solves (I5.3.9|l . then \\h\\oo < ll^lli sup^>2i; l-P'C?/)! < c(a;)||/i||i < oo, where c{x) ^ as 
a: ^ oo. Also \\h\\2 < dya^{^^)\\h\\l < ci{x)\\h\\l < oo, Ci(x) — > as x — > oo. So the first claim is 
proved. Also \\h\\i < \\h\\i sups>^ \F{s + y)\dy — C2{x)\\h\\i, C2{x) — > as a; — > oo. li h £ L^, then 

ll^lloo < ||^||2supy>^ (J^ \F{s + y)\'^ds)^ = C3{x)\\h\\2, csix) ^ as a; ^ oo. □ 
Lemma 5.3.3. If h £ solves (|5.3.9|l and x >0, then h — 0. 

Proof. By Lemma [5. 3. 21 h £ nL°". It is sufficient to give a proof assuming a: = 0. The function F{x) 
is real-valued, so one can assume that h is real- valued. Multiply (|5.3.9f) by h and integrate over (x,oo) 
to get 

WH^ + TT S{k)]h^{k)dk + J2sj ( e-'''''h{s)ds\ =0, 

h:^ e'^'h{s)ds, 

J X 

where \\h\\ = M ^^^^ one gets ^■°^^h\k)dk = 0. Also, \2- j^^S{k) h\k)dk\ < ^ J^, \h{k)\^dk = 
Therefore l|5.3.1()|) implies s^ = 0, = hj := he-^'^ds, l<j<J, and 

(h,h) = (h,S{-k)h{~k)), (5.3.11) 



where {h,g) := h{k)g{k)dk. Since h is real valued, one has h{—k). The unitarity of S implies 
S'-i(fc) = S{-k) = S{k), k £ M,jind \\ S {- k)h{- k)\\ = \\h{-k)\\. Because of J5.3.11|l . one has equality 
sign in the Cauchy inequality {h, S{—k)h{—k) < This means that h{k) — S{—k)h{—k), and 

H1.2.17|l implies 

h{k) h{-k) 



f{k) f{-~ky 



fceK. (5.3.12) 



Because hj — 0, one has h{ikj) = 0, and if /(O) =/= 0, then jpjy is analytic in C+ and vanishes as |fc| — > oo, 
k £ oo, k £ C+. Also is analytic in C_ and vanishes as \k\ — > oo, k £ C^. Therefore, by analytic 

continuation, is analytic in C and vanishes as — > oo. By Liouville theorem, jpy = 0, so h{k) — 

and h = 0. If /(O) = 0, then, by Theorem 13. 1.31 f{k) — ikAi(k), Ai{Q) ^ 0, and the above argument 
works. □ 
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Because is compact in L'^{'Rx), the Fredholm alternative is applicable to (|5.3.9|) . and Lemma 15. 3. 31 
implies that (|1.5.13l) has a solution in L^(M2:) for any x > 0, and this solution is unique. Note that the 
free term in H1.5.13|) is -~F{x + y), and this function of y belongs to L'^{Rx) (cf. H5.3.8|) 'l. Because 
is compact in L^(Rx), Lemma 15.3.31 and Lemma |5.3.2I imply existence and uniqueness of the solution 
to (|1.5.13|) in L'^{Rx) for any a; > 0, and F{x + y) e L'^{Rx) for any x > 0. Note that the solution to 
H1.5.13|l in L^(Ma;) is the same as its solution in L'^(Rx)- This is established by the argument used in the 
proof of Lemma 15.3.21 

We give a method for the derivation of the estimates (|5.3.6|l - (|5.3.8|) . Estimate H5.3.8I) is an immediate 
consequence of the first estimate (|5.3.6f) . Indeed, 

/•oo poo / poo \ 2 

\F{s + y)fdsdy<(j^ m^x\F{s + y)\dy] 

(pOC poo \ ^ / nOO \ 2 

J J \q\dtdyj <cyj t\q{t)\dtj < oo. 

Let us prove the first estimate H5.3.6|) . Put in H1.5.13|) x — y: 

/>oo 

A{x,x)+ A{x,s)F{s + x)dx + F{2x)=0. (5.3.13) 



Thus 

poo 

\F{2x)\<\A{x,x)\+ \A{x,s)F{s + x)\dx. (5.3.14) 

J X 

From 11.2.26|l and H5.3.14|l one gets 

\F(2x)\ < ca{x) +c r a f^^\ \F{s + x)\ds 

1 ^ (5.3.15) 

< ca{x) + ccr(x) / + x)\ds < ca{x), 

J X 

where c = const > stands for various constants and we have used the estimate 

poo pOO 

sup/ \F{s + x)\ds< \F{t)\dt = c<oo. 

x>oJx Jo 

This estimate can be derived from 11.2.26|) . Write H1.2.26|l as 

/•oo 

A{x,z-x)+ A{x,t + x- z)F{t)dt + F{z)=0, z>2x>0. (5.3.16) 



Let us prove that equation H5.3.16|l is uniquely solvable for F in Lp(Rn), P — oo, p = 1 for all x > 
where A'^ is a sufficiently large number. In fact, we prove that the operator in (|5.3.16|l has small norm 
in LP{M.isf) in is sufficiently large. Its norm in L°°(M.n) is not more than 

poo poo 

sup / \A{x,t + X — z)\dt < c / \q{s)\ds 
z>nJn Jx+i-^ 



<c dt \q{s)\ds^c {s - N)\q{s)\ds < I 
Jn Jt Jn 

because q S Li,i. We have used estimate H1.2.26|l above. The function A{x,~x) G L'^{M.n), so our 
claim is proved for p — oo. Consider the case p — I. One has the following upper estimate for the norm 
of the operator in 1)5.3.16(1 in L^{Wn)'. supojy \A{x,t + x — z)\dz < supo^ /J ^ \A{x,x + v)\dv < 
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dv /^i. I^l^s — 2 /^(s ~ — > as a; — > oo. Also \A{x,z — x)\dz < oo. Thus equation 
H5.3.1t)|l is uniquely solvable in L^(RAr) for all a; > if ^ sufficiently large. In order to finish the 
proof of the first estimate (|5.3.6|l it is sufficient to prove that ||-F||loo(o.jv) ^ c < oo. This estimate is 
obvious for Fdix) (cf. (|1.5.11^ V Let us prove it for Fs{x). Using (ESI, ^3.4.5'|l . one 

gets 



where all the Fourier transforms are taken of ■ (IR+) functions. Thus, one can conclude that Fs{x) € 

k 



L°°(K+) if one can prove that / :— ^(^)^^('') jg fj^^ Fourier transform of L°°(R+) function. One has 



I = J^°°dyA{y) -•^"7-"°" and 

e'''''I{k)dk= I dyA{y) 



00 ik(x+y) _ ik{x-y) 



k 

dyA{y)iTT[l - sgn{x - y)\ = 2i7r / dyA{y). 



(5.3.18) 



From (|1.2.2fi|) it follows that A{y)dy e i°°(M+). We have proved that |1F||l=o(r^) + |1F|1li(r+) < 00. 
Differentiate (|5.3.13|) to get 



(5.3.19) 



2F'{2x) + A{x,x) - A{x,x)F{2x) + / A^{x, s)F{s + x)ds 

dA(x, x) 



/•oo 

/ A{x,s)F'{s + x)ds^O, A:= 

J X 



dx 



or 

a(x) 1 
F'{2x)^^^ + A{x,x)F{2x)--j [A^{x,s) ~ As{x,s)]F(s + x)ds. (5.3.20) 

One has x\q\dx < 00, x\A{x,x)\\F{2x)\dx < sup^y,Q{x\A{x,x)\) ■ \F{2x)\dx < c. Let us check 
that / := x\ J^[Ax{x, s) — As{x, s)]F{s + x)ds\ dx < 00. Use (|1.2.27|l and get / < c xa{x) a 
(^^) \F{s + x)\dsdx < c a{x)dx \F{y)\dy ■ sup^yQ .,^^ xa (^7^) < c < 00. The desired estimate 
is derived. 

The third estimate ^J^, \F'i2x) - ^\ < ca^{x) follows from ()5.ij.!^0|) because \A{x,x)\ < 
ca{x), \F{2x)\ < caix), and |A,(x,s) - A,{x,s)\\F{s + x)ds < cct{x) a {^^) |F(s + < 
ca^j x) \F {s + x)\ds < ca'^{x). The estimate \F {2x) + A{x , x)\ < ca{x) follows similarly from (|5.3.13() 
and (|1.2.26|1 . Theorem 15X11 is proved. 



5.4 Invertibility of the steps of the reconstruction procedure 

The reconstruction procedure is (|1.5.1()|l . 1. The step 5 ^ F is done by formula (|1.5.11|l . 

To do the step F ^ S, one takes x — )■ —00 in p.5.11|) and finds Sj, kj, and J. Thus Fd{x) is found 
and Fs — F — Fd is found. From Fs{x) and finds 1 — S{k) by the inverse Fourier transform. So >5'(fc) is 
found and the data S (see (|1.2.17|l l is found 

2. The step F ^ A is done by solving equation H1.5.13|l . By Theorem 15 . 3 . II this equation is uniquely 
solvable in (M.^) L°° (R^) for all a; > if (7 G that is, if F came from S corresponding to g G 

To do the step A F, one finds f{k) = 1 + A{0,y)e^''ydy, then the numbers ikj, the zeros of 

/(fc) in C+, the number J, 1 < j < J, and S{k) = ^j^^- The numbers sj are found by formula 11.2.16|l . 
where 

/"OO 

/'(0,i%) = -fc,-A(0,0)+ / A,{0,y)e-''^ydy. (5.4.1) 

Jo 
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Thus S and S ^ F hy formula (|1.5.11|l . We also give a direct way to do the step F. 
Write equation 11.5.13() with z = x + y,v = s + y,as 

{I + B^)F -.^ F{z) + A{x,v + x-z)F{v)dv = -A{x,z-x), z>2x>0. (5.4.2) 

J z 

The norm of the operator in ija; is estimated as follows: 



r r f V ~ z\ f° 
\Bx\\ < sup / \A{x, V + X — z)\dz < csup / cr I .t H — \ dz < c I 

v>oJo v>oJq V 2 / 



a{t)dt, (5.4.3) 



where cr{x) = \q{t)\dt and the estimate (|1.2.26|l was used. If Xq is sufficiently large then < 1 

for X > xo because a{t)dt ^Oasx^ooifgG Therefore equation H5.4.2I) is uniquely solvable 

in for all x> xq (by the contraction mapping principle), and so F[z) is uniquely determined for all 
z > 2x0. 

Now rewrite (15. 4. 211 as 



F{z) + / A{x,v + X ~ z)F{v)dv = ~A{x,z ~ x) - / A{x,v + x ~ z)F{v)dv. (5.4.4) 

J z J 2x0 

This is a Volterra equation for F{z) on the finite interval (0, 2a:o). It is uniquely solvable since its 
kernel is a continuous function. One can put x = in H5.4.4|) and the kernel A{0, v ~ z) is a continuous 
function of v and z, and the right-hand side of (|5.4.4|l at x = is a continuous function of z. Thus F(z) 
is uniquely recovered for all 2; > from A{x, y), y > x > 0. Step S ^ F is done. 

3. The step A^ q is done by formula H1.5.12|l . The converse step q ^ A is done by solving Volterra 
equation p. 2. 28(1 . or, equivalently, the Goursat problem H5.1.1|l - 1(5.1.3(1 . 

We have proved: 

Theorem 5.4.1. If q € Li^i and S are the corresponding data ((1.2.17(1 . then each step in ((1.5.1U(I is 

invertible. In particular, the potential obtained by the procedure ((1.5.10() equals to the original potential 

q- 

Remark 5.4.2. If q E Li^i and Aq := Aq{x,y) is the solution to ((1.2.28(1 . then Aq satisfies equation 
((1.5.13(1 and, by the uniqueness of its solution, Aq — A, where A is the function obtained by the scheme 
((1.5.10(1 . Therefore, the q obtained by (11.5.10(1 equals to the original q. 

Remark 5.4.3. One can verify directly that the solution A{x, y) to ((1.5.13(1 solves the Goursat problem 
((5.1.1(1 - ((5.1.3(1 . This is done as in Section \4-.4\ Step 3. Therefore q{x), obtained by the scheme ((1.5.10(1 . 
generates the same A{x, y) which was obtained at the second step of this scheme, and therefore this q 
generates the original scattering data. 

Remark 5.4.4. The uniqueness Theorem \ 5.2.1\ does not imply that if one starts with a q^ = Li i, 
computes the corresponding scattering data ((1.2.17(1 . and applies inversion scheme ((1.5.10(1 . then the q 
is obtained by this scheme is equal to qo. Logically it is possible that this q generates data Si which 
generate by the ■sc/teme ((1.5.10(l potential qi, etc. To close this loop one has to check that q — q^. This is 
done in Theorem \5.4.1\ because qo = —2 '^"^^^'^-* = q{x). 



5.5 Characterization of the scattering data 

In this Section we give a necessary and sufficient condition for the data ((1.2.17(1 to be the scattering data 
corresponding to g G In Section IFTI we give such conditions on S for q to be compactly supported, 
or q e L2(R+). 
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Theorem 5.5.1. If q G Li,i, then the following conditions hold: 1) (|1.2.23l) ; 2) kj > 0, Sj > 0, 

i < j < J, S{k) = S{-k) = S-\k), k > 0, S{oo) = 1; 3) (l5T7jl hold. Conversely, if S satisfies 
conditions 1) -3), then S corresponds to a unique q £ ii.i- 

Proof. The necessity of conditions 1) ~ 3) has been proved in Theorem 15.3.11 Let us prove the suffi- 
ciency. If conditions 1) - 3) hold, then the scheme H1.5.1U|I yields a unique potential, as was proved in 
Remark 15.4.21 Indeed, equation (|1.5.13|) is of Fredholm type in L^{Rx) for every a; > if F satisfies 
(|5.3.7() . Moreover, equation (|5.3.9|l has only the trivial solution if conditions 1) - 3) hold. Every solution 
to (|5.3.9|l in L^iW^) is also a solution in ^^(Rj,) and in L°°(M^), and the proof of the uniqueness of the 
solution to (|5.3.9|l under the conditions 1) - 3) goes as in Theorem 15.3. II The role of f{k) is played by 
the unique solution of the Riemann problem: 

f+{k) = S{~k)f^{k) (5.5.1) 

which consists of finding two functions f+{k) and f_{k), satisfying (5.5.1), such that /+ is an analytic 
function in C_|_, f^{ikj) — 0, f^{ikj) ^ 0, I < j < J , /+(oo) = 1, and f^{k) is an analytic function in 
C_ such that f-{-ikj) = 0, f-{-ikj) ^0,l<j<J, /-(oo) = 1, and /+(0) = if indS'(fc) = -2J- 1, 
/+(0) 7^ if indS'(fc) — —2 J. Existence of a solution to (5.5.1) follows from the non-negativity of 
ind5(— fc) = — indS'(fc). Uniqueness of the solution to the above problem is proved as follows. Denote 
f+{k) := f{k) and f-{k) = f{—k). Assume that /i and /2 solve the above problem. Then (5.5.1) implies 



= TT-S ^^^^^ = /2(*fcj) = 0, h{ik,) + 0, A(*fc,) + 0, /i((X3) - /2((X3) - 1. 

(5.5.2) 

The function j.^ is analytic in C+ and tends to 1 at infinity in C+, The function is analytic in 

C_ and tends to 1 at infinity in C_. Both functions agree on R. Thus is analytic in C and tends to 1 
at infinity. Therefore /i(fc) = f2{k). To complete the proof we need to check that q, obtained by H1.5.10|l . 
belongs to Li^i. In other words, that q = — 2^^^^gj^ S To prove this, use H5.3.19|) and H5.3.2U|I . It 

is sufficient to check that F'(2x) € Li^i, A{x, x)F{2x) G and J^[Ax{x, s) — Ax{x, s)]F{s + x)ds e 
Li^i. The first inclusion follows from ||a;-F'||Li(R+) < oo- Let us prove that limj,^oo [2^-^(2;)] = 0. One 
has Jq sF'ds = xF{x) - Fds. Because xF' £ L'^0^+) and F e L^(M+) it follows that the limit 
Co := limx^ooxF exists. This limit has to be zero: if -F = ^ -I- o(i) as a; 00, and cq 7^ 0, then 
F ^ i^(R+). Now x\F{2x)A{x,x)\dx < c \A{x,x)\dx < 00. The last inequality follows from 
(|5.3.13|l : since F{2x) G i^(M+) it is sufficient to check that /^°° A{x, s)F{s + x)ds e L^{R+). One has 
dxj^ \A{x, s)\\F{s + x)\ds < dxcjF{2x) \A{x, s)m/< c. Here 

crF(a;) := sup|F(?/)|, ap £ L\R+). (5.5.3) 

y>x 

Note that limj^^oo x(Jf{x) = because <7f{x) is monotonically decreasing and belongs to L^(R+). 

□ 

5.6 A new equation of Mar chenko- type 

The basic result of this Section is: 
Theorem 5.6.1. Equation 

/oo 
A{t)F{t + y)dt^A{^y), -00 < y < 00, (5.6.1) 
-00 

holds, where A{y) :~ A{0^y), A{y) = for y < 0, A{x^y) is defined in H1.2.25|l and F(x) is defined in 
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Proof. Take the Fourier transform of H5.6.1|l in the sense of distributions and get: 



F{0 + AiO + A{-OF{0 = A{-0, (5.6.2) 

where, by (|1.5.11|l . 

J 

F(0 = + 2^^5,5(6 (5.6.3) 

Use p.5.16|l . the equation S{^)f{^) = /(-O, add 1 to both sides of ifO^ . and get: 

f{0 + f{~OF{0^f{~0- (5.6.4) 

From H5.6.3|) and H5.6.4|) one gets: 

,7 J 
/(?) = /(-O [5(-e) - 2^ ^ s,5{^ + tk,)] =f{0~2nY, + ^k,)f{-0 ^ fiO, (5-6.5) 

where the equation + ikj)J{—£) = was used. This equation holds because f{ikj) = 0, and the 
product (5(5 + ikj)f{~C) makes sense because /(5) is analytic in C+. Equation (|5.6.5|l holds obviously, 
and since each of our steps was invertible, equation (|5.6.1|l holds. □ 

Remark 5.6.2. Equation (|5.6.1|) has a unique solution A{y), such that A{y) G iy^(M-i-) and A{y) = 
for y <0. 

Proof. Equation 15.6.1|) for y > is identical with p.5.13|) because A{—y) = for y > 0. Equation 
(|1.5.13|l has a solution in L^(R+) and this solution is unique, see Theorem 15. 3. II Thus, equation (|5. 6.1(1 
cannot have more than one solution, because every solution A{y) e L^(M+), A{y) = for y < 0, of 
(|5.6.1|) solves (11.5. 13|l . and (|1.5.13|) has no more than one solution. On the other hand, the solution 
A{y) e iy^(M+) of (|1.5.13|) does exist, is unique, and solves H5.6.1f) . as was shown in the proof of 
Theorem l5.6.1l This proves Remark 15.6.21 □ 

5.7 Inequalities for the transformation operators and applica- 
tions 

5.7.1 Inequalities for A and F 

The scattering data (1.2.17) satisfy the following conditions: 

A) kj,sj > 0, S{-k) = 5(fc) = S'-i(/c), A: G R, S'(oo) = 1, 

B) := indS{k) := ^ dlogS{k) is a nonpositive integer, 

C) F e LP, p = 1 and p = oo, xF' e L\ LP := LP{0, oo). 

If one wants to study the characteristic properties of the scattering data, that is, a necessary and 
sufficient condition on these data to guarantee that the corresponding potential belongs to a prescribed 
functional class, then conditions A) and B) are always necessary for a real- valued q to be in Li,i, the 
usual class in the scattering theory, or in some other class for which the scattering theory is constructed, 
and a condition of the type C) determines actually the class of potentials q. Conditions A) and B) are 
consequences of the selfadjointness of the Hamiltonian, finiteness of its negative spectrum, and of the 
unitarity of the S — matrix. Our aim is to derive some inequalities for F and A from equation l(1.5.13|l . 
This allows one to describe the set of q, defined by (1.5.13). 

Let us assume: 

snp\Fiy)\:=aFix)eL\ F' e (5.7.1) 

y>x 
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The function ap is monotone decreasing, [^(a::)! < apix). Equation (|1.5.13|l is of Fredholm type in 
LP. := L'P{x, oo) Vx > and p — 1. The norm of the operator F :— in H1.5.13|l can be estimated : 

\\F,\\< aF{x + y)dy<aiF{2x), a^pix) := <TF{y)dy. (5.7.2) 

Therefore H1.5.13|) is uniquely solvable in for any x > xq if 

aiF{2xo) < 1. (5.7.3) 

This conclusion is vahd for any F satisfying H5.7.3|) . and conditions A), B), and C) are not used. Assuming 
H5.7.3|l and H5.7.1|) and taking x > xq, let us derive inequalities for A = A{x,y). Define 

aA{x) := sup \A{x,y)\ := \\A\\. 

y>x 

From H1.5.13|l one gets: 

POO 

crA{x) < (Tf {2x) + a a{x) snp / (Tf{s + y)ds < (TFi^x) + aAix)aiFi'2,x). 

y>x J X 

Thus, if if^T^ holds, then 

(^Aix) < c(Ti?(2x), x > xq. (5.7.4) 
By c > different constants depending on xq are denoted. Let 

/>oo 

<jia{x) := ||A||i / \A{x,s)\ds. 

J X 

Then 1)1.5. 13|l yields 0-1^(2;) < o-i_F(2a;) + aiA{x)aiF{2x). So 

<^ia{x) < caiF{2x), X > Xq. (5.7.5) 
Differentiate (|1.5.13|) with respect to x and y and get: 

{I + F^)A,{x,y)^A{x,x)F{x + y)-F'{x + y), y>x>0, (5.7.6) 

and 



Ay{x,y)+ A{x,s)F'{s + y)ds::^-F\x + y), y>x>0. (5.7.7) 

J X 

Denote 

/"OO 

a2F{x) := / \F\y)\dy, a2F{x) G L\ (5.7.8) 

J X 

Then, using (|5.7.7|) and (|5.7.4|) . one gets 
\\Ay\\i< / \F'{x + y)\dy + a,A{x)snp \F'{s + y)\dy < a2F{2x)[l + caM2x)] < caM^x), (5.7.9) 



and using 1)5.7.6(1 one gets: 

\\Ax\\i < A{x,x)aiF{'2x)+a2F{2x) + \\A I |icrii?(2a;), 

so 

\\A,\\i<c[a2F{2x)+aiF{2x)aF{2x)]. (5.7.10) 
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Let y = X in (|1.5.13|l . then differentiate (|1.5.13|l with respect to x and get: 



a;) = -2F\2x) + A{x, x)F{2x) 



Ax{x, s)F{x + s)ds ■ 



A{x,s)F'{s + x)ds. (5.7.11) 



From 15.7.4|) . H5.7.5|) . H5.7.1U|I and H5.7.11|l one gets: 

\A{x,x)\ < 2|F'(2x)| + cal{2x) + caF{2x)[a2F{2x) + aiFi2x)aFi'2x)] + caF{2x)a2F{2x). (5.7.12) 
Thus, 

x\A{x,x)\eL\ (5.7.13) 

provided that xF'{2x) E L^, x(jp{2x) e L^, and xaF{2x)a2F{'2'X) E L^. Assumption (|5.7.1|l implies 
xF'{2x) € . If aF{2x) G L^, and (Tf{2x) > decreases monotonically, then xaF{x) — ^ as a; ^ oo. 
Thus xap(2x) € L^, and cr2_F(2x) S because dx \F'{y)\dy — \F'{y)\ydy < oo, due to 
(EZHl- Thus, inZH) jniplics (inZH), (inZSJ, inZHl, (EZII), and (EZISi, while (|5.7.1:^|) and (1.5.13) 
imply q G Li^i where Li^i — {q : q — q, x\q{x)\dx < oo}, and xq > satisfies (|5.7.3|) . 

Let us assume now that (|5.7.4|) . H5.7.5() . H5.7.9I) . and H5.7.10|l hold, where ctf G and a2F & are 
some positive monotone decaying functions (which have nothing to do now with the function F, solving 
equation (|1.5.13(l . and derive estimates for this fmiction F. Let us rewrite H1.5.13|l as: 



(5.7.14) 



F{x + y)+ A{x,s)F{s + y)ds = -A{x,y), y>x>0. 

J X 

Let x + y = z, s + y — V. Then, 

/■oo 

F{z)+ A{x,v + X - z)F(v)dv = -A{x,z - x), z>2x. 

J z 

From H5.7.15|l one gets: 

POO 

aF{2x) < aA{x) + crF{2x) sup / \A{x,v + x — z)\dv < aA{x) + aF{2x) 

z>2x J z 

Thus, using H5.7.5() and H5.7.3() . one obtains: 

aF{2x) < caAix). 

Also from (|5.7.15|l it follows that: 

a^Fi2x) ■.= \\F\\,:^£\Fiv)\dv 

< \Aix, z - x)\dz + J^^ \A{x, v + x- z)\\F{v)\dvdz 

<\\A\\, + \\FMA\\„ 

so 

(JifC^x) < caiA{x). 

From 15.7.6(1 one gets: 

\F'{x + y)\dy = (T2f{'2.x) < caA{x)aiAix) + \\A.j:\\ + c\\A^\\iaiAix) . 
Let us summarize the results: 



(5.7.15) 



(5.7.16) 



(5.7.17) 



(5.7.18) 
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Theorem 5.7.1. If x > xo and (|5.7.1|l hold, then one has: 

aA{x) < caF{2x), aiA{x) < caipi^x), ||Aj^||i < cr2_F(2a;)(l + ccriF(2a;)), 
ll^xlli < c[(T2F(2a;)+c7iF(2x)aF(2x)]. 



(5.7.19) 



Conversely, if x > xq and 

aAix)+aiA{x) + \\A^\\i + \\Ay\\i < oo, (5.7.20) 

then 

ct_f(2x) < caAix), (Tif(2x) < caiAix), 
<y2F{x) < c[aA{x)aiA{x) + \\Aj;\\i{l + aiA{x))]. 



(5.7.21) 



In the next section we replace the assumption x > xq >Obya;>0. The argument in this case is 
based on the Fredholm alternative. 

5.7.2 Characterization of the scattering data revisited 

First, let us give necessary and sufficient conditions on S for q to be in Li i. These conditions are known 
(|M], 10, |R9j . Section 5.5) but we give a short new argument. We assume throughout that conditions 
A), B), and C) hold. These conditions are known to be necessary for q G Indeed, conditions A) 

and B) are obvious, and C) is proved in Theorem 15 . 7 . II and Theorem l5.7.4l Conditions A), B), and C) 
are also sufficient for q G Indeed if they hold, then we prove that equation l|1.5.13|l has a unique 

solution in L]. for all x > 0. This was proved in Theorem 5.3.1, but we give another proof. 

Theorem 5.7.2. If A), B), and C) hold, then (|1.5.13|l has a solution in L\ for any x > and this 
solution is unique. 

Proof. Since is compact in L\, \/x > 0, by the Fredholm alternative it is sufficient to prove that 

{I + F^)h = 0, heLl, (5.7.22) 

implies h = 0. Let us prove it for x = 0. The proof is similar for a; > 0. li h G L^, then h e L°° because 
\Moo < IMLiaFiO). lfh€ L^nL°^\ then h e because] |/i| < Thus, if /i G and 

solves (EZpi, then he L^nL^nL°°. 

Denote h = h{x)e'^^'^ dx, h G L^. Then, 

/>oo 

h'^dk = 0. (5.7.23) 
Since F(x) is real-valued, one can assume h real-valued. One has, using Parseval's equation: 







{{I + Fo)h,h)^—\\h\\^ + — [1- S{k)]h^{k)dk + y" sjh^ hj-.^ e~''^''h{x)dx 
2tt 2n J_^ Jo 



Thus, using H5.7.23|l . one gets 

hj = 0, 1 < J < J, (h,h) = {S{k)h, h{~k)), 

where we have used real-valuedness of h, i.e. h{—k) = h{k),\/k G R. 

Thus, (h,h) ^ (h,S{-k)h{-k)), where A) was used. Since ||5(-fc)|| = 1, one has \\h\\^ ^ \(h, 
S{-k)h{-~k))\ < so the equality sign is attained in the Cauchy inequality. Therefore, h(k) = 

S{-k)h{^k). 
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By condition B), the theory of Riemann problem guarantees existence and uniqueness of an analytic 
in C+ := {k : Irak > 0} function f{k) /+(fc), f{ikj) = 0, f{ikj) ^ 0, 1 < j < J, /(oo) = 1, such that 

f+{k) ^ S{-k)f^{k), keR, (5.7.24) 

and /_(fc) = /(-fc) is analytic in C_ := {k : Imk < 0}, /-(oo) = 1 in C_, f-{-ikj) = 0, f-{-ikj) 7^ 0. 
Here the property S{-k) = S~^{k), Vfc G R is used. 
One has _ _ 

h{k) _ h{-k) 



Vj(fc) := = fceM, hj-h{ikj)^0, l<j<J. 



The function ip{k) is analytic in C+ and ipi^k) is analytic in C_, they agree on E, so ip{k) is analytic 
in C. Since /(oo) — 1 and h{oo) = 0, it follows that ■0 = 0. 

Thus, h = and, consequently, h(x) = 0, as claimed. Theorem l5.7.2l is proved. □ 

The unique solution to equation (|1.5.14() satisfies the estimates given in Theorem 15. 7. II In the proof 
of Theorem 15. 7. li the estimate x|^(a;,a;)| S L^(a;o,oo) was established. So, by (1.5.13), xq G L^{xo,oo). 

The method developed in Section 5.7.1 gives accurate information about the behavior of q near 
infinity. An immediate consequence of Theorem 15.7.11 and Theorem 15. 7. 21 is: 

Theorem 5.7.3. If A). B), and C) hold, then q, obtained by the scheme (11.5.10(1 belongs to ii.i(xo,oo). 

Investigation of the behavior of q{x) on {0,xo) requires additional argument. Instead of using the 
contraction mapping principle and inequalities, one has to use the Fredholm theorem, which says that 
||(/ + F2;)~^|| < c for any a; > 0, where the operator norm is taken for acting in L^, p = 1 andp — 00, 
and the constant c does not depend on a; > 0. 

Such an analysis yields: 

Theorem 5.7.4. // and only if A), B), and C) hold, then q e £1,1. 

Proof. It is sufficient to check that Theorem 15.7.11 holds with x > replacing x > xq. To get (I5.7.4|) 
with Xq — 0, one uses (|1.5.14f) and the estimate: 

\\A{x,y)\\ < \\{I + F,y'\\\\F{x + y)\\<caF{2x), 1 1 • 1 1 = sup | • U > 0, (5.7.25) 

y>x 

where the constant c > does not depend on x. Similarly: 

/>oo 

||^(a;,2/)||i < csup / \F{s + y)\dy < c<jif{2x), x>0. (5.7.26) 

S>X J X 



From 1(5. 7. 6|) one gets: 



\\AAx,y)\\i<c[\\F'{x + y)\\i+A{x,x)\\F{x + ym 
< ca2F{2x) + caF{2x)aiF{2x), a; > 0. 



(5.7.27) 



From H5.7.7|l one gets: 

\\Ay{x,y)\\i < c[a2Fi2x) + aiF{2x)a2Fi2x)] <a2F{2x). (5.7.28) 

Similarly, from H5.7.11|) and 1(5. 7. 24(1 - H5.7.27|l one gets (|5.7.12() . Then one checks (|5.7.13|) as in the proof 
of Theorem 15. 7. II Consequently Theorem 15.7.11 holds with a;o = 0. Theorem l5.7.4l is proved. □ 
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5.7.3 Compactly supported potentials 

In this Section necessary and sufficient conditions are given for q E L'^ ^ := {q : q = q,q = if x > 
a, Jq x\q\dx < oo}. Recall that the Jost solution is: 

/•OO 

fix, k) = e^'^-^ + / A{x, y)e'''ydy, /(O, k) /(fc). (5.7.29) 

Lemma 5.7.5. If q € i, then f{x, k) — e^'^^ for x > a, A{x, y) = for y >x> a, F{x + y) = for 
y>x>a (cf (|1.5.13(l ), and F{x) = for x > 2a. 

Thus, H1.5.13|l with x = yields A{0, y) A{y) = for a; > 2a. The Jost function 

f{k) = 1+ r A{y)e'^ydy, A{y) e W^^\Q, a), (5.7.30) 



is an entire function of exponential type < 2a, that is, |/(fc)| < ce^"''"'!, k € C, and S{k) = f{—k)/f{k) 
IS a meromor phic function in C. In jS.V.aOl) W^^p is the Sobolev space, and the inclusion 15.7.30() follows 
from Theorem l5.7.1l 

Let us formulate the assumption D): 

D) the Jost function f{k) is an entire function of exponential type < 2a. 

Theorem 5.7.6. Assume A), B), C) and D). Then q e L^ ^. Conversely, if q & Ll ^, then A),B), C) 
and D) hold. 

Necessity. If g g then A), B) and C) hold by Theorem 15.7.41 and D) is proved in Lemma fS. 7 .51 

The necessity is proved. 

Sufficiency. If A), B) and C) hold, then q e Li i. One has to prove that q = Q for x > a. If D) holds, 
then from the proof of Lemma [5.7. 51 it follows that A{y) = for ?/ > 2a. 
We claim that F(x) = for x > 2a. 

If this is proved, then (|1.5.13|) yields A{x, y) = for y > x > a, and so q = for a; > a by (1.5.13). 

Let us prove the claim. 

Take a; > 2a in (1.5.12). The function 1 — S{k) is analytic in C+ except for J simple poles at the 
points ikj. If x > 2a then one can use the Jordan lemma and the residue theorem and get: 



1 ^ f(-ik \ 

Fs{x)^— [l-5(fc)]e''=-& = -z^A— ^e-^^-, x > 2a. (5.7.31) 
27rJ-oo t-^ f(lki) 



1 on..^.^k.,„,_ .y-/H^. 

U 

Since /(fc) is entire, the Wronskian formula 

/'(0,fc)/(-fc)-/'(0,-fc)/(fc) = 2zfc 
is valid on C, and at fc = ikj it yields: 

f'{0,ik,)f{~ik,)^~2k„ 
because f{ikj) = 0. This and (|5.7.31|) yield 

Oh 

Fsix) = y 6^^=^^ = -Ys^e-"^^ = -Fd{x), X > 2a. 

j^^nO,^kj)f{^k,) 

Thus, F{x) = Fs{x) + Fd{x) — for x > 2a. The sufficiency is proved. 

Theorem l5.7.6l is proved. □ 

In 1^ a condition on S, which guarantees that q = for x > a, is given under the assumption that 
there is no discrete spectrum, that is F = F^. 
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5.7.4 Square integrable potentials 

Let us introduce conditions (|5.7.32|l - (|5.7.34|) 

2ih{^{k)-\ + ^\(^I?i^):^L^, <9:= / qds, (5.7.32) 

2ik Jo 

k{\-S{k)^%\^L', (5.7.33) 
ik 

ei'- (5.7.34) 

Theorem 5.7.7. A), B), C), and any one of the conditions (|5.7.32|l ~ (|5. 7.34(1 hold, then q e L'^{R). 

Proof. We refer to ^ for the proof. □ 
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Chapter 6 

Inverse scattering problem with 
fixed-energy phase shifts as the data 



6.1 Introduction 



In Subsection 15.1.31 the scattering problem for spherically symmetric q was formulated, see (|5.1.15|l 
H5.1.17|l . The 6i are the fixed-energy {k = const > 0) phase shifts. Define 
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where (p = feir) is a regular solution to 



2 ^ I 2 2 / \ 
r -7^ + r - r qir) 



such that 



= ui+ Kir, p)ui{p)p-^dp, K{r, 0) = 0, 



(6.1.1) 



(6.1.2) 



(6.1.3) 



and ui = ^J^Jgj^i. (r), Jiir) is the Bessel function. In (|6.1.3|l K(r, p) is the transformation kernel, I + K 
is the transformation operator. In (|6.1.2|l we assume that fc = 1 without loss of generality. The ip£ is 
uniquely defined by its behavior near the origin: 



Mr) 



(2£+ 1)!! 

For ue we will use the known formula ([GR, 8.411.8]): 



(6.1.4) 



(6.1.5) 



where r(z) is the gamma-function. 

The inverse scattering problem with fixed-energy phase shifts {5£}^=o,i,2... as the data consists of 
finding q{r) from these data. We assume throughout this chapter that q{r) is a real-valued function, 
q{r) = for r > a, 

pa 

r2|g(r)pdr < 00. (6.1.6) 



Conditions H6. 1.6(1 imply that q e L'^{Ba), Ba := {x : x & M^, |a;| < a}. 
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In the literature there are books |CS| and |N] where the so called Newton-Sabatier (NS) theory is 
presented, and many papers were published on this theory, which attempts to solve the above inverse 
scattering problem with fixed-energy phase shifts as the data. In Section 6.4 it is proved that the NS 
theory is fundamentally wrong and is not an inversion method. The main results of this Chapter are 
Theorems 6.2.2, 6.3.1, 6.5.1, and the proof of the fact that the Newton-Sabatier theory is fundamentally 
wrong in the sense that its foundations are wrong. 

6.2 Existence and uniqueness of the transformation operators 
independent of angular momentum 

The existence and uniqueness of K{r,p) in H6.1.3|) we prove by deriving a Goursat problem for it, and 
investigating this problem. Substitute (|6.1.3|l into H6.1.3f) , drop index £ for notational simplicity and get 



= — r q(r)u + {r ~r q{r)) j K{r,p)up dp 

K{r,p)p-^Lopudp + r^d^ I K{r, p)up-~'^dp. (6.2.1) 







We assume first that K{r,p) is twice continuously diffcrcntiable with respect to its variables in 
the region < r < oo, < p < r. This assumption requires extra smoothness of q{r), q{r) G 
C^(0, a). If q{r) satisfies condition f (|6.1.6|l ). then equation H6.2.7|) below has to be understood in the 
sense of distributions. Eventually we will work with an integral equation ()6.2.34|l (see below) for which 
assumption f (|6.1.6|l 'l suffices. 

Note that 

I K{r, p)p^'^Lopudp = / LopK{r, p)up^^dp + K{r,r)ur - Kp{r,r)u, (6.2.2) 
Jo Jo 

provided that 

K{r,0) = 0. (6.2.3) 



We assume (|6.2.3f) to be valid. Denote 



Then 



k (6.2.4) 

dr 



r^df. j K{r, p)up ^dp = Ku + K{r, r)ur K{r, r)u+ 











Kr{r,r)u + r / Krr{r, p)up dp. (6.2.5) 
J{ 

Combining (|6.2.1|l - H6.2.5|l and writing again Ui in place of u, one gets 



0= / [LrK{r,p)~LopKir,p)]uiip)p-'dp + Ui{r)[-r\{r)+K- 

2Kr{r,r) 







+ Kr{r,r)+Kp{r,r)], Vr > 0, £ = 0,1,2,.... (6.2.6) 
Let us prove that (|6.2.6|l implies: 

LrK{r,p)=LopK{r,p), < p < r, (6.2.7) 



55 



2K 2K{r, r) _2 d K{r, r) 
r dr 

This proof requires a lemma. 



^ir)-^-^--^^—-^. (6.2.8) 



Lemma 6.2.1. Assume that pf{p) e L^(0, r) and pA{p) e L^{0,r). If 

= / f{p)u,{p)dp + ue{r)A{r) V£ = 0, 1, 2, (6.2.9) 



then 

f{p) EE and A{r) = 0. (6.2.10) 
Proof. Equations (|6.2.9|l and H6.1.5() imply: 

= £^ dt{l ij-^' [ dppf{p)e^P' + 



Therefore 



0=1 dt ^ ' [ dppf{p)e^P'+rA{r)e^^% ;-0,l,2,... (6.2.11) 



J -1 

Recall that the Lcgendre polynomials are defined by the formula 



and they form a complete system in iy'^(— 1, 1) 
Therefore (|6.2.11() implies 



Equation 16.2.13|) implies 



and 



/ dppf{p)e'P' + rA{r)e'^'' =0 Vt e [-1, 1]. (6.2.13) 

/ dppf{p)e'P' = 0, Vie [-1,1], (6.2.14) 
Jo 

rA{r) = 0. (6.2.15) 

Therefore A{r) = 0. Also f{p) = because the left-hand side of H6.2.14|l is an entire function of t, which 
vanishes on the interval [—1, 1] and, consequently, it vanishes identically, so that pf{p) = and therefore 

f[p) ^ 0. 

Lemma 16.2.11 is proved. □ 

We prove that the problem 1)6. 2. 7|) . (|6.2.8|l . (|6.2.3|l . which is a Goursat-type problem, has a solution 
and this solution is unique in the class of functions K{r,p), which are twice continuously differentiable 
with respect to p and r, < r < oo, < p < r. 

In this section we assume that q{r) G C^(0, a). This assumption implies that K{r,p) is twice 
continuously differentiable. If (|6.1.6|l holds, then the arguments in this section which deal with integral 
equation H6.2.34|) remain valid. Specifically, existence and uniqueness of the solution to equation (|6.2.34f) 
is proved under the only assumption r\q(r)\dr < oo as far as the smoothness of q{r) is concerned. 
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By a limiting argument one can reduce the smoothness requirements on q to the condition (|6.1.6|l . 
but in this case equation (|6.2.7|l has to be understood in distributional sense. 
Let us rewrite the problem we want to study: 

r^Krr - P^Kpp + [r^ - r^q{r) - p^]K{r, p) = 0, < p < r, (6.2.16) 

K{r,r) = ^ sqis)ds gir), (6.2.17) 

K{r,0) = 0. (6.2.18) 

The difficulty in the study of this Goursat-type problem comes from the fact that the coefficients in front 
of the second derivatives of the kernel K{r, p) are variable. 

Let us reduce problem H6.2.16|l - (|6.2.18|) to the one with constant coefficients. To do this, introduce 
the new variables: 

i^ = lnr + lnp, 77 = In r — In p. (6.2.19) 



Note that 



and 



Let 



r = e^, p = e^, (6.2.20) 



77 >0, -oo<^<oo, (6.2.21) 
dr = \{di + a„), dp = ^(9? - a„). (6.2.22) 



A routine calculation transforms equations H6.2.16|l - (|6.2.18|l to the following ones: 

B(^{C,v)-lB^iC,v) + Qitv)B ^0, ?7>0, -(X3<e<oo, (6.2.23) 

B{^, 0)=g (pi) := GiO, -00 < e < 00 (6.2.24) 

B(-oo,77) = 0, 77 >0, (6.2.25) 

where g{r) is defined in (|6. 2.17(1 . 



Here we have defined 



Q(^, ,7) 1 - e«+"g (e^) - e«-"] , (6.2.26) 



and took into account that p — r implies 77 = 0, while p = implies, for any fixed ?7 > 0, that ^ = — 00. 
Note that 

sup e"^G'(0 < c, (6.2.27) 

sup / \Q{s,rf)\ds <c{A,B), (6.2.28) 

for any A G K and B > 0, where c{A, i?) > is a constant. 

To get rid of the second term on the left-hand side of H6.2.23|l . let us introduce the new kernel L(^, 77) 
by the formula: 

Li^,r,):=Bi^,7j)e'^. (6.2.29) 
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Then (|6.2.23|) - 16.2.251) can be written as: 



—oo < ^ < oo, 



(6.2.30) 



L(e,0) = e-3G(O:-6(e):=^^ -oo < ^ < oo, (6.2.31) 

L(_oo,77) = 0, 77 >0. (6.2.32) 

We want to prove existence and uniqueness of the solution to (|6.2.30|l - 16.2.32|) . In order to choose a 
convenient Banach space in which to work, let us transform problem (|6.2.30|) - H6.2.32|) to an equivalent 
Volterra-type integral equation. 

Integrate H6.2.3U|I with respect to 77 from to ry and use (|6.2.31|) to get 

Ld^, V) - + /' Q(e, t)L{^, t)dt = 0. (6.2.33) 

Integrate ()6.2.33|l with respect to ^ from — 00 to ^ and use H6.2.33|l to get 

L{t 7]) = - ds r dtQ{s, t)L{s, t) + 6(0 := VL + &, (6.2.34) 

J-OD Jo 

where ^ 

VL :=- [ ds [ dtQ{s,t)L{s,t). (6.2.35) 



Consider the space X of continuous functions L{^,ri), defined in the half-plane rj > 0, — c»<i^<oo, 
such that for any B > and any — c» < ^ < c» one has 

||L|| := ||L|Ub := sup (e-^*|i(s, i)|) < 00, (6.2.36) 

— oc<s<A 
0<t<B 

where 7 > is a number which will be chosen later so that the operator V in H6.2.34(l will be a contraction 
mapping on the Banach space of functions with norm H6.2.36|l for a fixed pair A, B. To choose 7 > 0, 
let us estimate the norm of V . One has: 

\\VL\\< sup (t ds f dt\Q{s,t)\e~"''^'^-''^e-^'\L{s,t)\\ 

-ao<^<Afl<r,<B \J -oo Jo / 



sup rfs /"d^(2e^+*+e^+*|(7fe^) 1) e-T(''-*) < (6.2.37) 

<A,0<J)<S J-00 Jo ^ \ J J 7 



< Hill 

-ca<^<A,0<ri<B 



where c > is a constant depending on A, B and r\q{r)\dr. Indeed, one has: 

ds r dte^+*-^("-*) = 2e^ /" dte'-^'^^-'Ut < 2e^+^i^^^ = -, (6.2.38) 
) Jo Jo 7 7 
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and, using the substitution cr = e , one gets: 



ds / dte"+'\q(e' 

-oc Jo 



/ dte-''^"-*) / d(jG\q{a)\ = 
Jo Jo 

2(1 _e-^^) r ^ , , C2 



= 2 

Jo Jo 

da(j\q{(j)\ := — . (6.2.39) 

1 Jo 7 

From these estimates inequahty (6.2.38) follows. 

It follows from (6.2.38) that y is a contraction mapping in the space Xab of continuous functions 
in the region — oo <^<^, < V ^ with the norm H6.2.36|) provided that 

7 > c. (6.2.40) 

Therefore equation (|6.2.34|l has a unique solution L{^,rj) in the region 

~oo<i<A, 0<?]<B (6.2.41) 

for any real A and i? > if l|6.2.4()|l holds. This means that the above solution is defined for any ^ G M 
and any 77 > 0. 

Equation (|6. 2.34(1 is equivalent to problem H6.2.3U|I - (|6.2.32|l and, by ()6.2.29|l . one has: 

S(e,^)=i(e,'?)ei (6.2.42) 

Therefore we have proved the existence and uniqueness of i?(^, rj), that is, of the kernel K{r, p) = B{^, rj) 
of the transformation operator IjC. 1.3(1 . Recall that r and p are related to ^ and rj by formulas 1(6.2.20(1 . 
Let us formulate the result: 

Theorem 6.2.2. The kernel K{r,p) of the transformation operator (16.1.3(1 solves problem (16.2.16(1 
- ((6.2.181) . The solution to this problem does exist and is unique in the class of twice continuously 
differentiable functions for any potential qir) £ C^(0, a). If q{r) G L°°{0,a), then K(r,p) has first 
derivatives which are bounded and equation ((6.2.16(1 has to be understood in the sense of distributions. 
The following estimate holds for any r > 0: 

\K{r, p)\p-'^ dp < 00. (6.2.43) 

Proof of Theorem \6.2./^ We have already proved all the assertions of Theorem 16.2.21 except for the 
estimate ((6.2.43() . Let us prove this estimate. 
Note that 

\K{r,p)\p~^dp = r \L{2lnr - r], f])\e~i dr] < 00 (6.2.44) 



Indeed, if r > is fixed, then, by ((6.2.20(1 . ^ + 77 = 21n7' = const. Therefore = —drj, and p ^dp 
i(d^ — dri) = —drj, ^ = 21nr — 77. Thus: 



\K{r,p)\p^^dp^ I \L{2lnr - r],Tj)\e^'"^ ^ drj = r I |L(2 Inr - 77, 7;)|e"*fi7y. (6.2.45) 
Jq Jo 

The following estimate holds: 

\L{^,V)\ < ce(2+^i)['"'i (6.2.46) 
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where ej > 0, j ~ 1,2, are arbitrarily small numbers and /ii is defined in formula (|6.2.52|) below, see 
also formula (|6.2.49l) for the definition of /i. 

Estimate (|6.2.4t)|) is proved below, in Theorem 16. 2. 21 

From (|6.2.45|) and estimate (|5.2.56|l (see below) estimate (16.2. 43|) follows. Indeed, denote by / the 
integral on the right-hand side of H6.2.45|) . Then, by HB.2.5B|I one gets: 

/<2 + 2y [^A^i(^l°g^)]" ^2exp[2/ii(21ogr)] <oo. (6.2.47) 
n! 

Theorem 16. 2. 21 is proved. □ 
Theorem 6.2.3. Estimate (I6.2.46|) holds. 
Proof of Theorem \ti.2.'J[ From H6.2.34|) one gets: 

m(C,7?) <co + (W^to)(^,77), m(e,77) :=|i(e,r;)|, (6.2.48) 
where cq = sup_^<^.<^ |6(C)| < \ s\q{s)\ds (see (|6.2.31ll'). and 



Wm:=l' ds I dtn{s + t)m{s,t), fi{s) := ^e" {l + \q{ei)\) . (6.2.49) 



It is sufficient to consider inequality (|6.2.48|l with Cq = 1: if cq = 1 and the solution mo(^, rj) to (|6.2.48|) 
satisfies (|6.2.46|l with c = ci, then the solution to(^,7/) of H6.2.48|) with any cq > satisfies (|6.2.46ll with 

C = CoCi. 

Therefore, assume that cq = 1, then H6.2.48|l reduces to: 

m{£_,7j) <l + {Wm){^,r]). (6.2.50) 

Inequality (|6.2.46|l follows from H6.2.50|l by iterations. Let us give the details. 
Note that 

ds dt^i{s + t)^ dt ds^l{s + t)= diA'i(C + < Wi(C + ??)• (6.2.51) 

-oo Jo Jo J-oo Jo 

Here we have used the notation 



= / Ks)ds, (6.2.52) 

J —oc 

and the fact that /ii(s) is a monotonically increasing function, since /Lt(s) > 0. Note also that /Lti(s) < oo 
for any s, — oo < s < oo. 
Furthermore, 



W^l< I ds I dtfi{s + t)tfj.i{s + t) < dtt ds/x(s + t)/xi(s + t) = ^. 



(6.2.53) 



J —OO Jo Jo J —oo 

Let us prove by induction that 

nl n 

For n = 1 and n = 2 we have checked (|6.2.54|l . Suppose H6.2.54|l holds for some n, then 



< :i^^liAl_L^. (6.2.54) 



W"+'l < W ' 'M = / dt— dsfiis + t) "^'' , ' < . ' ,v ■ (6.2.55) 

\ „i „i ; / „\ I ' n\ ~ (n + 1)! (n + l)! ^ ' 
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By induction, estimate (|6.2.53|) is proved for all n = 1,2, 3, .... Therefore 16.2.5U|) implies 

^)<l + y !? Mi±^ < c^i2+^.Mv+0r^^'^ ^ (6.2.56) 
^-^ TV. n\ 



where we have used Theorem 2 from [Lev, section 1.2], namely the order of the entire function F{z) :— 
1 + Xl^i J^vp- is ^ and its type is 2. The constant c > in H6.2.46() depends on ej, j = 1, 2. 

Recall that the order of an entire function F{z) is the number p := limsup^^o^ '"'"^^ , where 
Mi?(r) := max\z\=r\F{z)\. The type of F{z) is the number a := limsup^^^g^ inMF(r) ^ known [Lev], 
that if F(z) = X!J^o ^n^^ is an entire function, then its order p and type a can be calculated by the 
formulas: 

,. '"'Inn limsup„^^(n|c„|") fROKv\ 
p = limsup- —, a = . (6.2.57) 

71 — ^CJO I^Tl — j- Gp 

If Cn — -^7^, then the above formulas yield p — and a = 2. Theorem 16. 2. 31 is proved. □ 

6.3 Uniqueness theorem. 

Denote by £ any fixed subset of the set N of integers {0, 1,2,...} with the property: 

5^1=00 (6.3.1) 

eec 

Theorem 6.3.1. f[R10]j Assume that q satisfies H6.1.6f) and ()6.3.1|l holds. Then the data {Si}\/£^c 
determine q uniquely. 

The idea of the proof is based on property C-type argument. 

Step 1 : If qi and q2 generate the same data {Si}\/i,£Ci then the following orthogonality relation holds 
for p:= 91-92: 

hi£) := [ p{r)^u{r)(j}2i{r)dr = "ileC, (6.3.2) 

where 4)ji is the scattering solution corresponding to qj, j — \, 2. 

Step 2 : Define /ii(£) := 2'^^[T{1 + l)Yh{l), where T is the Gamma-function. Check that hi{l) is 
holomorphic in 0+ :— {£ : Re£ > 0}, £ = <t + ir, cr > 0, and r are real numbers, hi{£) £ N (where 
N is the Ncvanlinna class in 11+ ), that is 



r , 1 - re'-'f 
sup / log+|/ii(— — —)\dip<(X), 



, I logx if log a; > 0, 

where log x = < If /ii e TV vanishes G C, then hi ~ in 11+, and, by property 

I if logx < 0. 

C^, p{r) — 0. Theorem 6.3.1 is proved. □ 

6.4 Why is the Newton-Sabatier (NS) procedure fundamentally 
wrong? 

The NS procedure is described in jN^ and (CS) . A vast bibliography of this topic is given in jCSj and 

El- 
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Below two cases are discussed. The first case deals with the inverse scattering problem with fixed- 
energy phase shifts as the data. This problem is understood as follows: an unknown spherically symmet- 
ric potential q from an a priori fixed class, say ii.i, a standard scattering class, generates fixed-energy 
phase shifts Si, I — 0,1,2, ... ,. The inverse scattering problem consists of recovery of q from these data. 

The second case deals with a different problem: given some numbers Si, I = 0, 1, 2, ... ,, which are 
assumed to be fixed-energy phase shifts of some potential q, from a class not specified, find some potential 
qi, which generates fixed-energy phase shifts equal to 5i, 1 = 0,1,2,...,. This potential qi may have no 
physical interest because of its non-physical" behavior at infinity or other undesirable properties. 

We first discuss NS procedure assuming that it is intended to solve the inverse scattering problem in 
case 1. Then we discuss NS procedure assuming that it is intended to solve the problem in case 2. 

Discussion of case 1: 

In |JN2| and ]lNj a procedure was proposed by R. Newton for inverting fixed-energy phase shifts 
Si, I = 0,1,2, ... , corresponding to an unknown spherically symmetric potential q{r). R. Newton did not 
specify the class of potentials for which he tried to develop an inversion theory and did not formulate 
and proved any results which would justify the inversion procedure he proposed (NS procedure). His 
arguments are based on the following claim, which is implicit in his works, but crucial for the validity 
of NS procedure: 



Claim Nl: The basic integral equation 

r 

Kir A f(i ^ 



K{r,s) = f[r,s)- I K{r,t)f{t,s)%, < s < r < oo, (6.4.1) 





is uniquely solvable for all r > 0. 
Here 

oo I 

f{r,s) ■.= '^ciui{r)ui{s), u; := J y J,+ i (r), (6.4.2) 

ci are real numbers, the energy is fixed: fc = 1 is taken without loss of generality, J;_|_i(r) are the 
Bessel functions. If equation IjG. 4.1(1 is uniquely solvable for all r > 0, then the potential qi, that NS 
procedure yields, is defined by the formula: 

2d K(r,r) 

qi{r) = —, 6.4.3 

r ar r 

The R. Newton's ansatz (|6.4.1|) - (|6.4.2|) for the transformation kernel K{r, s) of the Schroedinger oper- 
ator, corresponding to some q{r), namely, that K{r, s) is the unique solution to (|6.4.1|l - (|6.4.2|l . is not 
correct for a generic potential, as follows from our argument below (see the justification of Conclusions). 

// for some r > equation (|6.4.1|l is not uniquely solvable, then NS procedure breaks down: it leads 
to locally non-integrable potentials for which the scattering theory is, in general, not available (see \R9f 
for a proof of the above statement) . 

In the original paper |N2| and in his book IN] R. Newton did not study the question, fundamental 
for any inversion theory: does the reconstructed potential qi generate the data from which it was 
reconstructed? 

In |CS[ p. 205], there are two claims: 

Claim i) that qi{r) generates the original shifts {Si} "provided that {5i} are not "exceptional"", and 
Claim ii) that NS procedure "yields one (only one) potential which decays faster than r^a" and 
generates the original phase shifts {Si}. 

If one considers NS procedure as a solution to inverse scattering problem of finding an unknown 
potential q from a certain class, for example q(r) e Li.i := {q : q = 9i /q r\q(r)\dr < oo}, from 
the fixed-energy phase shifts, generated by this q, then the proof, given in |("S| . of Claim i) is not 
convincing: it is not clear why the potential qi, obtained by NS procedure, has the transformation 
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operator generated by the potential corresponding to the original data, that is, to the given fixed-energy 
phase shifts. In fact, as follows from Proposition l6.4.l1 below. the potential qi cannot generate the kernel 
K{r, s) of the transformation operator corresponding to a generic original potential q{r) £ Li.i := {q : 

Claim ii) is incorrect because the original generic potential q{r) e Li^i generates the phase shifts {Si}, 
and if qi{r), the potential obtained by NS procedure and therefore not equal to q{r) by Proposition l6.4.l1 
generates the same phase shifts {Si}, then one has two different potentials q{r) and qi{r), which both 
decay faster than r~2 and both generate the original phase shifts {Si}, contrary to Claim ii). 

Our aim is to formulate and justify the following 

Conclusions: Claim Nl and ansatz H6.4.1(l - (j6.4.2|l are not proved by R. Newton and, in general, are 
wrong. Moreover, one cannot approximate with a prescribed accuracy in the norm \ \q\ \ : = r\q{r)\dr 
a generic potential q{r) e Li^i by the potentials which might possibly be obtained by the NS procedure. 
Therefore NS procedure cannot be justified even as an approximate inversion procedure. The NS procedure 
is fundamentally wrong in the sense that its foundations are wrong. 

Let us justify these conclusions: 

Claim Nl formulated above and basic for NS procedure, is wrong, in general, for the following reason: 

Given fixed-energy phase shifts, corresponding to a generic potential q G one either cannot 

carry through NS procedure because: 

a) the system (12.2.5a) in |CS| . which should determine numbers c/ in formula (|6.4.2|l . given the 
phase shifts Si, may be not solvable, or 

b) if the above system is solvable, equation (|6.4.1|l may be not (uniquely) solvable for some r > 0, 
and in this case NS procedure breaks down since it yields a potential which is not locally integrable (see 
| R9 | for a proof). 

If equation H6.4.1() is solvable for all r > and yields a potential qi by formula (|6.4.3|l . then this 
potential is not equal to the original generic potential q g as follows from Proposition 16.4.11 which 
is proved in ^9! (see also |AR,Sp : 

Proposition 6.4.1. // equation 1)6.4.1(1 is solvable for all r > and yields a potential qi by formula 
1)6.4.3(1 . then this qi is a restriction to {0,oo) of a function analytic in a neighborhood o/(0,oo). 

Since a generic potential q e is not a restriction to (0, cxd) of an analytic function, one concludes 
that even if equation ((6.4.1(1 is solvable for all r > 0, the potential qi, defined by formula ((6.4.3(1 . is not 
equal to the original generic potential q e Li^i and therefore the inverse scattering problem of finding 
an unknown q G from its fixed-energy phase shifts is not solved by NS procedure. 

The ansatz 1(6.4.1(1 - ((6.4.2(1 for the transformation kernel is, in general, incorrect, as follows also from 
Proposition 16.4.11 

Indeed, if the ansatz 1(6.4. 1() - ((6.4.2(1 would be true and formula 1(6. 4. 3() would yield the original 
generic q, that is qi = q, this would contradict Proposition 16.4.11 If formula ((6.4.3(1 would yield a qi 
which is different from the original generic q, then NS procedure does not solve the inverse scattering 
problem formulated above. Note also that it is proved in |R10) that independent of the angular momenta 
I transformation operator, corresponding to a generic q G does exist, is unique, and is defined by a 
kernel K{r, s) which cannot have representation ((6.4.2(1 . since it yields by the formula similar to ((6.4.3(1 
the original generic potential q, which is not a restriction of an analytic in a neighborhood of (0, oo) 
function to (0, c»). 

The conclusion, concerning impossibility of approximation of a generic q G Li^i by potentials qi, which 
can possibly be obtained by NS procedure, is proved in Claim I?r4.3l section 2, see proof of Claim I(r4. 31 
there. 

Thus, our conclusions are justified. □ 
Let us give some additional comments concerning NS procedure. 

Uniqueness of the solution to the inverse problem in case 1 was first proved by A. G. Ramm in 1987 
(see (R,7) ) for a class of compactly supported potentials, while R. Newton's procedure was published in 
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| N2 | , when no uniqueness results for this inverse problem were known. It is still an open problem if for 
the standard in scattering theory class of Li i potentials the uniqueness theorem for the solution of the 
above inverse scattering problem holds. 

We discuss the inverse scattering problem with fixed-energy phase shifts (as the data) for potentials 
q G because only for this class of potentials a general theorem of existence and uniqueness of 

the transformation operators, independent of the angular momenta I, has been proved, see |R10| . In 
|N2| . |N). and in jCS| this result was not formulated and proved, and it was not clear for what class 
of potentials the transformation operators, independent of I, do exist. For slowly decaying potentials 
the existence of the transformation operators, independent of I, is not established, in general, and the 
potentials, discussed in jCSj and [N] in connection with NS procedure, are slowly decaying. 

Starting with |N2| . [N). and |CSj Claim Nl was not proved or the proofs given (see \CT\ ) were 
incorrect (see |Rllj '). This equation is uniquely solvable for sufficiently small r > 0, but, in general, 
it may be not solvable for some r > 0, and if it is solvable for all r > 0, then it yields by formula 
1)6.4.3(1 a potential qi, which is not equal to the original generic potential q G as follows from 

Proposition \6.4-l\ 

Existence of "transparent" potentials is often cited in the literature. A "transparent" potential is a 
potential which is not equal to zero identically, but generates the fixed-energy shifts which are all equal 
to zero. 

In fCAl p. 207], there is a remark concerning the existence of "transparent" potentials. This remark 
is not justified because it is not proved that for the values ci, used in JCAl p. 207], equation (|6.4.1|l is 
solvable for all r > 0. If it is not solvable even for one r > 0, then NS procedure breaks down and the 
existence of transparent potentials is not established. 

In the proof, given for the existence of the "transparent" potentials in |USI p. 197], formula (12.3.5), 
is used. This formula involves a certain infinite matrix M. It is claimed in |("SI p. 197], that this matrix 
M has the property MM = I, where / is the unit matrix, and on jCSI p. 198], formula (12.3.10), it is 
claimed that a vector v exists such that Mv = 0. However, then MMv = and at the same time 
MMv = w 7^ 0, which is a contradiction. The difficulties come from the claims about infinite matrices, 
which are not formulated clearly: it is not clear in what space M, as an operator, acts, what is the 
domain of definition of M, and on what set of vectors formula (12.3.5) in | CS | holds. 

The construction of the "transparent" potential in |CS| is based on the following logic: take all the 
fixed-energy shifts equal to zero and find the corresponding c; from the infinite linear algebraic system 
(12.2.7) in jCS,; then construct the kernel f{r,s) by formula 16.4.2(1 and solve equation 1(6.4.1(1 for all 
r > 0; finally construct the "transparent" potential by formula ((6.4.3(1 . As was noted above, it is not 
proved that equation 1(6.4.1(1 with the constructed above kernel f{r, s) is solvable for all r > 0. Therefore 
the existence of the "transparent" potentials is not established. 

The physicists have been using NS procedure without questioning its validity for several decades. 
Apparently the physicists still believe that NS procedure is "an analog of the Gel'fand-Levitan method" 
for inverse scattering problem with fixed-energy phase shifts as the data. In fact, the NS procedure is 
not a valid inversion method. Since modifications of NS procedure are still used by some physicists, who 
believe that this procedure is an inversion theory, the author pointed out some questions concerning this 
procedure in [AR^ and |R9| and wrote this paper. 

This concludes the discussion of case 1. □ 

Discussion of case 2: 

Suppose now that one wants just to construct a potential qi , which generates the phase shifts corre- 
sponding to some q. 

This problem is actually not an inverse scattering problem because one does not recover an original 
potential from the scattering data, but rather wants to construct some potential which generates these 
data and may have no physical meaning. Therefore this problem is much less interesting practically 
than the inverse scattering problem. 

However, NS procedure does not solve this problem either: there is no guarantee that this procedure 
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is applicable, that is, that the steps a) and h), described in the justification of the conclusions, can be 

done, in particular, that equation (|6.4.1|l is uniquely solvable for all r > 0. 

If these steps can be done, then one needs to check that the potential qi, obtained by formula ()6.4.3|l . 

generates the original phase shifts. This was not done in 'WT and |N|. 

This concludes the discussion of case 2. □ 

The rest of the paper contains formulation and proof of Remark l6 . 4 . 21 and Claim 1^4.31 

It was mentioned in |N3j that if Q rq{r)dr ^ 0, then the numbers q in formula (|6.4.2|) cannot 

satisfy the condition X^tT l*^'! ^ This observation can be obtained also from the following 

Remark 6.4.2. For any potential q{r) € Li i such that Q := rq{r)dr ^ the basic equation (|6.4.1() 
is not solvable for some r > and any choice of ci such that J^'iZo ^ 

Since generically, for q g one has Q 7^ 0, this gives an additional illustration to the conclusion 
that equation Ht).4.1l) . in general, is not solvable for some r > 0. Conditions X^i^o 1*^'! ^ °° Q 
are incompatible. 

In CS, p. 196], a weaker condition X^z^o '"^l^zl < 00 is used, but in the examples ( [CSI pp. 189-191]), 
ci — for all I > Iq > 0, so that J^lZo I*-'! < 00 in all of these examples. 

Claim 6.4.3. The set of the potentials v{r) g Li i, which can possibly be obtained by the NS procedure, 
is not dense (in the norm \\q\\ := r\q{r)\dr) in the set Li i. 

Let us prove Remark 16.4.21 and Claim IH?4. 31 

Proof of Remark \6.4-S\ Writing H6.4.3(l as K{r,r) = — § Jq sqi{s)ds and assuming Q ^ 0, one gets the 
following relation: 

Or 

K{r, r) = --^ [1 + 0(1)] -> cx) as r cx). (6.4.4) 

If l|6.4.1|l is solvable for all r > 0, then from 16.4.2(1 and (|6. 4.1(1 it follows that K{r, s) = X^i^o (^wir) 
ui{s), where ifiiir) := ui{r) — K{r,t)ui{t)jj, so that / — i^T is a transformation operator, where K is 

the operator with kernel K{r, s), ip'/ + ipi — '^'^"^^ (fii — qi{r)ipi — 0, qi{r) is given by 16.4.3() . (pi — 0(r'+^), 
as r ^ 0, 

ui{r) sin ( r - ) , ipi{r) ~ |F/| sin ( r - + (5; ) as r ^ 00, 



where 61 are the phase shifts at fc = 1 and Fi is the Jost function at k = I. One can prove that 
sup; \Fi\ < 00. Thus, if X^i^o 1"^'! ^^"^^ 

K{r,r)^0{l) asr-^oo. (6.4.5) 

If Q 7^ then ((6.4.5(1 contradicts 1(6.4.4(1 . It follows that if Q 7^ then equation ((6.4.1(1 cannot be uniquely 
solvable for all r > 0, so that NS procedure cannot be carried through if Q 7^ and Xi^o I*-'! < 00. This 
proves Remark (6.4.21 □ 

Proof of Claim \KJ~^ Suppose that v{r) g Li^i and Qy := rv{r)dr — 0, because otherwise NS 
procedure cannot be carried through as was proved in Remark 16. 4. 21 

If Qv — 0, then there is also no guarantee that NS procedure can be carried through. However, we 
claim that if one assumes that it can be carried through, then the set of potentials, which can possibly be 
obtained by NS procedure, is not dense in Li_i in the norm :— r\q(r)\dr. In fact, any potential q 
such that Q := rq{r)dr ^ 0, and the set of such potentials is dense in cannot be approximated 
with a prescribed accuracy by the potentials which can be possibly obtained by the NS procedure. 

Let us prove this. Suppose that q G 

/>oo 

Qq :— / rq{r)dr ^ 0, and ||t;„ — gjl ^ as n — > 00, 
Jo 
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where the potentials w„ € Li^i are obtained by the NS procedure, so that Q„ :— rvn{r)dr — 0. We 
assume v„ € Li^i because otherwise w„ obviously cannot converge in the norm || • || to g e Define 
a linear bounded on Li i functional 



/(?) := r 



where ||g|| := r\q{r)\dr. The potentials v £ ii.i, which can possibly be obtained by the NS procedure, 
belong to the null-space of /, that is f{v) — 0. 

If lim„^oo \\vn - q\\ = 0, then lim„^oo \ f{q - < lim^^oo Ik - «n|| = 0. Since / is a linear 
bounded functional and /(w„) = 0, one gets: f{q — w„) — f{q) — /(w„) = /(<?)■ So if f{q) ^ then 
lim„_^oo !/(<? ~ Vn)\ — \fiq)\ 7^ 0. Therefore, no potential q e Li i with Qq ^ can be approximated 
arbitrarily accurately by a potential v(r) € Li i which can possibly be obtained by the NS procedure. 
Claim IFX^ is proved. □ 

6.5 Formula for the radius of the support of the potential in 
terms of scattering data 

The aim of this section is to prove the formula for the radius of the support of the potential in terms of 
the phase shifts. Let us make the following assumption. 

Assumption (A): the potential q{r), r = \x\, is spherically symmetric, real-valued, \q\'^dr < oo, and 
q{r) — for r > a, but q{r) ^0 on [a — e, a) for all sufficiently small e > 0. 

The number a > we call the radius of compactness of the potential, or simply the radius of the 
potential. Let A{a' , a) denote the scattering amplitude corresponding to the potential q at a fixed energy 

> 0. Without loss of generality let us take fc = 1 in what follows. By a' ,a E 5^ the unit vectors in 
the direction of the scattered, respectively, incident wave, are meant, is the unit sphere in M^. Let 
us use formulas (|5.1.19|l and (|5.1.20|l . 

It is of interest to obtain some information about q from the (fixed-energy) scattering data, that is, 
from the scattering amplitude A{a' , a), or, equivalently, from the coefficients A£{a). Very few results of 
such type are known. 

A result of such type is a necessary and sufficient condition for q{x) — <z(|a^|): it was proved ^ p. 131], 
that q{x) — q{\x\) if and only if A{a',a) = A{a' ■ a). Of course, the necessity of this condition was a 
common knowledge, but the sufficiency, that is, the implication: A{a',a) — A{a' ■ a) ^ q{x) = ^da;]), 
is a new result |R2| . 

A (modified) conjecture from p. 356] says that if the potential q{x) is compactly supported, and 
a > is its radius (defined for non-spherically symmetric potentials in the same way as for the spherically 
symmetric), then 



a — \ime- 



21 



V 



sup |^i'm(Q!)| 



-l<m<l 



lim,^co ( -\5i\* ] , (6.5.1) 



e 



where 5^ are the fixed-energy [k = 1) phase shifts. We prove (|6.5.1|l for the spherically symmetric 
potentials q = q{r). 

li q = q{r) then Afm(a) = agYimia) where Ui depends only on £ and fc, but not on a or a' . Since 
/c = 1 is fixed, ai depends only on £ for q ~ q{r). Assuming q — q{r), one takes A[a' , a) = A{a' ■ a) and 

calculates Aimia) = /02 A{a' ■ a)Ytm{a') da' — >£„(«), where ai :— — If — j\ A{t)ci^\t) dt, I — 

C*2'(i) 

0, 1, 2, . . . Here we have used formula (14.4.46) in |RKI p. 413], and cY^(t) are the Gegenbauer polyno- 
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mials (see |RKI p. 408]). Since c]^'' = Pt{t), Pi{i) — 1, where Pi{t) are the Legendre polynomials (see, 
e.g., EH P-409]), one gets: ai = 27r j\ A{t)Pt{t) dt. 

Formula (|6.5.1|l for q = q{r) can be written as a = lim^^oo (^^^r^l^A^^ ■ 

Indeed, sup ^^52 

= O (£2), as is well known (see, e.g., [MP, p. 261]). Thus lim^^oo 
-i<m<e ^ ' 



sup ^gg2 |l£m(Q;)| I — 1, and formula for (|6.5.1|l yields: 



2 / 1 \ 

a ^ -Ih^^^^oo (^la^^j • (6.5.2) 



e 

Note that assumption (A) implies the following assumption: 

Assumption (A'): the 'potential q{r) does not change sign in some left neighborhood of the point a. 

This assumption in practice is not restrictive, however, as shown in ^ p. 282], the potentials which 
oscillate infinitely often in a neighborhood of the right end of their support, may have some new properties 
which the potentials without this property do not have. For example, it is proved in jR| p. 282], that 
such infinitely oscillating potentials may have infinitely many purely imaginary resonances, while the 
potentials which do not change sign in a neighborhood of the right end of their support cannot have 
infinitely many purely imaginary resonances. Therefore it is of interest to find out if assumption A' is 
necessary for the validity of IjG. 5.2(1 . 

The main result is: 

Theorem 6.5.1. Let assumption (A) hold. Then formula (|6.5.2|l holds with lim replaced by lim. 
This result can be stated equivalcntly in terms of the fixed-energy phase shift df. 

lim (^l±l\Se\^] =a. (6.5.3) 

e-*oc \ e J 

Below, we prove an auxiliary result: 

Lemma 6.5.2. If q — q{r) G L^{0,oo), q{r) is real-valued and does not change sign in some interval 
(ai,a] where ai < a, and a is the radius of q, then 



a = lim,, 



q{r)r"^ dr 



,m = l,2,.... (6.5.4) 



Below we prove (|6.5.3|) and, therefore, (|6.5.1|l for spherically symmetric potentials. 

Proof of Lemma \6.5.Sl First, we obtain a slightly different result than ((6.5.4|l as an immediate conse- 
quence of the Paley- Wiener theorem. Namely, we prove Lemma 16.5.21 with a continuous parameter t 
replacing the integer m and lim replacing lim. This is done for q{r) G L^(0, a) and without additional 
assumptions about q. However, we are not able to prove Lemma l6.5.2l assuming only that q{r) G L^{0, a). 
Since q(r) is compactly supported, one can write 

/•oo pa p\na 

I{t):^ q{rydr= q{r)e*^'^'' dr = g(e'')e"e*"dM. (6.5.5) 



Let us recall that Paley- Wiener theorem implies the following claim (see ILevp : 

If f{z) = JI^^ g(u)e~^"^du, [61, ^'2] is the smallest interval containing the support of g{u), and g{u) G 
L^ibi,b2), then 

, , , \n\ f'' g{u)e'"du\ 

62 =limt^+oo (r4n|/(ii)|) = hmt^+oo ) -■ (6.5.6) 
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Thus, using (|6.5.5|l and H(j.5.6|l . one gets: 



In a 



(6.5.7) 



Formula (|6.5.7|l is similar to H6.5.4|l with m replaced by t and lim replaced by lim. 

Remark 6.5.3. We have used formula Hfi.5.6(l with bi = — oo, while in the Paley-Wiener theorem it is 
assumed that bi > — oo. However, for bi < 62, .9^0 on [62 — e, ^'2] for any e > 0, one has: 

g{u)e''^du= / g(u)e*"rfu+ / 5(M)e*"dw := /ii(t) + /i2(i)- 

-co J —00 Jhx 



Thus linif 



fei(t) 
'12(0 



0, and 



ln|/ii(i) + /i2W| 
hmt^oo 

— \n\h2{t)\ ln|l + o(l)| _ ln|;i2(t)| , 

=limt^oo ' + lim , = hmt^oo ' = In a. 

i t^oo t t 

Therefore formula (|6.5.7|l follows. 

To prove (|6.5.4|1 , we use a different approach independent of the Paley-Wiener theorem. We will use 
1)6.5.41) below, in formula (|6.5.19() . In this formula the role of q{r) in (|6.5.4|l is played by rq{r)[l + e(r, £)], 
where e — 0{j). Let us prove (|6.5.4)l . 

Assume without loss of generality that q > near a. Let / := q{r)r"^dr = j^' q{ry'^dr + 
Jai <l{'r)r''^dr := Ii + l2- We have < ca™, ci(a — 77)™ < I2 < C2a™, where is an arbitrary small 
positive number. Thus, / > for all sufficiently large m, and /i/'" — /2^™(1 + j^)^/™. One has 

X / 777' r 1 / 

a — rj < I2 < a and -jA ^ as m ^ 00. Since 77 is arbitrary small, it follows that limm^oc / '™ = a. 
This completes the proof of 16.5.4|) . Lemma [6. 5. 21 is proved. □ 

Proof of formula (|6.5.3|l . From (|5.1.19|l and (|5.1.23|) denoting := e*''*sini5^, one gets A[a' ■ a) = 



Y^'^=QaeYi{a)Yi{a') := 47r X^^o where, ai ■= k = 1, and . 



drug{r)q{r)'ijji{r), 



e'*' sin^£. 



(6.5.8) 



where ue{r) — rje{r) ^ sin (r — -y") as r 00, jiir) are the spherical Bessel functions, ji{r) :— 
^Je+i{r), and tpeir) solves (|5.1.15|) - (|5.1.17|l . and the integral 



where 



tpi{r) ^ ut{r) + I gt{r,s)q{s)t})t{s)ds, fc = 1, 
Jo 



gt{r,s) = -ui{r)wt{s), r < s; gi{r, s) = ge{s,r), 



(6.5.9) 
(6.5.10) 



wi{s) 



and nj,^'' is the Hankel function. 
It is known [HKl p. 407] that 



«W Y^f+y'^''' W£(r) = WyJ^+i(r), 



\2u/ J2-KV V TTZ^ V2z^ 



Mr)Hl'\r) L,^^+oo, 



(6.5.11) 



(6.5.12) 
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and |ARI Appendix 4]: 



< V 



16 



V > 



(6.5.13) 



It follows from H6.5.12I) that Ui{r) does not have zeros on any fixed interval (0, a] if £ is sufficiently large. 

" gi{r,s)ui{s) 



Define Vi{r) := Then ^^T^ yields 



vi{r) = 1 + 
From H6.5.1U|I and H6.5.12|l one gets 



ui{r) 



-q{s)vi{s)ds. 



Thus 



Ui{r) Vr/ 



£ +00. 



(6.5.14) 

(6.5.15) 
(6.5.16) 
(6.5.17) 



Ui{r) 2e+l' 

This implies that for sufficiently large £ equation H6.5.14|l has small kernel and therefore is uniquely 
solvable in C(0, a) and one has 



ipi{r) = Ui{r) 



l + 0(- 



+00, < r < a, 



(6.5.18) 



uniformly with respect to r e [0, a]. 

In the book formula (12.180), which gives the asymptotic behavior of Se for large £, is misleading: 
the remainder in this formula is of order which is much greater, in general, than the order of the main 
term in this formula. That is why we had to find a different approach, which yielded formula H6.5.18|l . 

From mir^ . ^(i. 5.1111 . ^(i.5.12^ . and (16.5.1 SH one has: 



Oi 



dr q{r)uf (r) 
dr q(r)r^r'^^ 



1 



Therefore, using IjG. 5.4(1 . one gets: 



lim 



2£ + l, 



1 27 = lim 



£^00 



M + 2 \2£+l 



2£+l 



(6.5.19) 



dr q(r)r r 



(6.5.20) 



Theorem 16 . 5 . II is proved. 



□ 



Remark 6.5.4. Since Si ^ as i ~f +00, and sa\5i ^ Se, e^^" ^ 1, as Se 0, formulas (|6.5.2U|I and 

ai — e^^' sinSg imply lim£_+oo (^ '^^^^ I'^^l^^ = 1, where 61 is the phase shift at a fixed positive energy. 
This is formula (|6.5.3|l . 
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Chapter 7 

Inverse scattering with "incomplete 
data" 

7.1 Uniqueness results 

Consider equation H1.2.3|) on the interval [0, 1] with boundary conditions u{0) = u{l) = 1 (or some 
other selfadjoint homogeneous separated boundary conditions), and q = q, q £ L^[0, 1]. Fix < 6 < 1. 
Assume q{x) on [b, 1] is known and a subset {Am(„)}vn=i.2,3.... of the eigenvalues A„ = A;^ of the operator 
£ corresponding to the chosen boundary conditions is known. Here 

171(71] 1 

^ =-(! + £„)> (T= const >0, |£„| < 1, £„->0. (7.1.1) 



n a 
We assume sometimes that 

oc 

Theorem 7.1.1. // H7.1.1|) holds uTid a > 26, then the data {q{x),h < x < 1; {^m{n)}vn} determine 
q{x) on [0,5] uniquely. If H7.1.1|) and (|7.1.2I) hold, the saTne conclusion holds also if <J = 2b. 

The number a is "the percentage" of the spectrum of £ which is sufficient to determine q on [0, b] if 
a > 2b and H7.1.2I) holds. For example, if cr = 1 and b = ^, then "one spectrum" determines q on the 
half-interval [0, ^J. If 6 = i, cr = i, then "half of the spectrum" determines q on [0, j]. Of course, q is 
assumed known on [b, 1]. If 6 = 1, ct = 2, then "two spectra" determines q on the whole interval. By 
"two spectra" one means the set {A„} U {fin}, where {^„} is the set of eigenvalues of £ corresponding to 
the same boundary condition u{0) = at one end, say at a; = 0, and some other selfadjoint boundary 
condition at the other end, say u'(l) = or u'{l) + hu{l) — 0, h = const > 0. The last result is 
a well-known theorem of Borg, which was strengthened in |^, where it is proved that not only the 
potential but the boundary conditions as well are uniquely determined by two spectra. A version of 
"one spectrum" result was mentioned in ILlj p. 81]. 

Proof of Theorem 7.1.1. First, assume a > 26. If there are qi and <?2 which produce the same data, then 
as above, one gets 

i-b b 
G{X) :— g{k) / p{x)ipi{x, k)^p2{x, k) dx — {ipiw' — (piw) ^ {ipiw' — ip[w) , (7.1.3) 







where w := ipi — ip2, p ■— qi ~ q2, k ~ \/X. Thus 



gik) = at fc = ±JX,nin) ±kn- (7.1.4) 
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The function G{X) is an entire function of A of order i (see (|1.2.11f) with k = VX), and is an entire 
even function of k of exponential type < 2b. One has 

g26|/mfe| 

The indicator of g is defined by the formula 

h{9) := hgi0) Th^ri^M:^!!}!, (7.1.6) 

r — >oo T 

where k = re^^ . Since \Imk\ = r \ sin0|, one gets from (|7.1.5() and (|7. 1.6(1 the following estimate 

h{0) <2b\sme\. (7.1.7) 
It is known [Levi formula (4.16)] that for any entire function g{k) ^ of exponential type one has: 

n(r) 1 f^'^ 

M^^Tnl ^^^'^''^ ^'-'-'^ 

where n{r) is the number of zeros of g{k) in the disk < r. From H?. 1.7(1 one gets 

1 f'^'^ Ih f'^^ Ah 

TT hg{0)de<— \sme\d9 = - (7.1.9) 
27r Jq 2n J„ tt 

From 1(7.1.2(1 and the known asymptotics of the Dirichlet eigenvalues: 

A„ = (vrn)^ + c + o(l), rt — > oo, c — const, (7.1.10) 

one gets for the number of zeros the estimate 



n{r)>2 1 = 2 — [l + o(l)], r -> oo. (7.1.11) 



From fTTHj) . fTT^I) and (|7.1.11|) it follows that 

cr < 26. (7.1.12) 

Therefore, if cr > 26, then g(fc) = 0. If g{k) = then, by property C^,, p{x) — 0. Theorem 17. 1.11 is 
proved in the case a > 26. 

Assume now that a — 2b and 

oo 

^|e„|<oo. (7.1.13) 

We claim that if an entire function G(A) in 1(7.1.3(1 of order ^ vanishes at the points 

9 9 

A„ = ^-(l + e„), (7.1.14) 

and ((7.1.13(1 holds, then G(A) = 0. If this is proved, then Theorem 17. 1.11 is proved as above. 
Let us prove the claim. Define 

<i'(A) n f 1 - (7.1.15) 



n=l 
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and recall that 



fJ-n := 



Since G(A„) = 0, the function 



w{X) 



GjX) 
ci>(A) 



(7.1.16) 



(7.1.17) 



is entire, of order < i. Let us use a Phragmen-Lindclof lemma. 

Lemma 7.1.2. lLeiA Theorem 1.22] If an entire function wiX) of order < 1 tias ttie property sup-oo<y<oo 
\'w{iy)\ < c, t/ien w(A) = c. //, m addition w{iy) — > as y —f +oo, tfien w(A) = 0. 

We use this lemma to prove that w{X) = 0. If this is proved then G(A) = and Theorem 17. 1.11 is 
proved. 

The function w{X) is entire of order 5 < 1. 
Let us check that 

sup \w{iy)\ < 00, (7.1.18) 

— cx;<y<oo 

and that 

\w{iy)\ ^ as y ^ +00. (7.1.19) 
One has, using (|7.1.5|l . (|7. 1.15(1 . 1(7.1.16(1 and taking into account that a = 2b: 



\w{iy)\ = 



G{iy) $o(«y) 



$(iy) $o(«y) 



< 



-1 1 + 4^' 



nr 



A2 



<- 



■i + |y| 



n f 



< 



— ^ n (i+k«i)<7 



Cl 



Here we have used elementary inequalities: 



l + "<^ if a>d>0- < 1 if < a < d, 



1 + d - d 



1 + d 



(7.1.20) 



with a := d := and the assumption ((7.1.13(1 . 
We also used the relation: 



sm{ay/vy) 



^a\Im,/Ty\ 

2(j\Vw\ 



as y ^ +00. 



Estimate ((7.1.20(1 implies ((7.1.18() and ((7.1.19() . An estimate similar to ((7.1.20(1 has been used in the 
literature (see GS,). 

Theorem 17 .1.1 1 is proved. □ 



7.2 Uniqueness results: compactly supported potentials 

Consider the inverse scattering problem of Section [5. 21 and assume 

q = for a; > a > 0. 



(7.2.1) 



Theorem 7.2.1. Ifq G satisfies ((7.2.1(1 . then any one of the data S{k), S{k), f{k), f'{k), determine 
q uniquely. 
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Proof. We prove first S{k) ^ q. Note that without assumption 1)7. 2. or an assumption which imphes 
that f{k) is an entire function on C, the result does not hold. If (|7.2.1|l holds, or even a weaker 
assumption: 

<cie-^^l"l\ 7>1, ci,C2>0, (7.2.2) 

then /(fc), the Jost function 11.2.5(1 . is an entire function of fc, and S(k) is a meromorphic function on C 
with the only poles in C+ at the points ikj, 1 < j < J. Thus, kj and J are determined by S{k). Using 
H1.2.1t)|l and H1.2.11|) . which holds for all k G C, because /(fc) and f{x, k) are entire functions of fc, one 
finds Sj =iReSk=tk, S{k). Thus all the data l|l.!^.17|) are found from S{k) if ifTTT)) (or (|7T^ ) holds. If 
the data H1.2.17|l are known, then q is uniquely determined, see Theorem 15. 2. II 

If 5{k) is given, then S{k) = e2**('=), so 5{k) ^ q. If /(fc) is given then S{k) = so /(fc) ^ q. If 

/'(O, fc) is given, then one can uniquely find f{k) from H1.2.11|l . Indeed, assume there are two /(fc), /i 
and /2, corresponding to the given /'(O, fc). Subtract from H1.2.11|l with f — fi equation 11.2.11|l with 
/ = /2, denote /i - := w, and get (*) f{0,k)w{-k) = f'{0,-k)w{k) or = jl^^-^- Since 

w(oo) = 0, and /'(0,fc) ~ ik — A{0,0) + Ax{0,y)e^''ydy, one can conclude that w = if one can 
check that is analytic in C+. The function /'(O, fc) has at most finitely many zeros in C+, and 

these zeros are simple. From {*) one concludes that if /'(0,k) = 0, k G C+, then w{k) = 0, because 
if f'{0,K) = then /'(0,-k) ^ (see (ll.lllll 'l. Thus is analytic in C+. Similarly -f^^ is 

analytic in C+. These two functions agree on the real axis, so, by analytic continuation, the function 
y^g''],^ is analytic in C+ and vanishes at infinity. Thus it vanishes identically. So w{k) = 0, fi — /2, and 
/(fc) is uniquely determined by /'(O, fc). Thus Theorem 17. 2. II is proved. □ 

7.3 Inverse scattering on the full line by a potential vanishing 
on a half-line 

The scattering problem on the full line consists of finding the solution to: 

lu-k'^u = 0, xeR, (7.3.1) 

M = e*'="+r(fc)e-^'="+o(l), a; ^ -oo, (7.3.2) 

u = t(fc)e*'''^ +o(l), x^+oo, (7.3.3) 

where r(fc) and t{k) are, respectively, the reflection and transmission coefficients. The above scattering 
problem describes plane wave scattering by a potential, the plane wave is incident from — oo in the positive 
direction of the a:-axis. The inverse scattering problem consists of finding q(x) given the scattering data 

{r(fc),fc„s„l<j < J}, (7.3.4) 

where Sj > are norming constants, kj > 0, and — fc| are the negative eigenvalues of the operator io- 

It is known |M], that the data H7.3.4|l determine q e Li.i(M) := {q : q = q, J^^{1 + \x\)\q\dx < oo} 
uniquely. Assume that 

q{x) =0, a; < 0. (7.3.5) 

Theorem 7.3.1. If q & -Li^i(M) and (|7.3.5|l holds, then {?'(fc)}vi;>o determines q uniquely. 

Proof. If (|7.3.5|l holds, then u ~ e'*^^ + r{k)e^^^'^ for a; < 0, and u = t{k)f{x, fc) for x > 0, where /(fc, x) 
is the Jost solution (|1.2.5|l . Thus 

ifc(l-r(fc)) u'(-0,fc) u'(+0,fc) /'(0,fc) 



l + r(fc) u(-0,fc) u(+0,fc) /(fc) 



/(fc). (7.3.6) 
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Therefore r{k) determines I{k), so by Theorem 13.1. 21 g is uniquely determined. □ 
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Chapter 8 

Recovery of quarkonium systems 



8.1 Statement of the inverse problem 

The problem discussed in this Section is: to what extent does the spectrum of a quarkonium system 
together with other experimental data determines the interquark potential? This problem was discussed 
in |TQR| , where one can find further references. The method given in |TQR| for solving this problem 
is this: one has few scattering data Ej, Sj, which will be defined precisely later, and one constructs, 
using the known results of inverse scattering theory, a Bargmann potential with the same scattering data 
and considers this a solution to the problem. This approach is wrong because the scattering theory is 
applicable to the potentials which tend to zero at infinity, while our confining potentials grow to infinity 
at infinity, and no Bargmann potential can approximate a confining potential on the whole semiaxis 
(0, oo). Our aim is to give an algorithm which is consistent and yields a solution to the above problem. 
The algorithm is based on the Gel'fand-Levitan procedure of Section IT^ 

Let us formulate the problem precisely. Consider the Schroedinger equation 

-VVj + = Ejijj in R^ (8.1.1) 

where q{r) is a real-valued spherically symmetric potential, r := |a;|, 2: G M'^, 

q{r) = r + p{r), p{r) — o(l) as r — > 00. (8.1.2) 

The functions ^/'j(x), ||V'j||L2(K3) = 1, are the bound states, Ej are the energies of these states. We define 
Uj{r) := rij)j(r)^ which correspond to s- waves, and consider the resulting equation for Uj: 

(uj -.^ -Uj + q(r)uj ^ EjUj, r > 0, Uj(0) = 0, ||uj||i2(o,oo) 1- (8.1.3) 

One can measure the energies Ej of the bound states and the quantities sj = u'j (0) experimentally. 
Therefore the following inverse problem (IP) is of interest: 
(IP): given: 

{£;„s,}v,=i,2,... (8.1.4) 

can one recover p(r)7 

In |TQR this question was considered but the approach in |TQR| is inconsistent and no exact results 
are obtained. The inconsistency of the approach in |TQR| is the following: on the one hand |TQR| 
uses the inverse scattering theory which is applicable to the potentials decaying sufficiently rapidly at 
infinity, on the other hand, |TQR| is concerned with potentials which grow to infinity as r ^ +00. It 
is nevertheless of some interest that numerical results in |TQR| seem to give some approximation of the 
potentials in a neighborhood of the origin. 

Here we present a rigorous approach to IP and prove the following result: 
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Theorem 8.1.1. IP has at most one solution and the potential q{r) can he reconstructed from data 
(|8. 1.4(1 algorithmically. 

The reconstruction algorithm is based on the Gel'fand-Levitan procedure for the reconstruction of 
q{x) from the spectral function. We show that the data 1(8.1.41) allow one to write the spectral function 
of the selfadjoint in L^(0, oo) operator i, defined by the differential expression ((8.1.3(1 and the boundary 
condition l(8.1.3() at zero. 

In Section l8 . 21 proofs are given and the recovery procedure is described. 

Since in experiments one has only finitely many data {Ej^ Sj}i<j<J^ fli^ question arises: 

How does one use these data for the recovery of the potential? 

We give the following recipe: the unknown confining potential is assumed to be of the form ((8.1.2(1 . and 
it is assumed that for j > J the data {Ej, Sj}j^j for this potential are the same as for the unperturbed 
potential qoir) = r. In this case an easy algorithm is given for finding q{r). 

This algorithm is described in Section! 



8.2 Proofs 

We prove Theorem 18 . 1 . 1 1 bv reducing (IP) to problem of recovery of q{r) from the spectral function. 

Let us recall that the selfadjoint operator L has discrete spectrum since q{r) +oo. The formula 
for the number of eigenvalues (energies of the bound states), not exceeding A, is known: 



Vl:-7V(A)^i/ [X-q{r)]idr. 

,^ ^ Jqir)<\ 



This formula yields, under the assumption q(r) ?' as r — > oo, the following asymptotics of the eigen- 
values: 

7^ /Stt 2 

E]^\-^3)^ asj^+oo. 
The spectral function p(A) of the operator L is defined by the formula 

-Bj<A ■> 

where aj are the normalizing constants: 

/>oo 

:= / <j)]{r)dr. (8.2.2) 
Jo 

Here 0j(r) :— (j}{r, Ej) and (j){r, E) is the unique solution of the problem: 

:= -(j)" + q{r)(j) ^ E^, r > 0, (/)(0, E) = 0, 0'(O, E) = 1. (8.2.3) 

If i<; = Ej, then (pj ~ (p{r,Ej) e i^(0, oo). The function (p{r, E) is the unique solution to the Volterra 
integral equation: 

Hr, E) . + /" ^^"[^^-^\ (,)0(,, E)dy. (8.2.4) 

VE Jq Vh 

For any fixed r the function (j) is an entire function of E of order i, that is, < cexp(c|£^|"'^/^), where 
c denotes various positive constants. At E = Ej, where Ej are the eigenvalues of l(8.1.3() . one has 
4){r,Ej) :— (pj £ _L^(0,oo). In fact, if q{r) ^ cr°' , a > 0, then \(f)j\ < cexp(— 7r) for some 7 > 0. 
Let us relate aj and Sj . From 1(8. 2. 3() with E = Ej and from l(8.1.3() , it follows that 

u 



.2.5) 
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Therefore 



1 



%:=ll^.lli=(o,oo) = ^- (8.2.6) 
Thus data (|8.1.4|l define uniquely the spectral function of the operator L by the formula: 



p(A) := ^ ,s^2. (8.2.7) 

Ej<\ 

Given p{X), one can use the Gel'fand-Levitan (GL) method for recovery of q{r). According to this 
method, define 

fT(A) :=p(A)-po(A), (8.2.8) 

where po(A) is the spectral function of the unperturbed problem, which in our case is the problem with 
q{r) = r, then set 

/oo 
0o(x,A)0o(2/,A)da(A), (8.2.9) 
-oo 

where 0o(a;, A) are the eigenfunctions of the problem H8.2.3(l with q{r) — r, and solve the second kind 
Fredholm integral equation for the kernel K{x,y): 

K{x, y)+ [ K{x, t)L(t, y)dt = -L(x, y), 0<y<x. (8.2.10) 

The kernel L{x,y) in equation (8.2.10) is given by formula (8.2.9). If K{x,y) solves (8.2.10), then 

p(r) = 2^^^, r>0. (8.2.11) 

ar 

8.3 Reconstruction method 

Let us describe the algorithm we propose for recovery of the function q{x) from few experimental data 
{Ej, Sj}i<j<j. Denote by {Ej , s'j}i<j<,j the data corresponding to qo := r. These data are known 
and the corresponding eigenfunctions H8.1.3|l can be expressed in terms of Airy function Ai(r), which 
solves the equation w" — rw = and decays at +oo, see |Leb| . The spectral function of the operator Lq 
corresponding to q = qo := r is 

Po(A) := Yl i^' f- (8-3.1) 

£;o<A 



Define 



p(A) :=po(A)+a(A), (8.3.2) 

Ej<\ E°<X 



and 



J J 

L{x,y) ■.= Y,s]<j>{x,E,)cjy{y,E,) - E(^?)''^^(^)'^^(y)' (^-3-4) 

where (j){x,E) can be obtained by solving the Volterra equation (|8.2.5|l with q{r) — qo{r) :— r and 
represented in the form: 

sin(£;l/2^) ,sin(^l/2y) , ^ 

cf,{x,E) = ^^^j^+ K{x,y)^^^dy, (8.3.5) 
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where K{x,y) is the transformation kernel corresponding to the potential q(r) — qa{r) :— r, and (f>j are 
the eigenf unctions of the unperturbed problem: 

-cf>] + r4>j ^ Ej4>j r > 0, 0^(0) = 0, 0^(0) = 1. (8.3.6) 

Note that for E ^ E^ the functions H8.3.5|) do not belong to 2.2(0, oo), but (j){0,E) = 0. We denoted 
in this section the eigcnfunctions of the unperturbed problem by (f>j rather than (j)Qj for simplicity of 
notations, since the eigcnfunctions of the perturbed problem are not used in this section. One has: 
(f>j{r) — CjAi{r — Ej), where Cj — [Ai' {— Ej)]^^ , Ej>Ois the j~th positive root if the equation 
Ai{—E) = and, by formula (|8.2.6|l . one has s° = [c| Ai'^{r — Ej)dr]~^^'^ . These formulas make the 
calculation of (pjix), Ej and easy since the tables of Airy functions are available |Lebj . 
The equation analogous to (|8.2.10|) is: 

K{x,y) + Y,Cj-^Ay) K{x,t)^,{t)dt^-Y,^j-^3{^)'^M^ (8-3.7) 

1 = 1 1 = 1 

where *j(t) := 0(i, Ej), cj = , 1 < j < J, and *j(t) = (pj-jit), cj = (s°_j)^, J + 1 < j < 2 J. Equation 
H8.3.7|l has degenerate kernel and therefore can be reduced to a linear algebraic system. 
If K{x,y) is found from (|8.3.7|l . then 

p{r) ^2-^K{r,r), q{r)=r+p{r). (8.3.8) 
dr 

Equation (|8.2.10|l and, in particular 18.3. 7() . is uniquely solvable by the Fredholm alternative: the homo- 
geneous version of 18.2.1(J|I has only the trivial solution. Indeed, \i h + Jp L{t,y)h{t)dt = 0,0 < y < x, 

then + f^^\h\^[dp{X) - po{X)] = 0, so that, by Parseval equality, |/ipo?p(A) = 0. Here 

h :— h{t)cj){t, X)dt, where cj){t, A) are defined by (|8.3.5|l . This implies that h{Ej) ~ for aU j = 1, 2, .... 
Since h[X) is an entire function of exponential type < x, and since the density of the sequence Ej is 
infinite, i.e., limA^oo ^^-x^ — oo, because Ej — 0{j^^^), as was shown in the beginning of Section HT^ it 
follows that h = and consequently h{t) — 0, as claimed. 

In conclusion consider the case when Ej — Ej,Sj = for all j > 1, and {Eo,sq} is the new 
eigenvalue, Eq < E^, with the corresponding data sq. In this case L{t, y) — SQ(po{t, EQ)(po{y, Eq), so that 
equation (|8.2.10() takes the form 

K{x,y) + sl(j)oiy) / K{x,t)(l)o{t,Eo)dt ^ -sl(j)oix,Eo)(l)oiy,Eo). 
Jo 

Thus, one gets: 

' drl + slj;<l^l{t,Eo)dt 
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Chapter 9 

Krein's method in inverse scattering 



9.1 Introduction and description of the method 

Consider inverse scattering problem studied in Chapter [3 and for simplicity assume that there are no 
hound states. This assumption is removed in Section 19.41 

This chapter is a commentary to Krein's paper |Klj. It contains not only a detailed proof of the 
results announced in |K1| but also a proof of the new results not mentioned in |K1| . In particular, it 
contains an analysis of the invertibility of the steps in the inversion procedure based on Krein's results, 
and a proof of the consistency of this procedure, that is, a proof of the fact that the reconstructed 
potential generates the scattering data from which it was reconstructed. A numerical scheme for solving 
inverse scattering problem, based on Krein's inversion method, is proposed, and its advantages compared 
with the Marchenko and Gel'fand-Levitan methods are discussed. Some of the results are stated in 
Theorem inn - Theorem I^Tsl below. 

Consider the equation for a function Tx{t, s): 



[I + Hx)Tx ■.^Tx[t,s) + H{t-u)Tx{u,s)du^ H{t- s), 0<t,s<x. 



(9.1.1) 







Equation H9.1.1|) shows that = (I + H^y^H = I - (I + H^)-'^, so 



(9.1.2) 



in operator form, and 



i?. = (/-r,)-i-/. 

Let us assume that H{t) is a real- valued even function 

H{-t) = H{t), H{t) e L^R)nL^{R), 



(9.1.3) 




(9.1.4) 



Then (|9.1.1|) is uniquely solvable for any x > 0, and there exists a limit 



r(t, s) = lim Tx{t, s 



):=roo(t,s), i,s>0. 



(9.1.5) 



where T{t, s) solves the equation 




< t, s < oo. 



(9.1.6) 
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Given H{t), one solves (|9.1.1|) . finds r2x(s,0), then defines 



where 



E{x, k) := e 



ikx 



Jo 



(9.1.8) 



Formula H9.1.8|l gives a one-to-one correspondence between E{x, k) and r22;(s,0). 

Remark 9.1.1. In 'Kl'l r2a:(0, s) is used in place of T2x{s,0) in the definition of E{x, k) . By formula 
<|9.2.22|l (see SectionW^ below) one has Tx{0,x) = Tx{x,0), but Tx{OjS) ^ ra;(s,0) in general. The 
theory presented below cannot be constructed with r2i:(0, s) in place ofT2x{s,0) in formula (9.1.8). 



Note that 
where 

and 

Furthermore, 



E{x,k)^e'''''f{-k)+o{l), x^+oo, (9.1.9) 

/(fc):=l-/ T{s)e'"'ds, (9.1.10) 
Jo 

r(s) := lim r,(s,0) :=roo(s,0). (9.1.11) 

X — *+oo 

^(x,fc) = e''=V(-fc)-e-^''V(fc) ^^(^^^ x^+oo. (9.1.12) 
2i 

Note that il^ix.k) = |/(fc)| sin(fca; 5(fc)) + o(l), x +oo, where /(fc) = \f{k)\e-'^^^\ 5{k) = 
-S{-k), keM. 

The function S{k) is called the phase shift. One has S{k) = e'^^^^'^K 

We have changed the notations from ^Kl^ in order to show the physical meaning of the function 
(|9.1.9|) : f{k) is the Jost function of the scattering theory. The function ^^^j,^^ is the solution to the 
scattering problem: it solves equation (1.2.3), and satisfies the correct boundary conditions: ^ 
and = e^-^e^') sin(fcx + d{k)) + o(l) as a; c». 

Krein |K1| calls S{k) := ^-jj^ the S'-function, and S{k) is the S'-matrix used in physics. 
Assuming no bound states, one can solve the inverse scattering problem (ISP): Given S{k) Vfc > 0, 
find q{x). 

A solution of the ISP, based on the results of IKlf . consists of four steps: 

1) Given S{k), find f{k) by solving the Riemann problem H9. 2.38(1 . 

2) Given f{k), calculate H(t) using the formula 

f°° 1 
,,H^,,J_Ji,)e'-,,^^^ (9X13) 

3) Given H{t), solve (|^T^ for r^{t, s) and then find r2x(2x, 0), < a; < cx). 

4) Define 

a(a:) =2r2,(2a;,0), (9.1.14) 

where 

a(0) = 2iJ(0), (9.1.15) 

and calculate the potential 

q{x) ^ a^{x) +a'{x), a(0) 2iJ(0). (9.1.16) 
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One can also calculate q{x) by the formula: 

= 2-^[r2,(2a;,0)-r2,(0,0)]. (9.1.17) 
ax 

Indeed, 2T2,{2x,0)^a{x), see (Pmil . Ar2,(0,0) = -2r2,(2x, 0)r2,(0, 2a:), see (^223, and r2,(2a;, 0) 
= r2^(0,2a;), see 

There is an alternative way, based on the Wiener-Levy theorem, to do step 1). Namely, given S{k), 
find S{k), the phase shift, then calculate the function g(t) :— 5{k) sm{kt)dk, and finally calculate 

/(fc)=exp(/;^5We^'=*&)- 

The potential q € Li^i generates the S'-matrix S'(fc), with which we started, provided that the 
following conditions - hold: 

S{k)=^(Pk) = S'\k), keR, (9.1.18) 

the overbar stands for complex conjugation, and 

indRS{k) = 0, (9.1.19) 

II^^(x)||l^(k+) + \\F{x)\\LiiM^) + ||xF'(x)|Ui(K^) < oo, (9.1.20) 

where 

1 

F{x):^— [1~ S{k)]e'''''dk. (9.1.21) 
27r J-oo 

By the index (|9. 1.19(1 one means the increment of the argument of S{k) ( when k runs from —oo to +oo 

along the real axis) divided by 2tt. The function 1(9.1.7(1 satisfies equation ((1.2.5(1 . Recall that we have 
assumed that there are no bound states. 

In Section 19.21 the above method is justified and the following theorems are proved: 

Theorem 9.1.2. // ((9.1.18(1 - (19.1.20(1 hold, then q{x) defined by ((9.1.16(1 is the unique solution to ISP 
and this q{x) has S{k) as the scattering matrix. 

Theorem 9.1.3. The function f{k), defined by (19.1. 10() . is the J ost function corresponding to potential 

(imm . 

Theorem 9.1.4. Condition ((9.1.4() implies that equation ((9.1.1() is solvable for all x > and its solution 
is unique. 

Theorem 9.1.5. // condition ((9.1.4(1 holds, then relation ((9.1.11() holds and T{s) :— roo(s,0) is the 
unique solution to the equation 



/ H{s-u)T{u)du^ H{s), s>0. (9.1.22) 
Jo 



T{s) 

The diagram explaining the inversion method for solving ISP, based on Krein's results, can be shown 



now: 



(9.2^4) (9.2.13) (9J^l) (trivial) n\(^-h}'^') / X (9-1,16) 1 OQ^ 

S[k) ^ f[k) ^ H[t) ^ Lx[t,s) ^ i22;(2a;,0) ^ a[x) ^ q[x). (9.1.23) 

SI S2 S3 S4 S5 SQ 

In this diagram §,„ denotes step number m. Steps S2, S4, S5 and sg ai'e trivial. Step si is almost trivial: 
it requires solving a Riemann problem with index zero and can be done analytically, in closed form. 
Step S3 is the basic (non-trivial) step which requires solving a family of Fredholm-type linear integral 
equations ((9.1.1() . These equations are uniquely solvable if assumption l(9.1.4() holds, or if assumptions 

mrm - mrm hold. 

In Section 19.21 we analyze the invertibility of the steps in diagram ((9.1.23() . Note also that, if one 
assumes ((9.1.18(1 - ((9.1.20(1 . diagram 1(9.1.23(1 can be used for solving the inverse problems of finding q{x) 
from the following data: 
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a) from f{k), Vfc > 0, 



b) from |/(fc)P, Vfc > 0, or 

c) from the spectral function dp{X). 

Indeed, if (|9.1.18|) - (|9.1.20|l hold, then a) and b) are contained in diagram H9.1.23|l . and c) follows 

{ V^ d\ A > n 

^ |/(Va)P' Let A = P. Then (stiU assuming 

0, A < 0. 

one has: dp = ^jjjj^dk, k > 0. 

Note that the general case of the inverse scattering problem on the half-axis, when indR5'(fc) := 7^ 0, 
can be reduced to the case = by the procedure described in Section 9.4, provided that S{k) is the 
S*— matrix corresponding to a potential q ^ Li i(M_|_). Necessary and sufficient conditions for this are 

conditions (Uma - ^TM- 

Section f9 . 31 contains a discussion of the numerical aspects of the inversion procedure based on Krein's 
method. There are advantages in using this procedure (as compared with the Gel'fand-Levitan pro- 
cedure): integral equation H9. 1.1(1 . solving of which constitutes the basic step in the Krein inversion 
method, is a Fredholm convolution-type equation. Solving such an equation numerically leads to inver- 
sion of Toeplitz matrices, which can be done efficiently and with much less computer time than solving 
the Gel'fand-Levitan equation 1)1.5.4(1 . Combining Krein's and Marchenko's inversion methods yields an 
eflacient way to solve inverse scattering problems. 

Indeed, for small x equation ((9.1.1(1 can be solved by iterations since the norm of the integral operator 
in 1(9.1.1(1 is less than 1 for sufficiently small x, say < x < xq. Thus q{x) can be calculated for < a: < ^ 
by diagram 1(9.1.23(1 . 

For x > xq > one can solve by iterations Marchenko's equation ((1.5.13(1 for the kernel A{x^y)^ 
where, if ((9.1.19(1 holds, the function F{x) is defined by the ((1.5.11() with Fd = 0. 

Indeed, for a; > the norm of the operator in ((1.5.11(1 is less than 1 and it tends to as a; -f-oo. 

Finally let us discuss the following question: in the justification of both the Gel'fand-Levitan and 
Marchenko methods, the eigenfunction expansion theorem and the Parseval relation play the fundamental 
role. In contrast, the Krein method apparently does not use the eigenfunction expansion theorem and the 
Parseval relation. However, implicitly, this method is also based on such relations. Namely, assumption 
((9.1.41) implies that the S'- matrix corresponding to the potential 1(9.1.16(1 . has index 0. If, in addition, this 
potential is in Li^i(R+), then conditions ((9.1.18(1 and 1(9.1.20(1 are satisfied as well, and the eigenfunction 
expansion theorem and Parseval's equality hold. Necessary and sufficient conditions, imposed directly on 
the function H{t)^ which guarantee that conditions ((9.1.18(1 - ((9.1.20() hold, are not known. However, it 
follows that conditions l(9.1.18() - ((9.1.20(1 hold if and only if H{t) is such that the diagram l(9.1.23() leads to 
a q{x) e Li^i(IR+). Alternatively, conditions ((9.1.18() - ((9.1.20() hold (and consequently, q{x) € Li^i(R+)) 
if and only if condition 1(9.1.4(1 holds and the function /(fc), which is uniquely defined as the solution to 
the Riemann problem 

$+(/=) = [l + i?(fc)]"^$_(/c), /ceR, (9.1.24) 

by the formula /(fc) = $-|_(/c), generates the S'-matrix S{k) by formula ((9.1.15() . and this S{k) satisfies 
conditions 1(9.1.18(1 - l(9.1.20() . Although the above conditions are verifiable, they are not quite satisfactory 
because they are implicit, they are not formulated in terms of structural properties of the function H{t) 
(such as smoothness, rate of decay, etc.). 

In Section Theorem 19.1.^ - Theorem 19.1.51 are proved. In Section numerical aspects of the 
inversion method based on Krein's results are discussed. In Section WM the ISP with bound states is 
discussed. In Section 1^31 a relation between Krein's and Gel'fand-Levitan's methods is explained. 



82 



9.2 Proofs 

Proof of Theorem \UT^ li v G L'^{0,x), then 



(w + H^v, v) = ^[(«, ■f)i2(R) + {Hv, v)L^fL)] (9.2.1) 
where the Parseval equahty was used, v :— u(s)e*'^*(is, 

{v,v) = f |wpds= (v,w)i2(R). (9.2.2) 
Jo 

Thus I + Hx is a positive definite selfadjoint operator in the Hilbert space L^{Q,x) if (|9. 1.4(1 holds. Note 
that, since H{t) E L^(R), one has H{k) — + as |fc| — > oo, so ((9.1.4(1 imphes 

1 + H{k)>c>0. (9.2.3) 

A positive definite selfadjoint operator in a Hilbert space is boundedly invertible. Theorem 19.1.41 is 
proved. □ 

Our argument shows that 

sup||(/ + ff,)-i|U2(R) <c-i. (9.2.4) 

2;>0 

Before we prove Theorem 19. 1.51 let us prove a simple lemma. For results of this type, see |K2j . 
Lemma 9.2.1. // (9.1.4) holds, then the operator 

/■oo 

Hif-.^l H{t-u)ip{u)du (9.2.5) 
Jq 

is a bounded operator in LP(M+), p — 1, 2, oo. 
For Fx (it, s) e L"'"(R_|_) one has 



II / duH{t - u)Tx{u,s)\\l2(^o,x) < ci du\T^{u,s)\. (9.2.6) 

J X J X 

Proof Let \\ip\\p := ||(p||lp(r+)- One has 

poo poo poo 

||if<^||i < sup / dt\H{t~u)\ Mu)\du< |i?(s)|ds||^||i =2||iJ||i||^||i, (9.2.7) 

ueR+ JO Jo J -oo 

where we have used the assumption H{t) = H{~t). Similarly, 

||i?<p||oo <2||il||i Halloo. (9.2.8) 
Finally, using Parseval's equality, one gets: 

2^\\H^\\l = ||i/^+||i2(R) < sup|il(fc)n|^||2, (9.2.9) 

where 

a; > 0, 



^,(x):=|-- (9.2.10) 

Since |-ff(fc)| < 2||iJ||i one gets from ((9.2.9(1 the estimate: 

\\H^\\2<^\\H\\,M\2. (9.2.11) 
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To prove 1)9. 2. one notes that 



dt\ duH{t - u)T^{u, s)\'' < snp / dt\H{t - u)H{t - v)\{ j \T ^.{u, s)\duY 

J X u,v>x Jo 



Estimate (|9. 2.6(1 is obtained. Lemma [9. 2. II is proved. 



<ci(/ du\T4u,s)\y. 



□ 



Proof of Theorem \9.1.5\ Define r^(t, s) = for t or s greater tlian x. Let w := Tx{t, s) — T{t, s). Then 
(|nTT|) and l|^TB|) imply 



{I + Ha:)w= / H{t-u)T{u,s)du:=hx{t,s). 



(9.2.12) 



If condition (|9.1.4|l holds, then equations H9.1.()l) and 19.1.22|l have solutions in i^(R+), and, since 
suptgji < oo, it is clear that this solution belongs to L°°(M_|_) and consequently to i^(R+), because 

ll'/'lb < II'/'IIooII'pIIi- The proof of Theorem 19 . 1 . 41 shows that such a solution is unique and does exist. 
From 1(9.2.4(1 one gets 

sup||(/ + i?,)-i||L2(o,.) <c-i. (9.2.13) 

2;>0 

For any fixed s > one sees that sup2.>j^ ||^a;(^, s)|| — > as y — > oo, where the norm here stands for any 
of the three norms LP{0,x),p = l,2,oo^ Therefore 1(9.2.121) and 1(9.2.11(1 imply 



\W\\t,2 



L2(0,x) 



< c 



< c 



-2 



H{t — u)T{u, s)du 



H(t — u)T{u, s)du 



< const ||r(M, s)||^i(^ ~> as a; oo, 

since r(w, s) e L\R+) for any fixed s > and H{t) G L\R). 
Also 



L°°{0,x) 



\\wit,s)\\U^o.x)<Mh,\\U^,^^^ + \\HM\Uio.x)) 

< Ci||r(u,s)||^i(^_^) +C2 SUp||i7(t- u)||^2(o,^)||w||i2(o,^). 



where Cj > are some constants. Finally, by 1(9. 2. 6() . one has; 



W^it^ s)\\L2io,x) < csi / |r(u, s)\duy ^ as a; . 



-oo. 



From 1(9.2.15(1 and 1(9.2.17(1 relation ((9.1.11(1 follows. Theorem 19.1.51 is proved. 

Let us now prove Theorem 19.1.31 We need several lemmas. 
Lemma 9.2.2. The function ((9.1.8(1 satisfies the equations 

E' ^ ikE - a{x)E^, S(0, fc) = l, E^ := E{x,-k), 



(9.2.14) 



(9.2.15) 
(9.2.16) 



(9.2.17) 



□ 



(9.2.18) 



E'_ ^ -ikE^ - a{x)E, E^{Q,k)^l, 



(9.2.19) 



where E' = ^ , and a{x) is defined in l(9.1.14() . 
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Proof. Differentiate (|9.1.8|l and get 

d{2x) 



E' = ikE - e''^" ( 2r2,(22;, 0)6-'^^"= + 2 / ^^fe^e-^^^'ds ) . (9.2.20) 





We wifl clieclc befow tfiat 

^I^^ = -T,{t,x)r,{x,s), (9.2.21) 

and 

r,{t,s)^r,{x-t,x-s). (9.2.22) 

Tlius, by 

= -r2,(5, 2x)r2,(2x, 0). (9.2.23) 

Tlierefore p. 2. 20(1 can be written as 

E' = ikE - e-'''='a{x) + a{x)e'''^ / r2x{s, 2x)e-''"'ds. (9.2.24) 

Jo 

By ^'^.'2.'}'4 one gets 

r2,(s,2a:) =r2,(2a;-s,0). (9.2.25) 

Tlius 

r2x /.2x 



/ r2x{s,2x)e-''"'ds^ / r2^(2a;-s,0)e*'=(^~")(is 



2x 



T2x{y,0)e"'ydy. (9.2.26) 



From H9.2.24|l and H9.2.26II one gets (|9.2.18|l . 

Equation (|9.2.19(l can be obtained from H9.2.18|l by changing k to —k. Lemma [9.2.21 is proved if 
formulas (|9.2.21|l - (|9.2.22(l are checked. 

To check H9.2.22|l . use H{—t) — H{t) and compare the equation for Tx{x — t,x — s) :— (p, 

r;:,(a; -t,x- s)+ I H{x - t - u)T^{u, x - s)du = H{x -t-x + s)= H{t - s), (9.2.27) 

with equation H9. 1.1(1 . Let u = x — y. Then (|9.2.27|l can be written as 

ip+ [ Hit~y)ipdy ^ H{t- s), (9.2.28) 
Jo 

which is equation ((9.1.1|l for (p. Since 1(9.1.1(1 has at most one solution, as we have proved above (Theo- 
rem EHI), formula 1(9.2.22(1 is proved. 

To prove 1(9.2.21(1 . differentiate ((9.1.1(1 with respect to x and get: 

or 

T'^{t,s)+ H{t-u)T'^{u,s)du^ -H[t- x)T^{x,s), := — ^. (9.2.29) 
JO "-^ 

Set s = a; in 1(9.1.1(1 . multiply 1(9.1.1(1 by —Tx{x,s), compare with 1(9.2.29(1 and use again the uniqueness 
of the solution to (pnTT|l . This yields ((9.2.21(1 . 

Lemma [9.2.21 is proved. □ 

Lemma 9.2.3. Equation \i.2.b\ holds for ip defined in ((9.1.7(1 . 
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Proof. From l|nT7|) and - one gets 

^„ ^ - ^ jtkE - a{x)E_y - i-ikE^ - a{x)Ey ^ ^^ ^.30) 

Using I9.2.18|l - l|9.2.19|l again one gets 

iP" = -k^ii) + q{x)%l), q{x):=a^{x) + a{x). (9.2.31) 

Lemma [9.2.31 is proved. □ 

Proof of Theorem \9.1.y[ The function ^ defined in (|9.1.7|l solves equation p. 2. 511 and satisfies the con- 
ditions 

V'(0,/c) = 0, iP'{0,k)^k. (9.2.32) 

The first condition is obvious (in "Kl" there is a misprint: it is written that '0(0, k) — 1), and the second 
condition follows from (|9.1.7|l and (|9.2.15|l : 



V/(0 k) - ~ - ^ "'^^ ~ ^^^^^ ~ 



2i 2i 



2ik _ ^ 
'2i ^ 



Let f{x, k) be the Jost solution. Since f{x, k) and f{x, —k) are linearly independent, one has ip — 
Ci/(x, k) + C2f{x, —k), where ci, C2 are some constants independent of x but depending on k. 

From 1)9.2.321) one gets ci — ^^^p-, C2 = ~2i^' ■^(^) /(O, fc)- Indeed, the choice of ci and 
C2 guarantees that the first condition (|9.2.32ll is obviously satisfied, while the second follows from the 
Wronskian formula: /'(O, k)f{-k) - f{k)f'{Q, ~k) = 2ik. 

Comparing this with H9.1.12|l yields the conclusion of Theorem 19. 1.31 □ 



Invertibility of the steps of the inversion procedure and proof of 
Theorem 1.1 

Let us start with a discussion of the inversion steps 1) - 4) described in the introduction. 

Then we discuss the uniqueness of the solution to ISP and the consistency of the inversion method, 
that is, the fact that q{x), reconstructed from S{k) by steps 1) - 4), generates the original S{k). 

Let us go through steps 1) - 4) of the reconstruction method and prove their invertibility. The 
consistency of the inversion method follows from the invertibility of the steps of the inversion method. 

Step 1. S{k) ^ /(fc). 

Assume S{k) satisfying (|9.1.18|) - l|9.1.2U|) is given. Then solve the Ricmann problem 

f{k) = S{-k)f{~k), fceK. (9.2.33) 

Since indKS'(fc) = 0, one has indRS'(— fc) = 0. Therefore the problem (|9.2.33() of finding an analytic 
function f+{k) in C+ :— {k : Imfc > 0}, /(fc) f+{k) in C+, (and an analytic function f-{k) f{~k) 
in C_ :— {k : Imfc < 0},) from equation 19.2.3311 can be solved in closed form. Namely, define 

/(fc)=exp|-L r l^i^^L^l lnrfc>0. (9.2.34) 



. 2vri y-k j 

Then /(fc) solves (l9.2.33ll . /+(fc) = /(fc), /_(fc) = /(-fc). Indeed, 

In /+(fc) - In /_(fc) = In S'(-fc), fceM (9.2.35) 
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by the known jump formula for the Cauchy integral. Integral 19.2.34|l converges absolutely at infinity, 
In S{—y) is difFerentiable with respect to y for y ^ 0, and is bounded on the real axis, so the Cauchy 
integral in (|9.2.34|) is well defined. 

To justify the above claims, one uses the known properties of the Jost function 

f{k) = 1 + / A{0,y)e'>'ydy 1 + / A{y)e'''ydy, (9.2.36) 



where estimates (ll.2.26|l and (ll.2.27|) hold and A{y) is a real-valued function. Thus 

^(0) 1 



f{k) = 1 - - 4 / A'{t)e'^'dt, (9.2.37) 
ik ik Jq 

S{-k)^4MT= ^'4"?'^-^^^ -l + ofrV (9-2-38) 



Therefore 



Also 



lnS'(-fc) = oQ^ as |fc|-^oo, keR. (9.2.39) 

f{k)=z A{y)ye^^ydy, / (9.2.40) 

Estimate H1.2.26fl implies 

y\A{y)\dy<2 t\q{t)\dt < cx,, A{y) e L^R+), (9.2.41) 



so that f{k) is bounded for all k&R, f(k) - 1 e £^(M), S'(-fc) is differentiable for fc 7^ 0, and In S{~y) 
is bounded on the real axis, as claimed. Note that 



f{-k) = /(fc), keR. (9.2.42) 

The converse step f{k) =i> S{k) is trivial: S{k) = ^-j^^- If indRS* = then f{k) is analytic in C+, 
/(fc) 7^ in C+, /(fc) = 1 + O (i) as |fc| ^ 00, fc e C+, and (|5.^.4^| holds. 

Step 2. /(fc) ^ H{t). 

This step is done by formula (|9.1.13|l : 

\ r°° /I 



One has H e L'^{R). Indeed, it follows from l|9.2.44|l that 

|/(fc)|2- 1 = -^^ A'(i)sin(fci)di + (^y^^ , IfcHoo, fceM. (9.2.44) 

The function 

w{k):=-i A' (t) sm{kt)dt (9.2.45) 
^ ^0 

is continuous because A'{t) e ^^1^+) by (|1.2.27(l . and w £ ^^(R) since w = o (^-^^ as |fc| — > 00, fc e R. 
Thus, H e L^{R). 
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Also, H G L"'^(]R). Indeed, integrating by parts, one gets from 19.2.431) the relation: 2TrH{t) = 
f iZc (^~'''*[fik)f{-k)-fi-k)f{k)]jjf^ := fg(i), and g S L^{R), therefore H S L\R). To check that 

g G L2(K), one uses i\):2.M . - ii:2:27\, . and i^.'lAul - 14111 . to conclude that [f{k)f{-k) - 

j{-k)f{k)] e i2(M), and, since f{k) 7^ on R and /(oo) = 1, it follows that g e L^^R). The inclusion 
[f{k)f{-k) - f{-k)f{k)] e follows from (19.2.361) . il. 2.261) - jl.2.2'7ll . and (l9.2.4()l) - (l9.2.4ll) . 

By (9.2.43), the function H'{t) is the Fourier transform of ~ik{l - |/(fc)|2)| /(fc)|-2, and, by (9.2.44), 

fc(|/(fc)|2 - 1) = -2 A'(i) sm{kt)dt + O (^-^^ , as |fc| ^00, keR. Thus, i/'(i) behaves, essentially, 

as A'{t) plus a function, whose Fourier transform is O (^jj^^ ■ Estimate (1.2.27) shows how A'(t) behaves. 
Equation (9.1.1) shows that rx{t,0) is as smooth as H{t), so that formula (9.1.17) for q{x) shows that 
q is essentially is as smooth as A'{t). 
The converse step 

H{t) ^ f{k) (9.2.46) 

is also done by formula H9.1.13|) : Fourier inversion gives |/(fc)P ~ f{k)f{—k), and factorization yields 
the unique f{k), since f{k) does not vanish in C+ and tends to 1 at infinity. 

Step 3. H ^r^{s,0) ^T2x{2x,0). 

This step is done by solving equation (|9.1.1|l . By Theorem 19 . 1 . 41 eg nation 19.1.1|) is uniquely solvable 
since condition (|9.1.4|l is assumed. Formula H9.1.13|) holds and the known properties of the Jost function 
are used: /(fc) ^ 1 as fc ^ ±00, /(fc) 7^ for fc 7^ 0, fc e R, /(O) ^ since indES'(fc) = 0. 
The converse step ra;(s,0) H{t) is done by formula H9.1.3|) . The converse step 

r2,(2a;,0)^r,(s,0) (9.2.47) 

constitutes the essence of the inversion method. 
This step is done as follows: 

i 2x[2x, 0) ^ a[x) ^ t^yx, fc) ^ 1 2,(s, 0). (9.2.48) 

Given a{x), system H9.2.18() - H9.2.19|l is uniquely solvable for E{x,k). 

Note that the step q{x) ^ /(fc) can be done by solving the uniquely solvable integral equation (|1.2.6|l : 
with q e Li,i(R+), and then calculating f{k) ~ f{Q,k). 

Step 4. a{x) := 2T2xi2x, 0) ^ q{x). 

This step is done by formula (|9.1.16|l . The converse step 

q{x) =5> a{x) 

can be done by solving the Riccati problem H9.1.16|l for a{x) given q{x) and the initial condition 2H{0). 
Given q{x), one can find 2H{0) as follows: one finds f{x, k) by solving equation H1.2.6() . which is uniquely 
solvable if q € Li_i(R+), then one gets f{k) := /(O, k), and then calculates 2iJ(0) using formula H9.2.43|) 
with t = 0: 

Proof of Theorem \9.1.iA If (|9.1.18|l - H9.1.2U|I hold, then, as has been proved in Section [^31 there is a 
unique q{x) G ii^i(R+) which generates the given 5-matrix S{k). 

It is not proved in ]K1^ that q{x) defined in (1-19) (and obtained as a final result of steps 1) -4)) 
generates the scattering matrix S{k) with which we started the inversion. 

Let us now prove this. We have already discussed the following diagram: 

Sik) ^ f{k) ^ Hit) ^ T4s, 0) ^ T242X, 0) ^ a{x) ^ q{x). (9.2.49) 
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To close this diagram and therefore establish the basic one-to-one correspondence S{k) q{x), one 
needs to prove r2a;(2a;,0) =J> ra;(s,0). This is done by the scheme (19.2. 48(1 . 

Note that the step q{x) a{x) requires solving Riccati equation (|9.1.16|l with the boundary condition 
a(0) — 2H{0). Existence of the solution to this problem on all of M+ is guaranteed by the assumptions 
(|9.1.18|l - (|9.1.20|) . The fact that these assumptions imply q{x) G Li.i(R+) is proved in Section [5.51 
Theorem 19. 1.21 is proved. □ 

Uniqueness theorems for the inverse scattering problem are not given in |K1| . They can be found in 
Section [^31 

Remark 9.2.4. From our analysis one gets the following result: 

Proposition 9.2.5. If q{x) G _Li.i(K+) and has no hounds states and no resonance at zero, then Riccati 
equation ()9.1.16|l with the initial condition (|9.1.15|) has the solution a{x) defined for all x G M-|_. 

9.3 Numerical aspects of the Krein inversion procedure. 

The main step in this procedure from the numerical viewpoint is to solve equation (|9. 1.1(1 for all x > 
and all < s < a;, which are the parameters in equation (|9.1.1|l . 

Since equation 1)9.1.1(1 is an equation with the convolution kernel, its numerical solution involves 
inversion of a Toeplitz matrix, which is a well developed area of numerical analysis. Moreover, such an 
inversion requires much less computer memory and time than the inversion based on the Gel'fand-Levitan 
or Marchenko methods. This is the main advantage of Krein's inversion method. 

This method may become even more attractive if it is combined with the Marchenko method. In the 
Marchenko method the equation to be solved is ((1.5.13(1 where F{x) is defined in ((1.5.11(1 and is known if 
S(k) is known. The kernel A{x, y) is to be found from ((1.5.11() and if A{x^ y) is found then the potential 
is recovered by the formula: Equation ((1.5. 11() can be written in operator form: (/ -I- Fx)A = —F. The 
operator F^ is a contraction mapping in the Banach space L^{x, oo) for x > 0. The operator in ((9.1.1(1 
is a contraction mapping in L°°(0, x) for < x < xq, where xq is chosen to that " \H{t ~ u)\du < 1. 
Therefore it seems reasonable from the numerical point of view to use the following approach: 

1. Given S{k), calculate f{k) and H{t) as explained in Steps 1 and 2, and also F{x) by formula 1(1.5.11(1 . 

2. Solve by iterations equation l(9.1.1() for < x < xq, where xq is chosen so that the iteration method 

for solving ((9.1.6(1 converges rapidly. Then find q(x) as explained in Step 4. 

3. Solve equation (1.5.13) for x > xq hy iterations. Find q{x) for x > xq hy formula (1.5.12). 

9.4 Discussion of the ISP when the bound states are present. 

If the given data are l(9.1.15() . then one defines w{k) — Y[j=i k+ik- hidR5(x) = — 2J and W{k) — 
k^i^ 'w{k) if mdfjS{k) = — 2J — 1, where 7 > is arbitrary, and is chosen so that 7 ^ fcj, 1 < j < J. 

Then one defines Si{k) := S{k)w^{k) if indR5 = -2J or Si{k) := S{k)W^{k) if indR5 = 
-2 J - 1. Since ind^w'^lk) = 2 J and indjuW^ik) = 2 J -I- 1, one has indRS'i(fc) = 0. The theory of 
Section applies to Si{k) and yields qi{x). From qi{x) one gets q{x) by adding bound states — fcj and 
norming constants Sj using the known procedure (e.g. see 

9.5 Relation between Krein's and GL's methods. 

The GL (Gel'fand-Levitan) method in the case of absence of bound states of the following steps (see 
Chapter 21 for example): 
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Step 1. Given f{k), the Jost function, find 

2 n ,2 f 1 i\ sinfca: sinky 



2 f 

L{x, y) ■— — dkk 
Jo 



dk {\f{k)\-^ - 1) {cos[k{x - y)] - cos[k{x + y)]) 



_ 1 

TT 

:= M{x -y)- M{x + y) 



where M{x) i dk {\f{k)\-^ - l) cos(fca:). 

Step 2. Solve the integral equation (|1.5.4|l for K{x^ y): 

Step 3. Find q{x) = 2^^^^f^. Krein's function, H{t), see can be written as follows: 

Hit) = 7r (|/(fc)r2-l)e-^'=*& = - / {\f(k)r^~l)cos{kt)dk. (9.5.1) 

Thus, the relation between the two methods is given by the formula: 

M{x)^H{x). (9.5.2) 

,2 

In fact, the GL method deals with the inversion of the spectral foundation dp of the operator —-^ + 
q{x) defined in L^(R+) by the Dirichlet boundary condition at a; = 0. However, if indRS'(fc) = 
(in this case there are no bound states and no resonance at fc = 0), then (see (|1.2.21|l '): dp{X) = 
( 2k^dk A > A — A:^ 

; WUW' ^ ^ 'so dp{X) in this case is uniquely defined by f{k), k >0. 

I 0, A < 0, 



90 



Chapter 10 

Inverse problems for the heat and 
wave equations. 



10.1 Inverse problem for the heat equation 

Consider problem (|1.5.25|l - (|1.5.28|l . Assume 

a{t) = for or, / a{t)dt < oo, a{t) ^ 0. (10.1.1) 
Jo 

One can also take a{t) — d{t) where d{t) is the delta-function. We prove that the inverse problem of 
finding q{x) G L^[0, 1], q — q, from the conditions l|1.5.25|l - H1.5.28(l has at most one solution. 
If (|1.5.28|l is replaced by the condition 

u,{0,t) = bo{t), (10.1.2) 

then q{x), in general, is not uniquely defined by the conditions (ll.5.25|) . (|1.5.2()|l . (|1.5.27|) and (|10.1.2|l . 
but q is uniquely defined by these data if for example, q{^ — x) ~ q{^ + x), or if q{x) is known on 1] . 

Let us take the Laplace transform of H1.5.25|l - (|1.5.28l) and put v[x, A) :— u{x, t)e~^^dt, A{\) :— 
w(l, A), B{X) ■= Va:{l, A), Ba{X) := w^(0, A). Then H1.5.25|) - H1.5.28|l can be written as 

ev + Xv := ~v" + q{x)v + Xv ^ 0, < x < 1, t;(0,A)=0, u(l,A)=A(A) (10.1.3) 



v'{l,X)^B{X) (10.1.4) 

and (|10.1.2|l takes the form 

u'(0, A) Bo(A). (10.1.5) 

Theorem 10.1.1. The data {A(X), B{X)}, known on a set o/ A G (0, oo), which has a finite positive 
limit point, determine q uniquely. 

Proof. Since A{X) and B{X) are analytic in Y[+ ■= {A : 5RA > 0}, one can assume that A{X) and B{X) 

are known for all A > 0. If fc = iX^ and (p is defined in H1.2.3|l then v{x,X) = c{k)ip{x, k), c{k) ^ 0, 
A{X) = c{k)ip{l, k), B{X) = c{k)tp'{l, fc), so 

B{X) p'{l,k) 
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Thus the function is meromorphic in C, its zeros on the axis fc > are the eigenvalues of ^ = 

~lx^ + q{x), corresponding to the boundary conditions u{0) = — and its poles on the axis fc > 
are the eigenvalues of £ corresponding to u(0) = = 0. The knowledge of two spectra determines q 
uniquely (Section [OJ- 1^ 

An alternative proof of Theorem \10.1.1[ based on property Ctp, is: assume that qi and q2 generate 
the same data, p := qi — (72, w := vi — V2, where vj, j = 1, 2, solves - H1Q.1.4|I with q = q, and 

get (*) liw = pv2, w{0, A) = w{l, A) = A) = 0. Multiply (*) by tpi, inpi + \(pi = 0, v?i(0, A) = 0, 
93'(i(0, A) — 1, and integrate over [0, 1] to get 

I pv2^idx = VA > 0. (10.1.7) 
Jo 

By property Cip it follows from H10.1.7() that p = 0. Theorem 110. 1.11 is proved. □ 

Theorem 10.1.2. Data (I10.1.3() . (|10. 1.5(1 does not determine q uniquely in general. It does if q(x) is 
known on [i, 1], or if q{x + ^) = q{^ — x). 

Proof. Arguing as in the first proof of Theorem 110. 1.11 one concluded that the data (|10.1.3|l . (|10.1.5() 
yields only one (Dirichlet) spectrum of ^, since (/^'(O, fc) — 1. One spectrum determines q only on "a half 
of the interval", b = ^, see Section ITTI Theorem 1 1 . 1 . 21 is proved. □ 



10.2 What are the "correct" measurements? 

From Theorem llO.l.ll and Thcorcm llO. 1 .2l it follows that the measurements {mx(1, t)}vt>o are much more 
informative than {ux{0,t)}\/t>o for the problem (|1.5.25|) - (|1.5.27|) . In this section we state a similar 
result for the problem 

ut = (a(x)u')', < a; < 1, t > 0; u{x, 0) = 0, u(0, t) = 0, (10.2.1) 



u(l,t) = /(t). (10.2.2) 

The extra data, that is, measurements, are 

a(l)u'(l,i) = .g(t), (10.2.3) 

which is the fiux. Assume: 

/^O, f£L\0,l), a{x) eW'^^'^iO,!), a{x) > c> 0, (10.2.4) 

W^'P is the Sobolev space. Physically, a{x) is the conductivity, u is the temperature. We also consider 
in place of I|1U.2.3|I the following data: 

aiO)u{0,t) = h{t). (10.2.5) 

Our results are similar to those in Section llU.il 

data {f{t), g{t)}\/t>o determine q{x) uniquely, while data {/(i), ft-(i)}vt>o do not, in general, deter- 
mine a{x) uniquely. 

Therefore, the measurements {<?(i)}vt>o are much more informative than the measurements {/i(t)}vt>o- 
We refer the reader to |R9 |. 
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10.3 Inverse problem for the wave equation 

Consider inverse problem (|1.5.20() - H1.5.24() . Our result is 

Theorem 10.3.1. The above inverse problem has at most one solution. 

Proof. Take the Fourier transform of H1.5.20|l - (|1.5.24|l and get: 

/•OO 

fo-fc^w^O, a;>0, v{x,k):^ e'''\{x,t)dt, (10.3.1) 

Jo 

POO 

v{0,k) = 1, ^(l,^:) = A{k) = / a{t)e'^^dt. (10.3.2) 

Jo 

From 11U.3.1II one gets v{x,k) = c{k)f{x,k), where f{x,k) is the Jost solution, and from (|10.3.2|) one 
gets v{x, k) — ^jjj^ and A{k) = ^jjj^ — J(k)^ because g = for a; > 1. Thus /(fc) = -^jj^ is known. By 
Theorem lT. 2.11 a is uniquely determined. Theorem 111). 3. II is proved. □ 

Remark 10.3.2. The above method allows one to consider other boundary conditions at x = 0, such as 
u'(0, t) — or u'(0, t) = hu{0, t), h = const > 0, and different data at x = 1, for example, u'(l, t) — b{t). 
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Chapter 11 

Inverse problem for an 
inhomogeneous Schrodinger 
equation 



In this chapter an inverse problem is studied for an inhomogeneous Schrodinger equation. Most, if not 
all, of the earlier studies dealt with inverse problems for homogeneous equations. Let 

du 

tu- k'^u := -u" + q(x)u- k'^u ^ 5(x), a: G M\ ——- iku ^ 0, b| oo. (11.1.1) 

Assume that q{x) is a real- valued function, q{x) = for \x\ > 1, q E 1, 1]. Suppose that the data 

1, fc), u(l, fc)}, Vfc > are given. 
The inverse problem is: 
(IP) Given the data, find q{x). 

This problem is of practical interest: think about finding the properties of an inhomogeneous slab 
(the governing equation is plasma equation) from the boundary measurements of the field, generated 
by a point source inside the slab. Assume that the self-adjoint operator i — —-j^ + q{x) in L^(R) has 
no negative eigenvalues (this is the case when q{x) > 0, for example). The operator i is the closure in 
L^(M) of the symmetric operator £o defined on C^(M"'^) by the formula Iqu = —u" + q{x)u. Our result 
is: 

Theorem 11.1.3. Under the above assumptions IP has at most one solution. 
Proof of Theorem 11.1: The solution to (|11.1.1|) is 

f sp.,f{x,k), x>Q, , ^ 

Here /(x, k) and g{x, k) solve homogeneous version of equation Hll.l.l|l and have the following asymp- 
totics: 

/(x,A:) - 6*'==", a;->-|-oo, g{x,k) ^ e-''''' , x ^ -oo, (11.1.3) 
/(fc) :=/(0,A;), g{k) := g{0,k), (11.1.4) 

[f.g] ■■=fg' -f'g = -2ika{k), (11.1.5) 
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where the prime denotes differentiation with respect to x-variable, and a{k) is defined by the equation 

f{x, k) = b{k)g{x, k) + a{k)g{x, -k). (11.1.6) 

It is known that 

g{x, k) ^ -b{-k)f{x, k) + a{k)f{x, -k), (11.1.7) 



ai-k) = a{k), h{-'k)^h{k), |a(fc)|'' = 1 + |6(A:)|^ /c e M, (11.1.8) 
a(/c) = l + 0(i), fc^oo, /ceC+; 6(fc) = 0(i), |fc| ^ oo, fc e M, (11.1.9) 

[f{x,k),g{x,~k)] = 2ikb{k), [f{x,k),g{x,k)]=~2ika{k), (11.1.10) 

a{k) in analytic in C+, b{k) in general does not admit analytic continuation from K, but if q{x) is 
compactly supported, then a{k) and b{k) are analytic functions of fc € C \ 0. 
The functions 

.l,(/c) ,^ 9ik)fil,k) ^ ^^(fc) := /Wffhk^ (11.1.11) 
—2ika{k) ' ~2ika{k) 

are the data, they are known for all fc > 0. Therefore one can assume the functions 

hiik) -.^ h,{k):=l§ (11.1.12) 

a[k) a[k) 

to be known for all fc > because 

/(l,fc) = e^^ g{-l,k) = e'^ (11.1.13) 

as follows from the assumption q = if |x| > 1, and from (|11. 1.3(1 . 
From ^ll.l.ljjp . (|11.1.7| and (|11.1.6|l it follows that 

a{k)hi{k) = ~b{~k)f{k) + a{k)f{-k) = -b{~k)h2{k)a{k) + h2{-k)a{-k)a{k), (11.1.14) 

a(fc)/i2(fc) =^ b{k)a{k)hi{k) + a{k)hi{~k)a{~k). (11.1.15) 
From H11.1.14|l and (|11.1.15|) it follows: 

-6(-fc)/i2(fc) + /i2(-fc)a(-fc) = /ii(fc), (11.1.16) 

b{k)hi{k) + a(-fc)/ii(-fc) = h2{k). (11.1.17) 
Eliminating b{-k) from (|11.1.16|) and (|11.1.17|) . one gets: 

a{k)hi{k)h2{k) + a{-k)hi{-k)h2{-k) ^ hi{k)hi{-k) + h2{-k)h2{k), (11.1.18) 

or 

a(fc) = m(fc)a(-fc) + n(fc), fc G R (11.1.19) 

where 

Problem (|11. 1.19(1 is a Riemann problem for the pair {a(fc), a(— fc)}, the function a{k) is analytic in 
C+ := {fc : fc € C, Imk > 0} and a{—k) is analytic in C_. The functions a(fc) and a(— fc) tend to one as 
fc tends to infinity in C+ and, respectively, in C_ , see equation ((11.1.9(1 . 
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The function a{k) has finitely many simple zeros at the points ikj, ^ ^ j J , kj > 0, where — A:| 
are the negative eigenvalues of the operator £ defined by the differential expression £u = —u" + q{x)u in 

The zeros ikj are the only zeros of a(fc) in the upper half-plane k. 
Define 

inda{k):= / d\na{k). (11.1.21) 

One has 

inda^J, (11.1.22) 

where J is the number of negative eigenvalues of the operator £, and, using H11.1.12|l . H11.1.22|l and 
<|11.1.2U|I . one gets 

indm{k) = ~-2[indhi{k) + indh2{k)] = -2[ind g{k) + ind f{k) - 2J]. (11.1.23) 

Since £ has no negative eigenvalues, it follows that J = 0. 

In this case indf{k) — indg(k) — (see Lemma 1 below), so indm{k) — 0, and a{k) is uniquely 
recovered from the data as the solution of H11.1.19|l which tends to one at infinity, see equation H11.1.9|l . 
If a{k) is found, then h{k) is uniquely determined by equation (|11.1.17|) and so the reflection coefficient 
r{k) :— is found. The reflection coefficient determines a compactly supported q{x) uniquely |R9| . 
but we give a new proof. If qix) is compactly supported, then the reffection coefficient r(fc) : — ^^^j 
is meromorphic. Therefore, its values for all fc > determine uniquely r(k) in the whole complex k- 
plane as a meromorphic function. The poles of this function in the upper half-plane are the numbers 
ikj, 3 = 1, 2, J. They determine uniquely the numbers kj, ^ ^ j ^ J, which are a part of the standard 
scattering data {r(k),kj, Sj,l < j < J}, where Sj are the norming constants. Note that if a{ikj) — 
then b{ikj) ^ 0: otherwise equation Hll.l.t)|) would imply f{x,ikj) = in contradiction to the first 
relation (|11.1.3|l . If r(fc) is meromorphic, then the norming constants can be calculated by the formula 
Sj = = —iResk=ikjr{k), where the dot denotes differentiation with respect to fc, and Res denotes 

the residue. So, for compactly supported potential the values of r(fc) for all fc > determine uniquely 
the standard scattering data, that is, the reflection coefficient, the bound state s — and the norming 
constants Sj, 1 < j < J. These data determine the potential uniquely. Theorem 1 1 1 . 1 . 31 is proved. □ 

Lemma 11.1.4. // J = then indf = indg = 0. 

Proof. We prove indf = 0. The proof of the equation indg = is similar. Since indf(k) equals to the 
number of zeros of /(fc) in C+, we have to prove that /(fc) does not vanish in C-|-. If f{z) = 0, 2; € C+, 
then z — ik, k > 0, and — fc^ is an eigenvalue of the operator £ in L^(0,cx)) with the boundary condition 
u(0) = 0. 

From the variational principle one can find the negative eigenvalues of the operator £ in L^(M+) 
with the Dirichlet condition at a; = as consequitive minima of the quadratic functional. The minimal 
eigenvalue is: 

o 

-k^=inf [u'^ + qix)u^] dx := Ko, ueH\R+), \\u\\l2^m^) = 1, (11.1.24) 
Jo 

o 

where ff"'^(K+) is the Sobolev space of _ff^(]R+)-bmctions satisfying the condition u{0) ~ 0. 
On the other hand, if J = 0, then 

/•OO 

0<inf [u^ +q{x)u^]dx := Ki, ueH^{R), ||u||i2(R) = 1. (11.1.25) 



Since any element u of 7f ^(R+) can be considered as an element of iJ^(R) if one extends u to the whole 
axis by setting u = for a; < 0, it follows from the variational definitions H11.1.24|l and H11.1.25|l that 
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Ki < kq- Therefore, if J = 0, then ki > and therefore kq > 0. This means that the operator £ on 
L^(M+) with the Dirichlet condition at x = has no negative eigenvalues. This means that /(fc) does 
not have zeros in C+, if J = 0. Thus J = imphes indf{k) = 0. 

Lemma Fl 1.1. 41 is proved. □ 

Remark 11.1.5. The above argument shows that in general 

ind f < J and indg < J, (11.1.26) 

so that (|11.1.23| ) implies 

indm{k)>0. (11.1.27) 
Therefore the Riemann problem (|11. 1.19(1 is always solvable. 
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